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 12 results. Next

A007820 Stirling numbers of second kind S(2n,n).

Original entry on oeis.org

1, 1, 7, 90, 1701, 42525, 1323652, 49329280, 2141764053, 106175395755, 5917584964655, 366282500870286, 24930204590758260, 1850568574253550060, 148782988064375309400, 12879868072770626040000, 1194461517469807833782085, 118144018577011378596484455
Offset: 0

Views

Author

kemp(AT)sads.informatik.uni-frankfurt.de (Rainer Kemp)

Keywords

Comments

Chan and Manna prove that a(n) is odd if and only if n is in A003714. - Jason Kimberley, Sep 14 2009
The number of ways to partition a set of 2*n elements into n disjoint subsets. - Vladimir Reshetnikov, Oct 10 2016
Conjecture: a(2*n+1) is divisible by (2*n + 1)^2. - Peter Bala, Mar 30 2025

Examples

			G.f.: A(x) = 1 + x + 7*x^2 + 90*x^3 + 1701*x^4 + 42525*x^5 +...,
where A(x) = 1 + 1^2*x*exp(-1*x) + 2^4*exp(-2^2*x)*x^2/2! + 3^6*exp(-3^2*x)*x^3/3! + 4^8*exp(-4^2*x)*x^4/4! + 5^10*exp(-5^2*x)*x^5/5! + ... - _Paul D. Hanna_, Oct 17 2012
		

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. 835.

Crossrefs

Programs

  • Maple
    A007820 := proc(n) Stirling2(2*n,n) ; end proc:
    seq(A007820(n),n=0..20) ; # R. J. Mathar, Mar 15 2011
  • Mathematica
    Table[StirlingS2[2n, n], {n, 1, 12}] (* Emanuele Munarini, Mar 12 2011 *)
  • Maxima
    makelist(stirling2(2*n,n),n,0,12); /* Emanuele Munarini, Mar 12 2011 */
    
  • PARI
    a(n)=stirling(2*n,n,2); /* Joerg Arndt, Jul 01 2011 */
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(2*n))), n)} \\ Paul D. Hanna, Oct 17 2012
    
  • PARI
    {a(n)=polcoeff(sum(m=1,n,(m^2)^m*exp(-m^2*x+x*O(x^n))*x^m/m!),n)} \\ Paul D. Hanna, Oct 17 2012
    
  • Python
    from sympy.functions.combinatorial.numbers import stirling
    def A007820(n): return stirling(n<<1,n) # Chai Wah Wu, Jun 09 2025
  • Sage
    [stirling_number2(2*i,i) for i in range(1,20)] # Zerinvary Lajos, Jun 26 2008
    

Formula

a(n) = A048993(2n,n). - R. J. Mathar, Mar 15 2011
Asymptotic: a(n) ~ (4*n/(e*z*(2-z)))^n/sqrt(2*Pi*n*(z-1)), where z = A256500 = 1.59362426... is a root of the equation exp(z)*(2-z)=2. - Vaclav Kotesovec, May 30 2011
a(n) = 1/n! * Sum_{k = 0..n} binomial(n,k)*(-1)^k*(n-k)^(2*n). - Emanuele Munarini, Jul 01 2011
a(n) = [x^n] 1 / Product_{k=1..n} (1-k*x). - Paul D. Hanna, Oct 17 2012
O.g.f.: Sum_{n>=1} (n^2)^n * exp(-n^2*x) * x^n/n! = Sum_{n>=1} S2(2*n,n)*x^n. - Paul D. Hanna, Oct 17 2012
G.f.: Sum_{n > 0} (a(n)*n!/(2*n)!)*x^n = x*B'(x)/B(x)-1, where B(x) satisfies B(x)^2 = x*(exp(B(x))-1). - Vladimir Kruchinin, Mar 13 2013
a(n) = Sum_{j = 0..n} (-1)^(n-j)*n^j*binomial(2*n,j)*stirling2(2*n-j,n). - Vladimir Kruchinin, Jun 14 2013

Extensions

Typo in Mathematica program fixed by Vincenzo Librandi, May 04 2013
a(0)=1 prepended by Alois P. Heinz, Feb 01 2018

