OFFSET
0,2
FORMULA
a(n) = Sum_{k=0..n} Stirling1(n, k)*Bell(3*k).
a(n) = n!*exp(-1)*Sum_{k>=[n^(1/3)]} binomial(k^3,n)/k!.
EXAMPLE
E.g.f.: A(x) = 1 + 5*x + 198*x^2/2! + 20548*x^3/3! + 4088918*x^4/4! +...
such that
A(x) = exp(-1)*(1 + (1+x) + (1+x)^8/2! + (1+x)^27/3! + (1+x)^64/4! +...).
PROG
(PARI) {Stirling1(n, k)=n!*polcoeff(binomial(x, n), k)}
{Bell(n)=n!*polcoeff(exp(exp(x+x*O(x^n))-1), n)}
{a(n)=sum(k=0, n, Stirling1(n, k)*Bell(3*k))}
for(n=0, 15, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 07 2012
STATUS
approved