A248471 E.g.f.: Sum_{n>=0} x^n * (3*exp(n*x) - 2)^n.

1, 1, 8, 87, 1836, 49515, 1888218, 89221755, 5280079752, 373918542867, 31364236227630, 3055178634627603, 341669063554284348, 43366411236710160651, 6190420280148645465330, 985839231828100191119595, 173934377807560257682643088, 33790653299040390831598566051, 7189209864606264796455405368022

Offset

0,3

Links

Table of n, a(n) for n=0..18.

Formula

E.g.f.: Sum_{n>=0} (3*x)^n * exp(n^2*x) / (1 + 2*x*exp(n*x))^(n+1).

Example

E.g.f.: A(x) = 1 + x + 8*x^2/2! + 87*x^3/3! + 1836*x^4/4! + 49515*x^5/5! +...
where the e.g.f. satisfies following series identity:
A(x) = 1 + x*(3*exp(x)-2) + x^2*(3*exp(2*x)-2)^2 + x^3*(3*exp(3*x)-2)^3 + x^4*(3*exp(4*x)-2)^4 + x^5*(3*exp(5*x)-2)^5 + x^6*(3*exp(6*x)-2)^6 +...
A(x) = 1/(1+2*x) + (3*x)*exp(x)/(1+2*x*exp(x))^2 + (3*x)^2*exp(4*x)/(1+2*x*exp(2*x))^3 + (3*x)^3*exp(9*x)/(1+2*x*exp(3*x))^4 + (3*x)^4*exp(16*x)/(1+2*x*exp(4*x))^5 + (3*x)^5*exp(25*x)/(1+2*x*exp(5*x))^6 + (3*x)^6*exp(36*x)/(1+2*x*exp(6*x))^7 +...

Prog

(PARI) {a(n)=local(A=1); A=sum(k=0, n, x^k * (3*exp(k*x +x*O(x^n)) - 2)^k); n!*polcoeff(A, n)}
for(n=0, 25, print1(a(n), ", "))
(PARI) {a(n)=local(A=1); A=sum(k=0, n, (3*x)^k * exp(k^2*x +x*O(x^n)) / (1 + 2*x*exp(k*x +x*O(x^n)))^(k+1) ); n!*polcoeff(A, n)}
for(n=0, 25, print1(a(n), ", "))

Crossrefs

Cf. A193421, A248615, A248653, A248654, A248655, A249409.

Sequence in context: A375542 A246512 A366233 * A284069 A292458 A292738

Adjacent sequences: A248468 A248469 A248470 * A248472 A248473 A248474

Keyword

nonn

Author

Paul D. Hanna, Oct 27 2014

Status

approved