A217913 O.g.f.: Sum_{n>=0} (n^3)^n * exp(-n^3*x) * x^n / n!.

Original entry on oeis.org

1, 1, 31, 3025, 611501, 210766920, 110687251039, 82310957214948, 82318282158320505, 106563273280541795575, 173373343599189364594756, 346289681454731077633095526, 833091176987705031151553054843, 2376102520162485084539597049185710
Offset: 0

Views

Author

Paul D. Hanna, Oct 14 2012

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 31*x^2 + 3025*x^3 + 611501*x^4 + ... + Stirling2(3*n, n)*x^n + ...
where
A(x) = 1 + 1^3*x*exp(-1^3*x) + 2^6*exp(-2^3*x)*x^2/2! + 3^9*exp(-3^3*x)*x^3/3! + 4^12*exp(-4^3*x)*x^4/4! + 5^15*exp(-5^3*x)*x^5/5! + ...
simplifies to a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[3*n,n],{n,0,20}] (* Vaclav Kotesovec, Feb 28 2013 *)
  • Maxima
    makelist(stirling2(3*n, n), n, 0, 13); /* Martin Ettl, Oct 15 2012 */
  • PARI
    {a(n)=polcoeff(sum(k=0,n,(k^3)^k*exp(-k^3*x +x*O(x^n))*x^k/k!),n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^3)^k*x^k/(1+k^3*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(2*n))), 2*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(3*n, n)}
    for(n=0,20,print1(a(n),", "))
    

Formula

a(n) = Stirling2(3*n, n).
a(n) = [x^(3*n)] (3*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(2*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^3)^k*x^k / (1 + k^3*x)^(k+1).
a(n) ~ 9^n*exp(n*(c+1))*n^(2*n)/((c+3)^(2*n)*sqrt(2*Pi*(c+1)*n)), where c = -0.1785606278779211... = LambertW(-3/exp(3)) = A226750. - Vaclav Kotesovec, Feb 28 2013

A217915 O.g.f.: Sum_{n>=1} (n^5)^n * exp(-n^5*x) * x^n / n!.

Original entry on oeis.org

1, 1, 511, 2375101, 45232115901, 2436684974110751, 299310102746948685757, 72786959006434393367186463, 31712979422428631132831124895809, 22982258052528294182955639980819773510, 26154716515862881292012777396577993781727011
Offset: 0

Views

Author

Paul D. Hanna, Oct 14 2012

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 511*x^2 + 2375101*x^3 + 45232115901*x^4 +...+ Stirling2(5*n, n)*x^n +...
where
A(x) = 1 + 1^5*x*exp(-1^5*x) + 2^10*exp(-2^5*x)*x^2/2! + 3^15*exp(-3^5*x)*x^3/3! + 4^20*exp(-4^5*x)*x^4/4! + 5^25*exp(-5^5*x)*x^5/5! +...
is a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[5*n,n],{n,0,20}] (* Vaclav Kotesovec, May 23 2013 *)
  • Maxima
    makelist(stirling2(5*n, n), n, 0, 10); /* Martin Ettl, Oct 15 2012 */
  • PARI
    {a(n)=polcoeff(sum(k=0,n,(k^5)^k*exp(-k^5*x +x*O(x^n))*x^k/k!),n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^5)^k*x^k/(1+k^5*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(4*n))), 4*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(5*n, n)}
    for(n=0,12,print1(a(n),", "))
    

Formula

a(n) = Stirling2(5*n, n).
a(n) = [x^(5*n)] (5*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(4*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^5)^k*x^k / (1 + k^5*x)^(k+1).
a(n) ~ n^(4*n)*5^(5*n) / (sqrt(2*Pi*n*(1-c)) * exp(4*n) * (5-c)^(4*n) * c^n), where c = -LambertW(-5/exp(5)) = 0.0348857682557... - Vaclav Kotesovec, May 23 2013

A218141 a(n) = Stirling2(n^2, n).

Original entry on oeis.org

1, 1, 7, 3025, 171798901, 2436684974110751, 14204422416132896951197888, 50789872166903636182659702516635946082, 155440114706926165785630654089245708839702615196926765, 541500903058656141876322139677626107784896646583041951351456223689104719
Offset: 0

Views

Author

Paul D. Hanna, Oct 21 2012

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 7*x^2 + 3025*x^3 + 171798901*x^4 + 2436684974110751*x^5 +...
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[n^2, n],{n,0,10}] (* Vaclav Kotesovec, May 11 2014 *)
  • Maxima
    makelist(stirling2(n^2,n),n,0,30 ); /* Martin Ettl, Oct 21 2012 */
  • PARI
    {a(n)=polcoeff(sum(k=0,n,(k^n)^k*exp(-k^n*x +x*O(x^n))*x^k/k!),n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(n^2+1))), n^2-n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(n^2, n)}
    for(n=0, 10, print1(a(n), ", "))
    

