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-8 of 8 results.

A026898 a(n) = Sum_{k=0..n} (n-k+1)^k.

Original entry on oeis.org

1, 2, 4, 9, 23, 66, 210, 733, 2781, 11378, 49864, 232769, 1151915, 6018786, 33087206, 190780213, 1150653921, 7241710930, 47454745804, 323154696185, 2282779990495, 16700904488706, 126356632390298, 987303454928973, 7957133905608837, 66071772829247410
Offset: 0

Views

Author

Keywords

Comments

Row sums of A004248, A009998, A009999.
First differences are in A047970.
First differences of A103439.
Antidiagonal sums of array A003992.
a(n-1), for n>=1, is the number of length-n restricted growth strings (RGS) [s(0),s(1),...,s(n-1)] where s(0)=0 and s(k)<=1+max(prefix) for k>=1, that are simultaneously projections as maps f: [n] -> [n] where f(x)<=x and f(f(x))=f(x); see example and the two comments (Arndt, Apr 30 2011 Jan 04 2013) in A000110. - Joerg Arndt, Mar 07 2015
Number of finite sequences s of length n+1 whose discriminator sequence is s itself. Here the discriminator sequence of s is the one where the n-th term (n>=1) is the least positive integer k such that the first n terms are pairwise incongruent, modulo k. - Jeffrey Shallit, May 17 2016
From Gus Wiseman, Jan 08 2019: (Start)
Also the number of set partitions of {1,...,n+1} whose minima form an initial interval of positive integers. For example, the a(3) = 9 set partitions are:
{{1},{2},{3},{4}}
{{1},{2},{3,4}}
{{1},{2,4},{3}}
{{1,4},{2},{3}}
{{1},{2,3,4}}
{{1,3},{2,4}}
{{1,4},{2,3}}
{{1,3,4},{2}}
{{1,2,3,4}}
Missing from this list are:
{{1},{2,3},{4}}
{{1,2},{3},{4}}
{{1,3},{2},{4}}
{{1,2},{3,4}}
{{1,2,3},{4}}
{{1,2,4},{3}}
(End)
a(n) is the number of m-tuples of nonnegative integers less than or equal to n-m (including the "0-tuple"). - Mathew Englander, Apr 11 2021

Examples

			G.f.: A(x) = 1 + 2*x + 4*x^2 + 9*x^3 + 23*x^4 + 66*x^5 + 210*x^6 + ...
where we have the identity:
A(x) = 1/(1-x) + x/(1-2*x) + x^2/(1-3*x) + x^3/(1-4*x) + x^4/(1-5*x) + ...
is equal to
A(x) = 1/(1-x) + x/((1-x)^2*(1+x)) + 2!*x^2/((1-x)^3*(1+x)*(1+2*x)) + 3!*x^3/((1-x)^4*(1+x)*(1+2*x)*(1+3*x)) + 4!*x^4/((1-x)^5*(1+x)*(1+2*x)*(1+3*x)*(1+4*x)) + ...
From _Joerg Arndt_, Mar 07 2015: (Start)
The a(5-1) = 23 RGS described in the comment are (dots denote zeros):
01:  [ . . . . . ]
02:  [ . 1 . . . ]
03:  [ . 1 . . 1 ]
04:  [ . 1 . 1 . ]
05:  [ . 1 . 1 1 ]
06:  [ . 1 1 . . ]
07:  [ . 1 1 . 1 ]
08:  [ . 1 1 1 . ]
09:  [ . 1 1 1 1 ]
10:  [ . 1 2 . . ]
11:  [ . 1 2 . 1 ]
12:  [ . 1 2 . 2 ]
13:  [ . 1 2 1 . ]
14:  [ . 1 2 1 1 ]
15:  [ . 1 2 1 2 ]
16:  [ . 1 2 2 . ]
17:  [ . 1 2 2 1 ]
18:  [ . 1 2 2 2 ]
19:  [ . 1 2 3 . ]
20:  [ . 1 2 3 1 ]
21:  [ . 1 2 3 2 ]
22:  [ . 1 2 3 3 ]
23:  [ . 1 2 3 4 ]
(End)
		

Crossrefs

