OFFSET
1,3
EXAMPLE
E.g.f.: A(x) = x + x^2/2! + 8*x^3/3! + 144*x^4/4! + 4554*x^5/5! + 230940*x^6/6! + 18177900*x^7/7! + 2196712980*x^8/8! + 406693854000*x^9/9! + 115319921466960*x^10/10! + ...
such that
Sum_{n>=0} A(x)^n * exp(2^n*A(x)) / n! = exp(A(x)) + A(x)*exp(2*A(x)) + A(x)^2*exp(2^2*A(x))/2! + A(x)^3*exp(2^3*A(x))/3! + A(x)^4*exp(2^4*A(x))/4! +...
equals the sum
Sum_{n>=0} x^n/(1 - 2^n*x) = 1 + 2*x + 4*x^2 + 10*x^3 + 34*x^4 + 162*x^5 + 1090*x^6 + 10370*x^7 + 139522*x^8 + ... + A117402(n)*x^n + ...
RELATED SERIES.
exp(A(x)) = 1 + x + 2*x^2/2! + 12*x^3/3! + 186*x^4/4! + 5460*x^5/5! + 263940*x^6/6! + 20053740*x^7/7! + 2359326480*x^8/8! + 428122913400*x^9/9! + ...
PROG
(PARI) {a(n) = my(L=[0, 1]); for(i=1, n, L=concat(L, 0);
L[#L] = polcoeff( sum(n=0, #L, x^n/(1 - 2^n*x +x*O(x^#L))) - sum(n=0, #L, Ser(L)^n/n! * exp(2^n*Ser(L)) ) , #L-1)/2; ); n!*L[n+1]}
for(n=1, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 09 2021
STATUS
approved