Formula

a(n) = [x^n] Sum_{k>=0} k^(n*k) * exp(-k^n*x) * x^k / k!.
a(n) = [x^(n^2-n)] 1 / Product_{k=1..n} (1-k*x).
a(n) ~ n^(n^2)/n!. - Vaclav Kotesovec, May 11 2014

A222526 O.g.f.: Sum_{n>=0} (n^6)^n * exp(-n^6*x) * x^n / n!.

Original entry on oeis.org

1, 1, 2047, 64439010, 11681056634501, 7713000216608565075, 14204422416132896951197888, 61232072982330045410678351728440, 545827051514425992551826008968173372261, 9173647538352903119028122246836507680995590680
Offset: 0

Views

Author

Paul D. Hanna, Feb 23 2013

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 2047*x^2 + 64439010*x^3 + 11681056634501*x^4 +...+ Stirling2(6*n, n)*x^n +...
where
A(x) = 1 + 1^6*x*exp(-1^6*x) + 2^12*exp(-2^6*x)*x^2/2! + 3^18*exp(-3^6*x)*x^3/3! + 4^24*exp(-4^6*x)*x^4/4! + 5^30*exp(-5^6*x)*x^5/5! +...
is a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[6*n, n],{n,0,20}] (* Vaclav Kotesovec, May 11 2014 *)
  • PARI
    {a(n)=polcoeff(sum(k=0, n, (k^6)^k*exp(-k^6*x +x*O(x^n))*x^k/k!), n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^6)^k*x^k/(1+k^6*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(5*n))), 5*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(6*n, n)}
    for(n=0, 12, print1(a(n), ", "))

Formula

a(n) = Stirling2(6*n, n).
a(n) = [x^(6*n)] (6*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(5*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^6)^k*x^k / (1 + k^6*x)^(k+1).
a(n) ~ n^(5*n) * 6^(6*n) / (sqrt(2*Pi*(1-c)*n) * exp(5*n) * (6-c)^(5*n) * c^n), where c = -LambertW(-6*exp(-6)). - Vaclav Kotesovec, May 11 2014

A222527 O.g.f.: Sum_{n>=0} (n^7)^n * exp(-n^7*x) * x^n / n!.

Original entry on oeis.org

1, 1, 8191, 1742343625, 2998587019946701, 24204004899040755811870, 666480349285726891499539272955, 50789872166903636182659702516635946082, 9237419992097529135737293866043969707761346313, 3590622358224471993651445012122431990834934483552661750
Offset: 0

Views

Author

Paul D. Hanna, Feb 23 2013

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 8191*x^2 + 1742343625*x^3 + 2998587019946701*x^4 +...+ Stirling2(7*n, n)*x^n +...
where
A(x) = 1 + 1^7*x*exp(-1^7*x) + 2^14*exp(-2^7*x)*x^2/2! + 3^21*exp(-3^7*x)*x^3/3! + 4^28*exp(-4^7*x)*x^4/4! + 5^35*exp(-5^7*x)*x^5/5! +...
is a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[7*n, n],{n,0,20}] (* Vaclav Kotesovec, May 11 2014 *)
  • PARI
    {a(n)=polcoeff(sum(k=0, n, (k^7)^k*exp(-k^7*x +x*O(x^n))*x^k/k!), n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^7)^k*x^k/(1+k^7*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(6*n))), 6*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(7*n, n)}
    for(n=0, 12, print1(a(n), ", "))

Formula