Programs

  • Haskell
    a026898 n = sum $ zipWith (^) [n + 1, n .. 1] [0 ..]
    -- Reinhard Zumkeller, Sep 14 2014
    
  • Magma
    [(&+[(n-k+1)^k: k in [0..n]]): n in [0..50]]; // Stefano Spezia, Jan 09 2019
    
  • Maple
    a:= n-> add((n+1-j)^j, j=0..n): seq(a(n), n=0..23); # Zerinvary Lajos, Apr 18 2009
  • Mathematica
    Table[Sum[(n-k+1)^k, {k,0,n}], {n, 0, 25}] (* Michael De Vlieger, Apr 01 2015 *)
  • PARI
    {a(n)=polcoeff(sum(m=0,n,x^m/(1-(m+1)*x+x*O(x^n))),n)} /* Paul D. Hanna, Sep 13 2011 */
    
  • PARI
    {INTEGRATE(n,F)=local(G=F);for(i=1,n,G=intformal(G));G}
    {a(n)=local(A=1+x);A=sum(k=0,n,INTEGRATE(k,exp((k+1)*x+x*O(x^n))));n!*polcoeff(A,n)} \\ Paul D. Hanna, Dec 28 2013
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n)=polcoeff( sum(m=0, n, m!*x^m/(1-x +x*O(x^n))^(m+1)/prod(k=1, m, 1+k*x +x*O(x^n))), n)}  /* From o.g.f. (Paul D. Hanna, Jul 20 2014) */
    for(n=0, 25, print1(a(n), ", "))
    
  • Sage
    [sum((n-j+1)^j for j in (0..n)) for n in (0..30)] # G. C. Greubel, Jun 15 2021

Formula

