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

A167010 a(n) = Sum_{k=0..n} C(n,k)^n.

Original entry on oeis.org

1, 2, 6, 56, 1810, 206252, 86874564, 132282417920, 770670360699138, 16425660314368351892, 1367610300690018553312276, 419460465362069257397304825200, 509571049488109525160616367158261124, 2290638298071684282149128235413262383804352
Offset: 0

Views

Author

Paul D. Hanna, Nov 17 2009

Keywords

Comments

The number of n*n 0-1 matrices with equal numbers of nonzeros in every row. - David Eppstein, Jan 19 2012

Examples

			The triangle A209427 of coefficients C(n,k)^n, n>=k>=0, begins:
  1;
  1,     1;
  1,     4,        1;
  1,    27,       27,        1;
  1,   256,     1296,      256,        1;
  1,  3125,   100000,   100000,     3125,     1;
  1, 46656, 11390625, 64000000, 11390625, 46656,    1; ...
in which the row sums form this sequence.
		

Crossrefs

Programs

  • Magma
    [(&+[Binomial(n,j)^n: j in [0..n]]): n in [0..20]]; // G. C. Greubel, Aug 26 2022
    
  • Mathematica
    Table[Sum[Binomial[n, k]^n, {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Nov 05 2012 *)
  • PARI
    a(n)=sum(k=0,n,binomial(n,k)^n)
    
  • SageMath
    [sum(binomial(n,j)^n for j in (0..n)) for n in (0..20)] # G. C. Greubel, Aug 26 2022

Formula

Ignoring initial term, equals the logarithmic derivative of A167007. [Paul D. Hanna, Nov 18 2009]
If n is even then a(n) ~ c * exp(-1/4) * 2^(n^2 + n/2)/((Pi*n)^(n/2)), where c = Sum_{k = -oo..oo} exp(-2*k^2) = 1.271341522189... (see A218792). - Vaclav Kotesovec, Nov 05 2012
If n is odd then c = Sum_{k = -infinity..infinity} exp(-2*(k+1/2)^2) = 1.23528676585389... - Vaclav Kotesovec, Nov 06 2012
a(n) = (n!)^n * [x^n] (Sum_{k>=0} x^k / (k!)^n)^2. - Ilya Gutkovskiy, Jul 15 2020

A209424 Triangle defined by g.f.: A(x,y) = exp( Sum_{n>=1} x^n/n * Sum_{k=0..n} binomial(n, k)^n * y^k ), as read by rows.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 12, 12, 1, 1, 76, 347, 76, 1, 1, 701, 20429, 20429, 701, 1, 1, 8477, 1919660, 10707908, 1919660, 8477, 1, 1, 126126, 259227625, 9203978774, 9203978774, 259227625, 126126, 1, 1, 2223278, 47484618291, 12099129236936, 72078431500368
Offset: 0

Views

Author

Paul D. Hanna, Mar 08 2012

Keywords

Comments

Column 1 is A060946.
Column 2 is A209425.
Row sums equal A167007.
Antidiagonal sums equal A166894.
Central terms form A209426.

Examples

			This triangle begins:
1;
1, 1;
1, 3, 1;
1, 12, 12, 1;
1, 76, 347, 76, 1;
1, 701, 20429, 20429, 701, 1;
1, 8477, 1919660, 10707908, 1919660, 8477, 1;
1, 126126, 259227625, 9203978774, 9203978774, 259227625, 126126, 1;
1, 2223278, 47484618291, 12099129236936, 72078431500368, 12099129236936, 47484618291, 2223278, 1; ...
G.f.: A(x,y) = 1 + (1+y)*x + (1+3*y+y^2)*x^2 + (1+12*y+12*y^2+y^3)*x^3 + (1+76*y+20429*y^2+76*y^3+y^4)*x^4 +...
The logarithm of the g.f. equals the series:
log(A(x,y)) = (1 + y)*x
+ (1 + 2^2*y + y^2)*x^2/2
+ (1 + 3^3*y + 3^3*y^2 + y^3)*x^3/3
+ (1 + 4^4*y + 6^4*y^2 + 4^4*y^3 + y^4)*x^4/4
+ (1 + 5^5*y + 10^5*y^2 + 10^5*y^3 + 5^5*y^4 + y^5)*x^5/5 +...
in which the coefficients are found in triangle A209427.
		

Crossrefs

Programs

  • PARI
    {T(n,k)=polcoeff(polcoeff(exp(sum(m=1,n,x^m/m*sum(k=0,m,binomial(m,k)^m*y^k))+x*O(x^n)),n,x),k,y)}
    for(n=0,10,for(k=0,n,print1(T(n,k),", "));print(""))
Showing 1-2 of 2 results.