a(n) = Stirling2(7*n, n).
a(n) = [x^(7*n)] (7*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(6*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^7)^k*x^k / (1 + k^7*x)^(k+1).
a(n) ~ n^(6*n) * 7^(7*n) / (sqrt(2*Pi*(1-c)*n) * exp(6*n) * (7-c)^(6*n) * c^n), where c = -LambertW(-7*exp(-7)). - Vaclav Kotesovec, May 11 2014

A222528 O.g.f.: Sum_{n>=0} (n^8)^n * exp(-n^8*x) * x^n / n!.

Original entry on oeis.org

1, 1, 32767, 47063200806, 768305500780164501, 75740854251732106906082250, 31154086963475828638359480518580526, 41929298560838945526242744414099901692285884, 155440114706926165785630654089245708839702615196926765, 1396002062838446082394548660243302585983358463911636390911298400
Offset: 0

Views

Author

Paul D. Hanna, Feb 23 2013

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 32767*x^2 + 47063200806*x^3 + 768305500780164501*x^4 +...+ Stirling2(8*n, n)*x^n +...
where
A(x) = 1 + 1^8*x*exp(-1^8*x) + 2^16*exp(-2^8*x)*x^2/2! + 3^24*exp(-3^8*x)*x^3/3! + 4^32*exp(-4^8*x)*x^4/4! + 5^40*exp(-5^8*x)*x^5/5! +...
is a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[8*n, n],{n,0,20}] (* Vaclav Kotesovec, May 11 2014 *)
  • PARI
    {a(n)=polcoeff(sum(k=0, n, (k^8)^k*exp(-k^8*x +x*O(x^n))*x^k/k!), n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^8)^k*x^k/(1+k^8*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(7*n))), 7*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(8*n, n)}
    for(n=0, 12, print1(a(n), ", "))

Formula

a(n) = Stirling2(8*n, n).
a(n) = [x^(8*n)] (8*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(7*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^8)^k*x^k / (1 + k^8*x)^(k+1).
a(n) ~ n^(7*n) * 8^(8*n) / (sqrt(2*Pi*(1-c)*n) * exp(7*n) * (8-c)^(7*n) * c^n), where c = -LambertW(-8*exp(-8)). - Vaclav Kotesovec, May 11 2014

A222529 O.g.f.: Sum_{n>=0} (n^9)^n * exp(-n^9*x) * x^n / n!.

Original entry on oeis.org

1, 1, 131071, 1270865805301, 196740254364198919901, 236795997997922560392792426501, 1454443713270449746545892977574122129433, 34559048315358253352594346952765431711799794270765, 2610516895723221966171633379256064857587637240616032299710417
Offset: 0

Views

Author

Paul D. Hanna, Feb 23 2013

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 131071*x^2 + 1270865805301*x^3 + 196740254364198919901*x^4 +...+ Stirling2(9*n, n)*x^n +...
where
A(x) = 1 + 1^9*x*exp(-1^9*x) + 2^18*exp(-2^9*x)*x^2/2! + 3^27*exp(-3^9*x)*x^3/3! + 4^36*exp(-4^9*x)*x^4/4! + 5^45*exp(-5^9*x)*x^5/5! +...
is a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[9*n, n],{n,0,20}] (* Vaclav Kotesovec, May 11 2014 *)
  • PARI
    {a(n)=polcoeff(sum(k=0, n, (k^9)^k*exp(-k^9*x +x*O(x^n))*x^k/k!), n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^9)^k*x^k/(1+k^9*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(8*n))), 8*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(9*n, n)}
    for(n=0, 12, print1(a(n), ", "))

Formula

