OFFSET
0,2
COMMENTS
The e.g.f. A(x) of this sequence is motivated by the following identity:
Sum_{n>=0} p^n/n! * exp(q*r^n) = Sum_{n>=0} q^n/n! * exp(p*r^n) ;
here, p = x, q = 2*x, and r = exp(x).
FORMULA
E.g.f.: Sum_{n>=0} x^n * exp( 2*x*exp(n*x) ) / n!.
E.g.f.: Sum_{n>=0} 2^n*x^n * exp( x*exp(n*x) ) / n!.
EXAMPLE
E.g.f.: A(x) = 1 + 3*x + 9*x^2/2! + 39*x^3/3! + 249*x^4/4! + 2083*x^5/5! + 21309*x^6/6! + 258891*x^7/7! + 3676465*x^8/8! + 60188355*x^9/9! + ...
where
A(x) = exp(2*x) + x*exp(2*x*exp(x)) + x^2*exp(2*x*exp(2*x))/2! + x^3*exp(2*x*exp(3*x))/3! + x^4*exp(2*x*exp(4*x))/4! + x^5*exp(2*x*exp(5*x))/5! + ...
also
A(x) = exp(x) + 2*x*exp(x*exp(x)) + 2^2*x^2*exp(x*exp(2*x))/2! + 2^3*x^3*exp(x*exp(3*x))/3! + 2^4*x^4*exp(x*exp(4*x))/4! + 2^5*x^5*exp(x*exp(5*x))/5! + ...
PROG
(PARI) {a(n) = my(A=1); A = sum(m=0, n, (x^m/m!)*exp(2*x*exp(m*x +x*O(x^n)))); n!*polcoeff(A, n)}
for(n=0, 30, print1(a(n), ", "))
(PARI) {a(n) = my(A=1); A = sum(m=0, n, (2^m*x^m/m!)*exp(x*exp(m*x +x*O(x^n)))); n!*polcoeff(A, n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 30 2021
STATUS
approved