a(n) = A003101(n) + 1.
G.f.: Sum_{n>=0} x^n/(1 - (n+1)*x). - Paul D. Hanna, Sep 13 2011
G.f.: G(0) where G(k) = 1 + x*(2*k*x-1)/((2*k*x+x-1) - x*(2*k*x+x-1)^2/(x*(2*k*x+x-1) + (2*k*x+2*x-1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 26 2013
E.g.f.: Sum_{n>=0} Integral^n exp((n+1)*x) dx^n, where Integral^n F(x) dx^n is the n-th integration of F(x) with no constant of integration. - Paul D. Hanna, Dec 28 2013
O.g.f.: Sum_{n>=0} n! * x^n/(1-x)^(n+1) / Product_{k=1..n} (1 + k*x). - Paul D. Hanna, Jul 20 2014
a(n) = A101494(n+1,0). - Vladimir Kruchinin, Apr 01 2015
a(n-1) = Sum_{k = 1..n} k^(n-k). - Gus Wiseman, Jan 08 2019
log(a(n)) ~ (1 - 1/LambertW(exp(1)*n)) * n * log(1 + n/LambertW(exp(1)*n)). - Vaclav Kotesovec, Jun 15 2021
a(n) ~ sqrt(2*Pi/(n+1 + (n+1)/w(n))) * ((n+1)/w(n))^(n+2 - (n+1)/w(n)), where w(n) = LambertW(exp(1)*(n+1)). - Vaclav Kotesovec, Jun 25 2021, after user "leonbloy", see Mathematics Stack Exchange link.

Extensions

a(23)-a(25) from Paul D. Hanna, Dec 28 2013

A009998 Triangle in which j-th entry in i-th row is (j+1)^(i-j).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 4, 1, 1, 16, 27, 16, 5, 1, 1, 32, 81, 64, 25, 6, 1, 1, 64, 243, 256, 125, 36, 7, 1, 1, 128, 729, 1024, 625, 216, 49, 8, 1, 1, 256, 2187, 4096, 3125, 1296, 343, 64, 9, 1, 1, 512, 6561, 16384, 15625, 7776, 2401, 512, 81, 10, 1
Offset: 0

Views

Author

Keywords

Comments

Read as a square array this is the Hilbert transform of triangle A123125 (see A145905 for the definition of this term). For example, the fourth row of A123125 is (0,1,4,1) and the expansion (x + 4*x^2 + x^3)/(1-x)^4 = x + 8*x^2 + 27*x^3 + 64*x^4 + ... generates the entries in the fourth row of this array read as a square. - Peter Bala, Oct 28 2008

Examples

			Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  3,  1;
  1,  8,  9,  4,  1;
  1, 16, 27, 16,  5,  1;
  1, 32, 81, 64, 25,  6,  1;
  ...
From _Gus Wiseman_, May 01 2021: (Start)
The rows of the triangle are obtained by reading antidiagonals upward in the following table of A(k,n) = n^k, with offset k = 0, n = 1:
         n=1:     n=2:     n=3:     n=4:     n=5:     n=6:
   k=0:   1        1        1        1        1        1
   k=1:   1        2        3        4        5        6
   k=2:   1        4        9       16       25       36
   k=3:   1        8       27       64      125      216
   k=4:   1       16       81      256      625     1296
   k=5:   1       32      243     1024     3125     7776
   k=6:   1       64      729     4096    15625    46656
   k=7:   1      128     2187    16384    78125   279936
   k=8:   1      256     6561    65536   390625  1679616
   k=9:   1      512    19683   262144  1953125 10077696
  k=10:   1     1024    59049  1048576  9765625 60466176
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 24.

Crossrefs

Row sums give A026898.
Column n = 2 of the array is A000079.
Column n = 3 of the array is A000244.
Row k = 2 of the array is A000290.
Row k = 3 of the array is A000578.
Diagonal n = k of the array is A000312.
Diagonal n = k + 1 of the array is A000169.
Diagonal n = k + 2 of the array is A000272.
The transpose of the array is A009999.
The numbers of divisors of the entries are A343656 (row sums: A343657).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Haskell
    a009998 n k = (k + 1) ^ (n - k)
    a009998_row n = a009998_tabl !! n
    a009998_tabl = map reverse a009999_tabl
    -- Reinhard Zumkeller, Feb 02 2014
    
  • Maple
    E := (n,x) -> `if`(n=0,1,x*(1-x)*diff(E(n-1,x),x)+E(n-1,x)*(1+(n-1)*x));
    G := (n,x) -> E(n,x)/(1-x)^(n+1);
    A009998 := (n,k) -> coeff(series(G(n-k,x),x,18),x,k);
    seq(print(seq(A009998(n,k),k=0..n)),n=0..6);
    # Peter Luschny, Aug 02 2010
  • Mathematica
    Flatten[Table[(j+1)^(i-j),{i,0,20},{j,0,i}]] (* Harvey P. Dale, Dec 25 2012 *)
  • PARI
    T(i,j)=(j+1)^(i-j) \\ Charles R Greathouse IV, Feb 06 2017

Formula

T(n,n) = 1; T(n,k) = (k+1)*T(n-1,k) for k=0..n-1. - Reinhard Zumkeller, Feb 02 2014
T(n,m) = (m+1)*Sum_{k=0..n-m}((n+1)^(k-1)*(n-m)^(n-m-k)*(-1)^(n-m-k)*binomial(n-m-1,k-1)). - Vladimir Kruchinin, Sep 12 2015

Extensions

a(62) corrected to 512 by T. D. Noe, Dec 20 2007

A003992 Square array read by upwards antidiagonals: T(n,k) = n^k for n >= 0, k >= 0.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 4, 1, 0, 1, 4, 9, 8, 1, 0, 1, 5, 16, 27, 16, 1, 0, 1, 6, 25, 64, 81, 32, 1, 0, 1, 7, 36, 125, 256, 243, 64, 1, 0, 1, 8, 49, 216, 625, 1024, 729, 128, 1, 0, 1, 9, 64, 343, 1296, 3125, 4096, 2187, 256, 1, 0, 1, 10, 81, 512, 2401, 7776, 15625, 16384, 6561, 512, 1, 0
Offset: 0

Views

Author

Keywords

Comments

If the array is transposed, T(n,k) is the number of oriented rows of n colors using up to k different colors. The formula would be T(n,k) = [n==0] + [n>0]*k^n. The generating function for column k would be 1/(1-k*x). For T(3,2)=8, the rows are AAA, AAB, ABA, ABB, BAA, BAB, BBA, and BBB. - Robert A. Russell, Nov 08 2018
T(n,k) is the number of multichains of length n from {} to [k] in the Boolean lattice B_k. - Geoffrey Critzer, Apr 03 2020

Examples

			Rows begin:
[1, 0,  0,   0,    0,     0,      0,      0, ...],
[1, 1,  1,   1,    1,     1,      1,      1, ...],
[1, 2,  4,   8,   16,    32,     64,    128, ...],
[1, 3,  9,  27,   81,   243,    729,   2187, ...],
[1, 4, 16,  64,  256,  1024,   4096,  16384, ...],
[1, 5, 25, 125,  625,  3125,  15625,  78125, ...],
[1, 6, 36, 216, 1296,  7776,  46656, 279936, ...],
[1, 7, 49, 343, 2401, 16807, 117649, 823543, ...], ...
		

Crossrefs

Main diagonal is A000312. Other diagonals include A000169, A007778, A000272, A008788. Antidiagonal sums are in A026898.
Cf. A099555.
Transpose is A004248. See A051128, A095884, A009999 for other versions.
Cf. A277504 (unoriented), A293500 (chiral).

Programs

  • Magma
    [[(n-k)^k: k in [0..n]]: n in [0..10]]; // G. C. Greubel, Nov 08 2018
  • Mathematica
    Table[If[k == 0, 1, (n - k)^k], {n, 0, 11}, {k, 0, n}]//Flatten
  • PARI
    T(n,k) = (n-k)^k \\ Charles R Greathouse IV, Feb 07 2017
    

Formula

E.g.f.: Sum T(n,k)*x^n*y^k/k! = 1/(1-x*exp(y)). - Paul D. Hanna, Oct 22 2004
E.g.f.: Sum T(n,k)*x^n/n!*y^k/k! = e^(x*e^y). - Franklin T. Adams-Watters, Jun 23 2006

Extensions

More terms from David W. Wilson
Edited by Paul D. Hanna, Oct 22 2004

A333988 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where column k is the expansion of (1-(k+1)*x) / (1-2*(k+1)*x+((k-1)*x)^2).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 8, 1, 1, 4, 17, 32, 1, 1, 5, 28, 99, 128, 1, 1, 6, 41, 208, 577, 512, 1, 1, 7, 56, 365, 1552, 3363, 2048, 1, 1, 8, 73, 576, 3281, 11584, 19601, 8192, 1, 1, 9, 92, 847, 6016, 29525, 86464, 114243, 32768, 1, 1, 10, 113, 1184, 10033, 62976, 265721, 645376, 665857, 131072, 1
Offset: 0

Views

Author

Seiichi Manyama, Sep 04 2020

Keywords

Examples

			Square array begins:
  1,   1,    1,     1,     1,     1, ...
  1,   2,    3,     4,     5,     6, ...
  1,   8,   17,    28,    41,    56, ...
  1,  32,   99,   208,   365,   576, ...
  1, 128,  577,  1552,  3281,  6016, ...
  1, 512, 3363, 11584, 29525, 62976, ...
		

Crossrefs

Main diagonal gives A333990.

Programs

  • Mathematica
    T[n_, 0] := 1; T[n_, k_] := Sum[k^j * Binomial[2*n, 2*j], {j, 0, n}]; Table[T[k, n - k], {n, 0, 10}, {k, 0, n}] // Flatten (* Amiram Eldar, Sep 04 2020 *)
  • PARI
    {T(n, k) = sum(j=0, n, k^j*binomial(2*n, 2*j))}

Formula

T(n,k) = Sum_{j=0..n} k^j * binomial(2*n,2*j).
T(0,k) = 1, T(1,k) = k+1 and T(n,k) = 2 * (k+1) * T(n-1,k) - (k-1)^2 * T(n-2,k) for n>1.

A095884 Triangle read by rows: T(n,k) = (n-k)^k, n>=1, 1<=k<=n.

Original entry on oeis.org

0, 1, 0, 2, 1, 0, 3, 4, 1, 0, 4, 9, 8, 1, 0, 5, 16, 27, 16, 1, 0, 6, 25, 64, 81, 32, 1, 0, 7, 36, 125, 256, 243, 64, 1, 0, 8, 49, 216, 625, 1024, 729, 128, 1, 0, 9, 64, 343, 1296, 3125, 4096, 2187, 256, 1, 0, 10, 81, 512, 2401, 7776, 15625, 16384, 6561, 512, 1, 0, 11, 100
Offset: 1

Views

Author

Herman Jamke (hermanjamke(AT)fastmail.fm), Jul 10 2004

Keywords

Examples

			0
1 0
2 1 0
3 4 1 0
4 9 8 1 0
5 16 27 16 1 0
6 25 64 81 32 1 0
7 36 125 256 243 64 1 0
8 49 216 625 1024 729 128 1 0
9 64 343 1296 3125 4096 2187 256 1 0
		

Crossrefs

A046688 Antidiagonals of square array in which k-th row (k>0) is an arithmetic progression of difference 2^(k-1).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 5, 1, 1, 5, 7, 9, 9, 1, 1, 6, 9, 13, 17, 17, 1, 1, 7, 11, 17, 25, 33, 33, 1, 1, 8, 13, 21, 33, 49, 65, 65, 1, 1, 9, 15, 25, 41, 65, 97, 129, 129, 1, 1, 10, 17, 29, 49, 81, 129, 193, 257, 257, 1, 1, 11, 19, 33, 57, 97, 161, 257, 385, 513, 513, 1
Offset: 0

Views

Author

Keywords

Examples

			From _Gus Wiseman_, May 08 2021: (Start):
Array A(m,n) = 1 + n*2^(m-1) begins:
       n=0: n=1: n=2: n=3: n=4: n=5: n=6: n=7: n=8: n=9:
  m=0:   1    1    1    1    1    1    1    1    1    1
  m=1:   1    2    3    5    9   17   33   65  129  257
  m=2:   1    3    5    9   17   33   65  129  257  513
  m=3:   1    4    7   13   25   49   97  193  385  769
  m=4:   1    5    9   17   33   65  129  257  513 1025
  m=5:   1    6   11   21   41   81  161  321  641 1281
  m=6:   1    7   13   25   49   97  193  385  769 1537
  m=7:   1    8   15   29   57  113  225  449  897 1793
  m=8:   1    9   17   33   65  129  257  513 1025 2049
  m=9:   1   10   19   37   73  145  289  577 1153 2305
Triangle T(n,k) = 1 + (n-k)*2^(k-1) begins:
   1
   1   1
   1   2   1
   1   3   3   1
   1   4   5   5   1
   1   5   7   9   9   1
   1   6   9  13  17  17   1
   1   7  11  17  25  33  33   1
   1   8  13  21  33  49  65  65   1
   1   9  15  25  41  65  97 129 129   1
   1  10  17  29  49  81 129 193 257 257   1
   1  11  19  33  57  97 161 257 385 513 513   1
(End)
		

References

  • G. H. Hardy, A Theorem Concerning the Infinite Cardinal Numbers, Quart. J. Math., 35 (1904), p. 90 = Collected Papers, Vol. VII, p. 430.

Crossrefs

Row sums are A000079.
Diagonal n = m + 1 of the array is A002064.
Diagonal n = m of the array is A005183.
Column m = 1 of the array is A094373.
Diagonal n = m - 1 of the array is A131056.
A002109 gives hyperfactorials (sigma: A260146, omega: A303281).
A009998(k,n) = n^k.
A009999(n,k) = n^k.
A057156 = (2^n)^(2^n).
A062319 counts divisors of n^n.

Programs

  • Mathematica
    Table[If[k==0,1,n*2^(k-1)+1],{n,0,9},{k,0,9}] (* ARRAY, Gus Wiseman, May 08 2021 *)
    Table[If[k==0,1,1+(n-k)*2^(k-1)],{n,0,10},{k,0,n}] (* TRIANGLE, Gus Wiseman, May 08 2021 *)
  • PARI
    A(m,n)={if(m>0, 1+n*2^(m-1), 1)}
    { for(m=0, 10, for(n=0, 10, print1(A(m,n), ", ")); print) } \\ Andrew Howroyd, Mar 07 2020

Formula

A(m,n) = 1 + n*2^(m-1) for m > 1. - Andrew Howroyd, Mar 07 2020
As a triangle, T(n,k) = A(k,n-k) = 1 + (n-k)*2^(k-1). - Gus Wiseman, May 08 2021

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Apr 06 2000

A329940 Square array read by antidiagonals upwards: T(n,k) is the number of right unique relations between set A with n elements and set B with k elements.

Original entry on oeis.org

1, 3, 2, 7, 8, 3, 15, 26, 15, 4, 31, 80, 63, 24, 5, 63, 242, 255, 124, 35, 6, 127, 728, 1023, 624, 215, 48, 7, 255, 2186, 4095, 3124, 1295, 342, 63, 8, 511, 6560, 16383, 15624, 7775, 2400, 511, 80, 9, 1023, 19682, 65535, 78124, 46655, 16806, 4095, 728, 99, 10
Offset: 1

Views

Author

Roy S. Freedman, Nov 24 2019

Keywords

Comments

A relation R between set A with n elements and set B with k elements is a subset of the Cartesian product A x B. A relation R is right unique if (a, b1) in R and (a,b2) in R implies b1=b2. T(n,k) is the number of right unique relations and T(k,n) is the number of left unique relations: relation R is left unique if (a1,b) in R and (a2,b) in R implies a1=a2.

Examples

			T(n,k) begins:
    1,    2,     3,      4,       5,       6,        7,        8, ...
    3,    8,    15,     24,      35,      48,       63,       80, ...
    7,   26,    63,    124,     215,     342,      511,      728, ...
   15,   80,   255,    624,    1295,    2400,     4095,     6560, ...
   31,  242,  1023,   3124,    7775,   16806,    32767,    59048, ...
   63,  728,  4095,  15624,   46655,  117648,   262143,   531440, ...
  127, 2186, 16383,  78124,  279935,  823542,  2097151,  4782968, ...
  255, 6560, 65535, 390624, 1679615, 5764800, 16777215, 43046720, ...
		

Crossrefs

Cf. A037205 (main diagonal).

Programs

  • Maple
    T:= (n, k)-> (k+1)^n-1:
    seq(seq(T(1+d-k, k), k=1..d), d=1..12);
  • Mathematica
    T[n_, k_] := (k + 1)^n - 1; Table[T[n - k + 1, k], {n, 1, 10}, {k, 1, n}] // Flatten (* Amiram Eldar, Nov 25 2019 *)
  • MuPAD
    T:=(n,k)->(k+1)^n-1:

Formula

T(n,k) = (k+1)^n - 1.

A355395 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals downwards, where T(n,k) = Sum_{j=0..n} k^(j*(n-j)) * binomial(n,j).

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 1, 2, 4, 2, 1, 2, 6, 8, 2, 1, 2, 8, 26, 16, 2, 1, 2, 10, 56, 162, 32, 2, 1, 2, 12, 98, 704, 1442, 64, 2, 1, 2, 14, 152, 2050, 15392, 18306, 128, 2, 1, 2, 16, 218, 4752, 84482, 593408, 330626, 256, 2, 1, 2, 18, 296, 9506, 318752, 7221250, 39691136, 8488962, 512, 2
Offset: 0

Views

Author

Seiichi Manyama, Jul 02 2022

Keywords

Comments

The Stanley reference below describes a family of binomial posets whose elements are two colored graphs with vertices labeled on [n] and with edges labeled on [k-1]. T(n,k) is the number of elements in an n-interval of such a binomial poset. - Geoffrey Critzer, Aug 21 2023

Examples

			Square array begins:
  1,  1,    1,     1,     1,      1, ...
  2,  2,    2,     2,     2,      2, ...
  2,  4,    6,     8,    10,     12, ...
  2,  8,   26,    56,    98,    152, ...
  2, 16,  162,   704,  2050,   4752, ...
  2, 32, 1442, 15392, 84482, 318752, ...
		

References

  • R. P. Stanley, Enumerative Combinatorics, Volume 1, Second Edition, Example 3.18.3(e), page 323.

Crossrefs

Columns k=0..4 give A040000, A000079, A047863, A135079, A355440.
Main diagonal gives A320287.
Cf. A009999.

Programs

  • PARI
    T(n, k) = sum(j=0, n, k^(j*(n-j))*binomial(n, j));

Formula

E.g.f. of column k: Sum_{j>=0} exp(k^j * x) * x^j/j!.
G.f. of column k: Sum_{j>=0} x^j/(1 - k^j * x)^(j+1).
For k>=1, E(x)^2 = Sum_{n>=0} T(n,k)*x^n/B_k(n) where B_k(n) = n!*k^binomial(n,2) and E(x) = Sum_{n>=0} x^n/b_k(n). - Geoffrey Critzer, Aug 21 2023
Showing 1-8 of 8 results.