a(n) = Stirling2(9*n, n).
a(n) = [x^(9*n)] (9*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(8*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^9)^k*x^k / (1 + k^9*x)^(k+1).
a(n) ~ n^(8*n) * 9^(9*n) / (sqrt(2*Pi*(1-c)*n) * exp(8*n) * (9-c)^(8*n) * c^n), where c = -LambertW(-9*exp(-9)). - Vaclav Kotesovec, May 11 2014

A222530 O.g.f.: Sum_{n>=1} (n^10)^n * exp(-n^10*x) * x^n / n!.

Original entry on oeis.org

1, 1, 524287, 34314651811530, 50369882873307917364901, 740095864368253016271188139587625, 67872880319721869662486234870635119906757244, 28468832412072117193931250482560479429446507352468258480, 43812568949824405485262661429905291482204531455805230631187460302069
Offset: 0

Views

Author

Paul D. Hanna, Feb 23 2013

Keywords

Comments

Generally, for p>=2 is StirlingS2(p*n,n) asymptotic to n^((p-1)*n) * p^(p*n) / (sqrt(2*Pi*(1-c)*n) * exp((p-1)*n) * (p-c)^((p-1)*n) * c^n), where c = -LambertW(-p*exp(-p)). - Vaclav Kotesovec, May 11 2014

Examples

			O.g.f.: A(x) = 1 + x + 524287*x^2 + 34314651811530*x^3 + 50369882873307917364901*x^4 +...+ Stirling2(10*n, n)*x^n +...
where
A(x) = 1 + 1^10*x*exp(-1^10*x) + 2^20*exp(-2^10*x)*x^2/2! + 3^30*exp(-3^10*x)*x^3/3! + 4^40*exp(-4^10*x)*x^4/4! + 5^50*exp(-5^10*x)*x^5/5! +...
is a power series in x with integer coefficients.
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[10*n, n],{n,0,20}] (* Vaclav Kotesovec, May 11 2014 *)
  • PARI
    {a(n)=polcoeff(sum(k=0, n, (k^10)^k*exp(-k^10*x +x*O(x^n))*x^k/k!), n)}
    
  • PARI
    {a(n)=1/n!*polcoeff(sum(k=0, n, (k^1)^k*x^k/(1+k^10*x +x*O(x^n))^(k+1)), n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(9*n))), 9*n)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(10*n, n)}
    for(n=0, 12, print1(a(n), ", "))

Formula

a(n) = Stirling2(10*n, n).
a(n) = [x^(10*n)] (10*n)! * (exp(x) - 1)^n / n!.
a(n) = [x^(9*n)] 1 / Product_{k=1..n} (1-k*x).
a(n) = 1/n! * [x^n] Sum_{k>=0} (k^10)^k*x^k / (1 + k^10*x)^(k+1).
a(n) ~ n^(9*n) * 10^(10*n) / (sqrt(2*Pi*(1-c)*n) * exp(9*n) * (10-c)^(9*n) * c^n), where c = -LambertW(-10*exp(-10)). - Vaclav Kotesovec, May 11 2014

A218142 a(n) = Stirling2(n^2+n, n).

Original entry on oeis.org

1, 1, 31, 86526, 45232115901, 7713000216608565075, 666480349285726891499539272955, 41929298560838945526242744414099901692285884, 2610516895723221966171633379256064857587637240616032299710417
Offset: 0

Views

Author

Paul D. Hanna, Oct 21 2012

Keywords

Examples

			O.g.f.: A(x) = 1 + x + 31*x^2 + 86526*x^3 + 45232115901*x^4 +...
		

Crossrefs

Programs

  • Mathematica
    Table[StirlingS2[n^2+n, n],{n,0,10}] (* Vaclav Kotesovec, May 11 2014 *)
  • Maxima
    makelist(stirling2(n^2+n,n),n,0,30 ); /* Martin Ettl, Oct 21 2012 */
  • PARI
    {a(n)=polcoeff(sum(k=0,n,(k^(n+1))^k*exp(-k^(n+1)*x +x*O(x^n))*x^k/k!),n)}
    
  • PARI
    {a(n)=polcoeff(1/prod(k=1, n, 1-k*x +x*O(x^(n^2))), n^2)}
    
  • PARI
    {Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
    {a(n) = Stirling2(n^2+n, n)}
    for(n=0, 10, print1(a(n), ", "))
    

Formula

a(n) = [x^n] Sum_{k>=0} k^((n+1)*k) * exp(-k^(n+1)*x) * x^k / k!.
a(n) = [x^(n^2)] 1 / Product_{k=1..n} (1-k*x).
a(n) ~ n^(n^2+n)/n!. - Vaclav Kotesovec, May 11 2014
Showing 1-10 of 12 results. Next