OFFSET
0,2
COMMENTS
More generally, the following sums are equal:
(1) Sum_{n>=0} (q^n + p)^n * r^n/n!,
(2) Sum_{n>=0} q^(n^2) * exp(p*q^n*r) * r^n/n!;
here, q = (1+x) and p = -1, r = 4.
In general, let F(x) be a formal power series in x such that F(0)=1, then
Sum_{n>=0} m^n * F(q^n*r)^p * log( F(q^n*r) )^n / n! =
Sum_{n>=0} r^n * [y^n] F(y)^(m*q^n + p);
here, F(x) = exp(x), q = 1+x, p = -1, r = 4, m = 1.
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..300
FORMULA
E.g.f. may be expressed by the following sums.
(1) Sum_{n>=0} ((1+x)^n - 1)^n * 4^n / n!.
(2) Sum_{n>=0} (1+x)^(n^2) * exp(-4*(1+x)^n) * 4^n / n!.
EXAMPLE
E.g.f: A(x) = 1 + 4*x + 64*x^2/2! + 1920*x^3/3! + 86464*x^4/4! + 5304320*x^5/5! + 418131456*x^6/6! + 40727959552*x^7/7! + 4765747597312*x^8/8! + 655794545577984*x^9/9! + 104360850604687360*x^10/10! +...
such that
A(x) = 1 + 4*((1+x) - 1) + 4^2*((1+x)^2 - 1)^2/2! + 4^3*((1+x)^3 - 1)^3/3! + 4^4*((1+x)^4 - 1)^4/4! + 4^5*((1+x)^5 - 1)^5/5! + 4^6*((1+x)^6 - 1)^6/6! + 4^7*((1+x)^7 - 1)^7/7! + ...
also
A(x) = 1 + 4*(1+x)*exp(-4*(1+x)) + 4^2*(1+x)^4*exp(-4*(1+x)^2)/2! + 4^3*(1+x)^9*exp(-4*(1+x)^3)/3! + 4^4*(1+x)^16*exp(-4*(1+x)^4)/4! + 4^5*(1+x)^25*exp(-4*(1+x)^5)/5! + 4^6*(1+x)^36*exp(-4*(1+x)^6)/6! + 4^7*(1+x)^49*exp(-4*(1+x)^7)/7! + ...
PROG
(PARI) {a(n)=n!*polcoeff(sum(m=0, n, 4^m*((1+x+x*O(x^n))^m-1)^m/m!), n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 22 2019
STATUS
approved