OFFSET
0,2
FORMULA
a(n) = Sum_{k=0..n} Stirling1(n, k)*Bell(4*k).
a(n) = n!*exp(-1)*Sum_{k>=[n^(1/4)]} binomial(k^4,n)/k!.
EXAMPLE
E.g.f.: A(x) = 1 + 15*x + 4125*x^2/2! + 4201207*x^3/3! + 10454906015*x^4/4! +...
such that
A(x) = exp(-1)*(1 + (1+x) + (1+x)^16/2! + (1+x)^81/3! + (1+x)^256/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(4*k))}
for(n=0, 15, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 07 2012
STATUS
approved