OFFSET
0,2
COMMENTS
E.g.f. series identity: Sum_{n>=0} m^n * q^(n^2) * exp(b*q^n*x) * x^n/n! = Sum_{n>=0} (m*q^n + b)^n * x^n/n! for all q, m, b.
O.g.f. series identity: Sum_{n>=0} m^n * q^(n^2) * x^n/(1-b*q^n*x)^(n+1) = Sum_{n>=0} (m*q^n + b)^n * x^n for all q, m, b.
FORMULA
a(n) = (3^n + 1)^n.
O.g.f.: Sum_{n>=0} 3^(n^2) * x^n/(1 - 3^n*x)^(n+1).
EXAMPLE
E.g.f.: A(x) = 1 + 4*x + 100*x^2/2! + 21952*x^3/3! + 45212176*x^4/4! +..
By the series identity, the e.g.f.:
A(x) = exp(x) + 3*exp(3*x)*x + 3^4*exp(3^2*x)*x^2/2! + 3^9*exp(3^3*x)*x^3/3! +...
expands into:
A(x) = 1 + 4*x + 10^2*x^2/2! + 28^3*x^3/3! + 82^4*x^4/4! + 244^5*x^5/5! +...+ (3^n+1)^n*x^n/n! +...
PROG
(PARI) {a(n, q=3, m=1, b=1)=(m*q^n + b)^n}
(PARI) {a(n, q=3, m=1, b=1)=n!*polcoeff(sum(k=0, n, m^k*q^(k^2)*exp(b*q^k*x+x*O(x^n))*x^k/k!), n)}
(PARI) {a(n, q=3, m=1, b=1)=polcoeff(sum(k=0, n, m^k*q^(k^2)*x^k/(1-b*q^k*x +x*O(x^n))^(k+1)), n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 26 2011
STATUS
approved