cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 26 results. Next

A003462 a(n) = (3^n - 1)/2.

Original entry on oeis.org

0, 1, 4, 13, 40, 121, 364, 1093, 3280, 9841, 29524, 88573, 265720, 797161, 2391484, 7174453, 21523360, 64570081, 193710244, 581130733, 1743392200, 5230176601, 15690529804, 47071589413, 141214768240, 423644304721, 1270932914164
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A000244. Values of base 3 strings of 1's.
a(n) = (3^n-1)/2 is also the number of different nonparallel lines determined by pair of vertices in the n dimensional hypercube. Example: when n = 2 the square has 4 vertices and then the relevant lines are: x = 0, y = 0, x = 1, y = 1, y = x, y = 1-x and when we identify parallel lines only 4 remain: x = 0, y = 0, y = x, y = 1 - x so a(2) = 4. - Noam Katz (noamkj(AT)hotmail.com), Feb 11 2001
Also number of 3-block bicoverings of an n-set (if offset is 1, cf. A059443). - Vladeta Jovovic, Feb 14 2001
3^a(n) is the highest power of 3 dividing (3^n)!. - Benoit Cloitre, Feb 04 2002
Apart from the a(0) and a(1) terms, maximum number of coins among which a lighter or heavier counterfeit coin can be identified (but not necessarily labeled as heavier or lighter) by n weighings. - Tom Verhoeff, Jun 22 2002, updated Mar 23 2017
n such that A001764(n) is not divisible by 3. - Benoit Cloitre, Jan 14 2003
Consider the mapping f(a/b) = (a + 2b)/(2a + b). Taking a = 1, b = 2 to start with and carrying out this mapping repeatedly on each new (reduced) rational number gives the sequence 1/2, 4/5, 13/14, 40/41, ... converging to 1. Sequence contains the numerators = (3^n-1)/2. The same mapping for N, i.e., f(a/b) = (a + Nb)/(a+b) gives fractions converging to N^(1/2). - Amarnath Murthy, Mar 22 2003
Binomial transform of A000079 (with leading zero). - Paul Barry, Apr 11 2003
With leading zero, inverse binomial transform of A006095. - Paul Barry, Aug 19 2003
Number of walks of length 2*n + 2 in the path graph P_5 from one end to the other one. Example: a(2) = 4 because in the path ABCDE we have ABABCDE, ABCBCDE, ABCDCDE and ABCDEDE. - Emeric Deutsch, Apr 02 2004
The number of triangles of all sizes (not counting holes) in Sierpiński's triangle after n inscriptions. - Lee Reeves (leereeves(AT)fastmail.fm), May 10 2004
Number of (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < 6 and |s(i) - s(i-1)| = 1 for i = 1, 2, ..., 2*n + 1, s(0) = 1, s(2n+1) = 4. - Herbert Kociemba, Jun 10 2004
Number of non-degenerate right-angled incongruent integer-edged Heron triangles whose circumdiameter is the product of n distinct primes of shape 4k + 1. - Alex Fink and R. K. Guy, Aug 18 2005
Also numerator of the sum of the reciprocals of the first n powers of 3, with A000244 being the sequence of denominators. With the exception of n < 2, the base 10 digital root of a(n) is always 4. In base 3 the digital root of a(n) is the same as the digital root of n. - Alonso del Arte, Jan 24 2006
The sequence 3*a(n), n >= 1, gives the number of edges of the Hanoi graph H_3^{n} with 3 pegs and n >= 1 discs. - Daniele Parisse, Jul 28 2006
Numbers n such that a(n) is prime are listed in A028491 = {3, 7, 13, 71, 103, 541, 1091, ...}. 2^(m+1) divides a(2^m*k) for m > 0. 5 divides a(4k). 5^2 divides a(20k). 7 divides a(6k). 7^2 divides a(42k). 11^2 divides a(5k). 13 divides a(3k). 17 divides a(16k). 19 divides a(18k). 1093 divides a(7k). 41 divides a(8k). p divides a((p-1)/5) for prime p = {41, 431, 491, 661, 761, 1021, 1051, 1091, 1171, ...}. p divides a((p-1)/4) for prime p = {13, 109, 181, 193, 229, 277, 313, 421, 433, 541, ...}. p divides a((p-1)/3) for prime p = {61, 67, 73, 103, 151, 193, 271, 307, 367, ...} = A014753, 3 and -3 are both cubes (one implies other) mod these primes p = 1 mod 6. p divides a((p-1)/2) for prime p = {11, 13, 23, 37, 47, 59, 61, 71, 73, 83, 97, ...} = A097933(n). p divides a(p-1) for prime p > 7. p^2 divides a(p*(p-1)k) for all prime p except p = 3. p^3 divides a(p*(p-1)*(p-2)k) for prime p = 11. - Alexander Adamchuk, Jan 22 2007
Let P(A) be the power set of an n-element set A. Then a(n) = the number of [unordered] pairs of elements {x,y} of P(A) for which x and y are disjoint [and both nonempty]. Wieder calls these "disjoint usual 2-combinations". - Ross La Haye, Jan 10 2008 [This is because each of the elements of {1, 2, ..., n} can be in the first subset, in the second or in neither. Because there are three options for each, the total number of options is 3^n. However, since the sets being empty is not an option we subtract 1 and since the subsets are unordered we then divide by 2! (The number of ways two objects can be arranged.) Thus we obtain (3^n-1)/2 = a(n). - Chayim Lowen, Mar 03 2015]
Also, still with P(A) being the power set of a n-element set A, a(n) is the number of 2-element subsets {x,y} of P(A) such that the union of x and y is equal to A. Cf. A341590. - Fabio Visonà, Feb 20 2021
Starting with offset 1 = binomial transform of A003945: (1, 3, 6, 12, 24, ...) and double bt of (1, 2, 1, 2, 1, 2, ...); equals polcoeff inverse of (1, -4, 3, 0, 0, 0, ...). - Gary W. Adamson, May 28 2009
Also the constant of the polynomials C(x) = 3x + 1 that form a sequence by performing this operation repeatedly and taking the result at each step as the input at the next. - Nishant Shukla (n.shukla722(AT)gmail.com), Jul 11 2009
It appears that this is A120444(3^n-1) = A004125(3^n) - A004125(3^n-1), where A004125 is the sum of remainders of n mod k for k = 1, 2, 3, ..., n. - John W. Layman, Jul 29 2009
Subsequence of A134025; A171960(a(n)) = a(n). - Reinhard Zumkeller, Jan 20 2010
Let A be the Hessenberg matrix of order n, defined by: A[1,j] = 1, A[i, i] := 3, (i > 1), A[i, i-1] = -1, and A[i, j] = 0 otherwise. Then, for n >= 1, a(n) = det(A). - Milan Janjic, Jan 27 2010
This is the sequence A(0, 1; 2, 3; 2) = A(0, 1; 4, -3; 0) of the family of sequences [a, b:c, d:k] considered by Gary Detlefs, and treated as A(a, b; c, d; k) in the Wolfdieter Lang link given below. - Wolfdieter Lang, Oct 18 2010
It appears that if s(n) is a first order rational sequence of the form s(0) = 0, s(n) = (2*s(n-1)+1)/(s(n-1)+2), n > 0, then s(n)= a(n)/(a(n)+1). - Gary Detlefs, Nov 16 2010
This sequence also describes the total number of moves to solve the [RED ; BLUE ; BLUE] or [RED ; RED ; BLUE] pre-colored Magnetic Towers of Hanoi puzzle (cf. A183111 - A183125).
From Adi Dani, Jun 08 2011: (Start)
a(n) is number of compositions of odd numbers into n parts less than 3. For example, a(3) = 13 and there are 13 compositions odd numbers into 3 parts < 3:
1: (0, 0, 1), (0, 1, 0), (1, 0, 0);
3: (0, 1, 2), (0, 2, 1), (1, 0, 2), (1, 2, 0), (2, 0, 1), (2, 1, 0), (1, 1, 1);
5: (1, 2, 2), (2, 1, 2), (2, 2, 1).
(End)
Pisano period lengths: 1, 2, 1, 2, 4, 2, 6, 4, 1, 4, 5, 2, 3, 6, 4, 8, 16, 2, 18, 4, ... . - R. J. Mathar, Aug 10 2012
a(n) is the total number of holes (triangles removed) after the n-th step of a Sierpiński triangle production. - Ivan N. Ianakiev, Oct 29 2013
a(n) solves Sum_{j = a(n) + 1 .. a(n+1)} j = k^2 for some integer k, given a(0) = 0 and requiring smallest a(n+1) > a(n). Corresponding k = 3^n. - Richard R. Forberg, Mar 11 2015
a(n+1) equals the number of words of length n over {0, 1, 2, 3} avoiding 01, 02 and 03. - Milan Janjic, Dec 17 2015
For n >= 1, a(n) is also the total number of words of length n, over an alphabet of three letters, such that one of the letters appears an odd number of times (See A006516 for 4 letter words, and the Balakrishnan reference there). - Wolfdieter Lang, Jul 16 2017
Also, the number of maximal cliques, maximum cliques, and cliques of size 4 in the n-Apollonian network. - Andrew Howroyd, Sep 02 2017
For n > 1, the number of triangles (cliques of size 3) in the (n-1)-Apollonian network. - Andrew Howroyd, Sep 02 2017
a(n) is the largest number that can be represented with n trits in balanced ternary. Correspondingly, -a(n) is the smallest number that can be represented with n trits in balanced ternary. - Thomas König, Apr 26 2020
These form Sierpinski nesting-stars, which alternate pattern on 3^n+1/2 star numbers A003154, based on the square configurations of 9^n. The partial sums of 3^n are delineated according to the geometry of a hexagram, see illustrations in links. (3*a(n-1) + 1) create Sierpinski-anti-triangles, representing the number of holes in a (n+1) Sierpinski triangle (see illustrations). - John Elias, Oct 18 2021
For n > 1, a(n) is the number of iterations necessary to calculate the hyperbolic functions with CORDIC. - Mathias Zechmeister, Jul 26 2022
a(n) is the least number k such that A065363(k) = n. - Amiram Eldar, Sep 03 2022
For all n >= 0, Sum_{k=a(n)+1..a(n+1)} 1/k < Sum_{j=a(n+1)+1..a(n+2)} 1/j. These are the minimal points which partition the infinite harmonic series into a monotonically increasing sequence. Each partition approximates log(3) from below as n tends to infinity. - Joseph Wheat, Apr 15 2023
a(n) is also the number of 3-cycles in the n-Dorogovtsev-Goltsev-Mendes graph (using the convention the 0-Dorogovtsev-Goltsev-Mendes graph is P_2). - Eric W. Weisstein, Dec 06 2023

Examples

			There are 4 3-block bicoverings of a 3-set: {{1, 2, 3}, {1, 2}, {3}}, {{1, 2, 3}, {1, 3}, {2}}, {{1, 2, 3}, {1}, {2, 3}} and {{1, 2}, {1, 3}, {2, 3}}.
Ternary........Decimal
0.................0
1.................1
11................4
111..............13
1111.............40 etc. - _Zerinvary Lajos_, Jan 14 2007
There are altogether a(3) = 13 three letter words over {A,B,C} with say, A, appearing an odd number of times: AAA; ABC, ACB, ABB, ACC; BAC, CAB, BAB, CAC; BCA, CBA, BBA, CCA. - _Wolfdieter Lang_, Jul 16 2017
		

References

  • J. G. Mauldon, Strong solutions for the counterfeit coin problem, IBM Research Report RC 7476 (#31437) 9/15/78, IBM Thomas J. Watson Research Center, P. O. Box 218, Yorktown Heights, N. Y. 10598.
  • Paulo Ribenboim, The Book of Prime Number Records, Springer-Verlag, NY, 2nd ed., 1989, p. 60.
  • Paulo Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, p. 53.
  • Amir Sapir, The Tower of Hanoi with Forbidden Moves, The Computer J. 47 (1) (2004) 20, case three-in-a row, sequence a(n).
  • Robert Sedgewick, Algorithms, 1992, pp. 109.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Sequences used for Shell sort: A033622, A003462, A036562, A036564, A036569, A055875.
Cf. A179526 (repeats), A113047 (characteristic function).
Cf. A000225, A000392, A004125, A014753, A028491 (indices of primes), A059443 (column k = 3), A065363, A097933, A120444, A321872 (sum reciprocals).
Cf. A064099 (minimal number of weightings to detect lighter or heavier coin among n coins).
Cf. A039755 (column k = 1).
Cf. A006516 (binomial transform, and special 4 letter words).
Cf. A341590.
Cf. A003462(n) (3-cycles), A367967(n) (5-cycles), A367968(n) (6-cycles).

Programs

Formula

G.f.: x/((1-x)*(1-3*x)).
a(n) = 4*a(n-1) - 3*a(n-2), n > 1. a(0) = 0, a(1) = 1.
a(n) = 3*a(n-1) + 1, a(0) = 0.
E.g.f.: (exp(3*x) - exp(x))/2. - Paul Barry, Apr 11 2003
a(n+1) = Sum_{k = 0..n} binomial(n+1, k+1)*2^k. - Paul Barry, Aug 20 2004
a(n) = Sum_{i = 0..n-1} 3^i, for n > 0; a(0) = 0.
a(n) = A125118(n, 2) for n > 1. - Reinhard Zumkeller, Nov 21 2006
a(n) = StirlingS2(n+1, 3) + StirlingS2(n+1, 2). - Ross La Haye, Jan 10 2008
a(n) = Sum_{k = 0..n} A106566(n, k)*A106233(k). - Philippe Deléham, Oct 30 2008
a(n) = 2*a(n-1) + 3*a(n-2) + 2, n > 1. - Gary Detlefs, Jun 21 2010
a(n) = 3*a(n-1) + a(n-2) - 3*a(n-3) = 5*a(n-1) - 7*a(n-2) + 3*a(n-3), a(0) = 0, a(1) = 1, a(2) = 4. Observation by G. Detlefs. See the W. Lang comment and link. - Wolfdieter Lang, Oct 18 2010
A008344(a(n)) = 0, for n > 1. - Reinhard Zumkeller, May 09 2012
A085059(a(n)) = 1 for n > 0. - Reinhard Zumkeller, Jan 31 2013
G.f.: Q(0)/2 where Q(k) = 1 - 1/(9^k - 3*x*81^k/(3*x*9^k - 1/(1 - 1/(3*9^k - 27*x*81^k/(9*x*9^k - 1/Q(k+1)))))); (continued fraction ). - Sergei N. Gladkovskii, Apr 12 2013
a(n) = A001065(3^n) where A001065(m) is the sum of the proper divisors of m for positive integer m. - Chayim Lowen, Mar 03 2015
a(n) = A000244(n) - A007051(n) = A007051(n)-1. - Yuchun Ji, Oct 23 2018
Sum_{n>=1} 1/a(n) = A321872. - Amiram Eldar, Nov 18 2020

Extensions

More terms from Michael Somos
Corrected my comment of Jan 10 2008. - Ross La Haye, Oct 29 2008
Removed comment that duplicated a formula. - Joerg Arndt, Mar 11 2010

A002718 Number of bicoverings of an n-set.

Original entry on oeis.org

1, 0, 1, 8, 80, 1088, 19232, 424400, 11361786, 361058000, 13386003873, 570886397340, 27681861184474, 1511143062540976, 92091641176725504, 6219762391554815200, 462595509951068027741, 37676170944802047077248, 3343539821715571537772071, 321874499078487207168905840
Offset: 0

Views

Author

Keywords

Comments

Another description: number of proper 2-covers of [1,...,n].

Examples

			For n=3, there are 8 collections of distinct subsets of {1,2,3} with the property that each of 1, 2, and 3 appears in exactly two subsets:
  {1,2,3},{1,2},{3}
  {1,2,3},{1,3},{2}
  {1,2,3},{2,3},{1}
  {1,2,3},{1},{2},{3}
  {1,2},{1,3},{2,3}
  {1,2},{1,3},{2},{3}
  {1,2},{2,3},{1},{3}
  {1,3},{2,3},{1},{2}
Therefore a(3) = 8. - _Michael B. Porter_, Jul 16 2016
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 303, #40.
  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    nmax = 16; imax = 2*(nmax - 2); egf := E^(-x - 1/2*x^2*(E^y - 1))*Sum[(x^i/i!)*E^(Binomial[i, 2]*y), {i, 0, imax}]; fx = CoefficientList[Series[egf, {y, 0, imax}], y]*Range[0, imax]!; a[n_] := Drop[ CoefficientList[ Series[fx[[n + 1]], {x, 0, imax}], x], 3] // Total; Table[ a[n] , {n, 2, nmax}] (* Jean-François Alcover, Apr 04 2013 *)

Formula

E.g.f. for k-block bicoverings of an n-set is exp(-x-1/2*x^2*(exp(y)-1))*Sum_{i=0..inf} x^i/i!*exp(binomial(i, 2)*y).
Stirling_2 transform of A060053.
The e.g.f.'s of A002718 (T(x)) and A060053 (V(x)) are related by T(x) = V(e^x-1).
a(n) = Sum_{m=0..n + floor(n/2); k=0..n; s=0..min(m/2,k); t=0..m-2s} Stirling2(n,k) * k!/m! * binomial(m,2s) * A001147(s) * (-1)^(m+s+t) * binomial(m-2s,t) * binomial(t*(t-1)/2,k-s). Interpret m as the number of blocks in a bicovering, k the number of clumps of points that are always all together in blocks. This formula counts bicoverings by quotienting them to the clumpless case (an operation which preserves degree) and counting incidence matrices of those, and counts those matrices as the transposes of incidence matrices of labeled graphs with no isolated points and no isolated edges. - David Pasino, Jul 09 2016

Extensions

More terms from Vladeta Jovovic, Feb 18 2001
a(0), a(1) prepended by Alois P. Heinz, Jul 29 2016

A060053 Number of r-bicoverings (or restricted proper 2-covers) of an n-set.

Original entry on oeis.org

1, 0, 1, 5, 43, 518, 8186, 163356, 3988342, 116396952, 3985947805, 157783127673, 7131072006829, 364166073164914, 20827961078794845, 1323968417981743817, 92917890994442697487, 7157607311779373890120, 602043767970637640566684
Offset: 0

Views

Author

Vladeta Jovovic, Feb 15 2001

Keywords

Comments

A bicovering is called an r-bicovering if the intersection of every two blocks contains at most one element.
Another name for this sequence is the number of restricted proper 2-covers of [1,...,n].
Number of T_0 2-regular set-systems on an n-set. - Andrew Howroyd, Jan 08 2020

Examples

			There are 5 r-bicoverings of a 3-set: 1 3-block bicovering {{1, 2}, {1, 3}, {2, 3}} and 4 4-block bicoverings {{1}, {2}, {3}, {1, 2, 3}}, {{2}, {3}, {1, 2}, {1, 3}}, {{1}, {3}, {1, 2}, {2, 3}}, {{1}, {2}, {1, 3}, {2, 3}}.
G.f. = 1 + x^2 + 5*x^3 + 43*x^4 + 518*x^5 + 8186*x^6 + 163356*x^7 + ...
		

References

  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983. (See p. 203.)

Crossrefs

Row 2 of A331039.
Row sums of A060052.

Programs

  • Maple
    A060053 := proc(n) local h, m; h := (m,n) -> add((-1/2)^k*binomial(m*(m-1)/2,n-k)/k!, k=0..n); n!*add(h(m,n)/m!, m=0..3*n); ceil(evalf(%/exp(1),99)) end: seq(A060053(i), i=0..18);
    # Caveat computator! Limited accuracy. Do not use it for n > 50. - Peter Luschny, Jul 06 2011
  • Mathematica
    f[n_] := FullSimplify[(n!/E)*Sum[(1/m!)*Sum[(-1/2)^k*Binomial[m*(m - 1)/2,
    n - k]/k!, {k, 0, n}], {m, 0, Infinity}]] (* Robert G. Wilson v, Jul 03 2011 *)
  • PARI
    a(n)=round(n!/exp(1)*sum(m=0,3*n+1,1/m!*sum(k=0,n,(-1/2)^k*binomial(m*(m-1)/2,n-k)/k!)))
    
  • PARI
    \\ here egf1 is A020556 as e.g.f.
    egf1(n)={my(bell=serlaplace(exp(exp(x + O(x^(2*n+1)))-1))); sum(i=0, n, sum(k=0, i, (-1)^k*binomial(i,k)*polcoef(bell, 2*i-k))*x^i/i!) + O(x*x^n)}
    seq(n)={my(A=egf1(n), B=log(1+x + O(x*x^n))/2); Vec(serlaplace(exp(-x/2 + O(x*x^n))*sum(k=0, n, polcoef(A,k)*B^k)))} \\ Andrew Howroyd, Jan 13 2020

Formula

E.g.f. for number of k-block r-bicoverings of an n-set is exp(-x-x^2*y/2)*Sum_{i=0..inf} (1+y)^binomial(i, 2)*x^i/i!.
a(n) = row sums of A060052.
Inverse binomial transform of A014500. - Vladeta Jovovic, Aug 22 2006
The e.g.f.'s of A002718 (T(x)) and A060053 (V(x)) are related by T(x) = V(e^x-1).
The e.g.f.'s of A014500 (U(x)) and A060053 (V(x)) are related by U(x) = e^x*V(x).
E.g.f.: exp(-x/2)*(Sum_{k>=0} A020556(k)*(log(1 + x)/2)^k/k!). - Andrew Howroyd, Jan 13 2020

A188445 T(n,k) is the number of (n*k) X k binary arrays with nonzero rows in decreasing order and n ones in every column.

Original entry on oeis.org

1, 2, 0, 5, 1, 0, 15, 8, 0, 0, 52, 80, 5, 0, 0, 203, 1088, 205, 1, 0, 0, 877, 19232, 11301, 278, 0, 0, 0, 4140, 424400, 904580, 67198, 205, 0, 0, 0, 21147, 11361786, 101173251, 24537905, 250735, 80, 0, 0, 0, 115975, 361058000, 15207243828, 13744869502
Offset: 1

Views

Author

R. H. Hardin, Mar 31 2011

Keywords

Examples

			Array begins:
============================================================================
n\k| 1 2 3   4       5          6             7              8             9
---+------------------------------------------------------------------------
1  | 1 2 5  15      52        203           877           4140         21147
2  | 0 1 8  80    1088      19232        424400       11361786     361058000
3  | 0 0 5 205   11301     904580     101173251    15207243828 2975725761202
4  | 0 0 1 278   67198   24537905   13744869502 11385203921707 ...
5  | 0 0 0 205  250735  425677958 1184910460297 ...
6  | 0 0 0  80  621348 5064948309 ...
7  | 0 0 0  15 1058139 ...
8  | 0 0 0   1 ...
...
Some solutions for 16 X 4:
  1 1 1 0    1 1 1 1    1 1 1 1    1 1 1 0    1 1 1 1
  1 0 1 1    1 1 0 1    1 1 0 0    1 0 1 1    1 1 0 0
  1 0 1 0    1 0 1 1    1 0 1 1    1 0 0 1    1 0 1 1
  1 0 0 1    1 0 0 0    1 0 0 0    1 0 0 0    1 0 0 0
  0 1 1 1    0 1 1 0    0 1 1 1    0 1 1 0    0 1 1 1
  0 1 0 1    0 1 0 0    0 1 0 0    0 1 0 1    0 1 0 0
  0 1 0 0    0 0 1 1    0 0 1 1    0 1 0 0    0 0 1 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 1 1    0 0 0 1
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
  0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0    0 0 0 0
		

Crossrefs

Columns 5..6 are A331127, A331129.
Column sums are A319190.

Programs

  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
    D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); WeighT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)}
    T(n, k)={my(m=n*k+1, q=Vec(exp(intformal(O(x^m) - x^n/(1-x)))/(1+x))); if(n==0, 1, (-1)^m*sum(j=0, m, my(s=0); forpart(p=j, s+=(-1)^#p*D(p, n, k), [1, n]); s*q[#q-j])/2)} \\ Andrew Howroyd, Dec 16 2018

Formula

A(n,k) = 0 for n > 2^(k-1). - Andrew Howroyd, Jan 24 2020

A060090 Number of ordered bicoverings of an unlabeled n-set.

Original entry on oeis.org

1, 0, 3, 23, 290, 4298, 79143, 1702923, 42299820, 1188147639, 37276597020, 1291633545897, 48995506718702, 2019395409175529, 89864601931874318, 4294295828157319651, 219321170795303112118, 11922219151375200468886
Offset: 0

Views

Author

Vladeta Jovovic, Feb 25 2001

Keywords

Examples

			There are 23 ordered bicoverings of an unlabeled 3-set, 7 3-block bicoverings:
1 ( { 3 }, { 1, 2 }, { 1, 2, 3 } )
2 ( { 3 }, { 1, 2, 3 }, { 1, 2 } )
3 ( { 2, 3 }, { 1 }, { 1, 2, 3 } )
4 ( { 2, 3 }, { 1, 3 }, { 1, 2 } )
5 ( { 2, 3 }, { 1, 2, 3 }, { 1 } )
6 ( { 1, 2, 3 }, { 3 }, { 1, 2 } )
7 ( { 1, 2, 3 }, { 2, 3 }, { 1 } )
and 16 4-block bicoverings:
1 ( { 3 }, { 2 }, { 1 }, { 1, 2, 3 } )
2 ( { 3 }, { 2 }, { 1, 3 }, { 1, 2 } )
3 ( { 3 }, { 2 }, { 1, 2 }, { 1, 3 } )
4 ( { 3 }, { 2 }, { 1, 2, 3 }, { 1 } )
5 ( { 3 }, { 2, 3 }, { 1 }, { 1, 2 } )
6 ( { 3 }, { 2, 3 }, { 1, 2 }, { 1 } )
7 ( { 3 }, { 1, 2 }, { 2 }, { 1, 3 } )
8 ( { 3 }, { 1, 2 }, { 2, 3 }, { 1 } )
9 ( { 3 }, { 1, 2, 3 }, { 2 }, { 1 } )
10 ( { 2, 3 }, { 3 }, { 1 }, { 1, 2 } )
11 ( { 2, 3 }, { 3 }, { 1, 2 }, { 1 } )
12 ( { 2, 3 }, { 1 }, { 3 }, { 1, 2 } )
13 ( { 2, 3 }, { 1 }, { 1, 3 }, { 2 } )
14 ( { 2, 3 }, { 1, 3 }, { 2 }, { 1 } )
15 ( { 2, 3 }, { 1, 3 }, { 1 }, { 2 } )
16 ( { 1, 2, 3 }, { 3 }, { 2 }, { 1 } )
		

Crossrefs

Row n=2 of A331571.
Row sums of A060092.

Programs

  • PARI
    seq(n)={my(m=3*n\2, y='y + O('y^(n+1))); Vec(subst(Pol(serlaplace(exp(-x - x^2*y/(2*(1-y)) + O(x*x^m))*sum(k=0, m, 1/(1-y)^binomial(k, 2)*x^k/k!))), x, 1))} \\ Andrew Howroyd, Jan 30 2020

Formula

E.g.f. for ordered k-block bicoverings of an unlabeled n-set is exp(-x-x^2/2*y/(1-y)) * Sum_{k>=0} 1/(1-y)^binomial(k,2)*x^k/k!.

A060052 Triangle read by rows: T(n,k) gives number of r-bicoverings of an n-set with k blocks, n >= 2, k = 3..n+floor(n/2).

Original entry on oeis.org

1, 1, 4, 0, 15, 25, 3, 0, 30, 222, 226, 40, 0, 30, 1230, 3670, 2706, 535, 15, 0, 0, 5040, 39900, 69450, 40405, 8141, 420, 0, 0, 15120, 345240, 1254960, 1498035, 722275, 142877, 9730, 105, 0, 0, 30240, 2492280, 18587520, 40701780, 36450820, 15031204, 2871240, 226828, 5040
Offset: 2

Views

Author

Vladeta Jovovic, Feb 15 2001

Keywords

Comments

A bicovering is r-bicovering if intersection of every two blocks contains at most one element.

Examples

			Triangle starts:
[1],
[1, 4],
[0, 15, 25, 3],
[0, 30, 222, 226, 40],
[0, 30, 1230, 3670, 2706, 535, 15],
[0, 0, 5040, 39900, 69450, 40405, 8141, 420],
[0, 0, 15120, 345240, 1254960, 1498035, 722275, 142877, 9730, 105],
[0, 0, 30240, 2492280, 18587520, 40701780, 36450820, 15031204, 2871240, 226828, 5040],
...
		

References

  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.

Crossrefs

Row sums are A060053.
Column sums are A060051.

Programs

  • PARI
    \\ returns k-th column as vector.
    C(k)=if(k<3, [], Vecrev(serlaplace(polcoef(exp(-x-1/2*x^2*y + O(x*x^k))*sum(i=0, 3*k\2, (1+y)^binomial(i, 2)*x^i/i!), k))/y)) \\ Andrew Howroyd, Jan 30 2020
    
  • PARI
    T(n)={my(m=(3*n\2), y='y + O('y^(n+1))); my(g=exp(-x-1/2*x^2*y + O(x*x^m))*sum(k=0, m, (1+y)^binomial(k, 2)*x^k/k!)); Mat([Col(serlaplace(p), -n) | p<-Vec(g)[2..m+1]])}
    { my(A=T(8)); for(n=2, matsize(A)[1], print(A[n, 3..3*n\2])) } \\ Andrew Howroyd, Jan 30 2020

Formula

E.g.f.: A(x, y) = exp(-x-1/2*x^2*y)*Sum_{i>=0} (1+y)^binomial(i, 2)*x^i/i!.
T(n, k) = (n!/k!) * A276640(k, n). - David Pasino, Sep 22 2016
T(n,k) = 0 for n > binomial(k,2). - Andrew Howroyd, Jan 30 2020

Extensions

Zeros inserted into data by Andrew Howroyd, Jan 30 2020

A060487 Triangle T(n,k) of k-block tricoverings of an n-set (n >= 3, k >= 4).

Original entry on oeis.org

1, 3, 1, 7, 57, 95, 43, 3, 35, 717, 3107, 4520, 2465, 445, 12, 155, 7845, 75835, 244035, 325890, 195215, 50825, 4710, 70, 651, 81333, 1653771, 10418070, 27074575, 33453959, 20891962, 6580070, 965965, 52430, 465
Offset: 3

Views

Author

Vladeta Jovovic, Mar 20 2001

Keywords

Comments

A covering of a set is a tricovering if every element of the set is covered by exactly three blocks of the covering.

Examples

			Triangle begins:
  [1, 3, 1];
  [7, 57, 95, 43, 3];
  [35, 717, 3107, 4520, 2465, 445, 12];
  [155, 7845, 75835, 244035, 325890, 195215, 50825, 4710, 70];
  [651, 81333, 1653771, 10418070, 27074575, 33453959, 20891962, 6580070, 965965, 52430, 465];
   ...
There are 205 tricoverings of a 4-set(cf. A060486): 7 4-block, 57 5-block, 95 6-block, 43 7-block and 3 8-block tricoverings.
		

Crossrefs

Programs

  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
    D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); WeighT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)}
    row(n, k)={my(m=n*k+1, q=Vec(exp(intformal(O(x^m) - x^n/(1-x)))/(y+x))); if(n==0, 1, (-1)^m*sum(j=0, m, my(s=0); forpart(p=j, s+=(-1)^#p*D(p, n, k), [1, n]); s*q[#q-j])*y^(m-n)/(1+y))}
    for(n=3, 8, print(Vecrev(row(3,n)))); \\ Andrew Howroyd, Dec 23 2018

Formula

E.g.f. for k-block tricoverings of an n-set is exp(-x+x^2/2+(exp(y)-1)*x^3/3)*Sum_{k=0..inf}x^k/k!*exp(-1/2*x^2*exp(k*y))*exp(binomial(k, 3)*y).

A059945 Number of 4-block bicoverings of an n-set.

Original entry on oeis.org

0, 0, 4, 39, 280, 1815, 11284, 68859, 416560, 2509455, 15086764, 90610179, 543928840, 3264374295, 19588645444, 117539063499, 705255937120, 4231600258335, 25389795391324, 152339353740819, 914037866361400, 5484232429393575, 32905410268988404, 197432508689714139
Offset: 1

Views

Author

Vladeta Jovovic, Feb 14 2001

Keywords

Examples

			There are 4 4-block bicoverings of a 3-set: {{1},{2},{3},{1,2,3}}, {{2},{3},{1,2},{1,3}}, {{1},{3},{1,2},{2,3}} and {{1},{2},{1,3},{2,3}}.
		

References

  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.

Crossrefs

Column k=4 of A059443.
Cf. A002718.

Programs

  • Mathematica
    With[{c=1/4!},Table[c(6^n-4 3^n-3 2^n+12),{n,20}]] (* or *) LinearRecurrence[ {12,-47,72,-36},{0,0,4,39},20] (* Harvey P. Dale, Aug 10 2011 *)
  • PARI
    a(n) = {(1/4!)*(6^n - 4*3^n - 3*2^n + 12)} \\ Andrew Howroyd, Jan 29 2020

Formula

a(n) = (1/4!)*(6^n - 4*3^n - 3*2^n + 12).
E.g.f. for m-block bicoverings of an n-set is exp(-x-1/2*x^2*(exp(y)-1))*Sum_{i=0..inf} x^i/i!*exp(binomial(i, 2)*y).
a(n) = 12*a(n-1) - 47*a(n-2) + 72*a(n-3) - 36*a(n-4) for n > 4. - Harvey P. Dale, Aug 10 2011
G.f.: -x^3*(9*x-4) / ((x-1)*(2*x-1)*(3*x-1)*(6*x-1)). - Colin Barker, Jan 11 2013

Extensions

More terms from Colin Barker, Jan 11 2013

A060092 Triangle T(n,k) of k-block ordered bicoverings of an unlabeled n-set, n >= 2, k = 3..n+floor(n/2).

Original entry on oeis.org

3, 7, 16, 12, 63, 125, 90, 18, 162, 722, 1716, 1680, 25, 341, 2565, 11350, 27342, 29960, 7560, 33, 636, 7180, 49860, 208302, 503000, 631512, 302400, 42, 1092, 17335, 173745, 1099602, 4389875, 10762299, 14975730, 9632700, 1247400
Offset: 2

Views

Author

Vladeta Jovovic, Feb 26 2001

Keywords

Comments

All columns are polynomials of order binomial(k, 2). - Andrew Howroyd, Jan 30 2020

Examples

			[3],
[7, 16],
[12, 63, 125, 90],
[18, 162, 722, 1716, 1680],
[25, 341, 2565, 11350, 27342, 29960, 7560],
[33, 636, 7180, 49860, 208302, 503000, 631512, 302400],
[42, 1092, 17335, 173745, 1099602, 4389875, 10762299, 14975730, 9632700, 1247400], ...
There are 23=7+16 ordered bicoverings of an unlabeled 3-set: 7 3-block bicoverings and 16 4-block bicoverings, cf. A060090.
		

Crossrefs

Row sums are A060090.
Columns k=3..7 are A055998(n-1), A060091, A060093, A060094, A060095.

Programs

  • PARI
    \\ gives g.f. of k-th column.
    ColGf(k) = k!*polcoef(exp(-x - x^2*y/(2*(1-y)) + O(x*x^k))*sum(j=0, k, 1/(1-y)^binomial(j, 2)*x^j/j!), k) \\ Andrew Howroyd, Jan 30 2020
    
  • PARI
    T(n)={my(m=(3*n\2), y='y + O('y^(n+1))); my(g=serlaplace(exp(-x - x^2*y/(2*(1-y)) + O(x*x^m))*sum(k=0, m, 1/(1-y)^binomial(k, 2)*x^k/k!))); Mat([Col(p/y^2, -n) | p<-Vec(g)[2..m+1]])}
    { my(A=T(8)); for(n=2, matsize(A)[1], print(A[n, 3..3*n\2])) } \\ Andrew Howroyd, Jan 30 2020

Formula

E.g.f. for k-block ordered bicoverings of an unlabeled n-set is exp(-x-x^2/2*y/(1-y))*Sum_{k=0..inf} 1/(1-y)^binomial(k, 2)*x^k/k!.

A059951 Number of 10-block bicoverings of an n-set.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 420, 154637, 20368816, 1775801814, 124151410020, 7596257673279, 426319554841752, 22564352299016528, 1146221298547133380, 56531610963314602401, 2728475248127447671008, 129586638359127411410442, 6080467290450346517206500, 282689089820505452872162403
Offset: 1

Views

Author

Vladeta Jovovic, Feb 14 2001

Keywords

References

  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.

Crossrefs

Column k=10 of A059443.
Cf. A002718.

Formula

a(n) = (1/10!)*(45^n - 10*36^n - 45*29^n + 90*28^n + 360*22^n - 480*21^n + 630*17^n - 2520*16^n + 2100*15^n - 3780*12^n + 10080*11^n - 6552*10^n - 3150*9^n + 18900*8^n - 31500*7^n + 28560*6^n - 46620*5^n + 27720*4^n + 85560*3^n - 146160*2^n + 83520).
E.g.f. for m-block bicoverings of an n-set is exp(-x-1/2*x^2*(exp(y)-1))*Sum_{i=0..inf} x^i/i!*exp(binomial(i, 2)*y).
G.f.: -x^7*(5467233152463667200*x^14 -6460773223081605120*x^13 +3312489509664336576*x^12 -965946275708647680*x^11 +175045400422088532*x^10 -19853467917718628*x^9 +1255863452001343*x^8 -11591551437545*x^7 -5424120630669*x^6 +520759916751*x^5 -24697320639*x^4 +659527325*x^3 -8843563*x^2 +25697*x +420) / ((x -1)*(2*x -1)*(3*x -1)*(4*x -1)*(5*x -1)*(6*x -1)*(7*x -1)*(8*x -1)*(9*x -1)*(10*x -1)*(11*x -1)*(12*x -1)*(15*x -1)*(16*x -1)*(17*x -1)*(21*x -1)*(22*x -1)*(28*x -1)*(29*x -1)*(36*x -1)*(45*x -1)). - Colin Barker, Jul 09 2013

Extensions

More terms from Colin Barker, Jul 09 2013
Showing 1-10 of 26 results. Next