OFFSET
0,2
COMMENTS
More generally, for m integer, exp( Sum_{n>=1} (m^n + y)^n * x^n/n ) is a power series in x and y with integer coefficients.
EXAMPLE
G.f.: A(x) = 1 + 5*x + 157*x^2 + 92285*x^3 + 1091087581*x^4 +...
log(A(x)) = 5*x + 17^2*x^2/2 + 65^3*x^3/3 + 257^4*x^4/4 + 1025^5*x^5/5 +...
PROG
(PARI) {a(n)=polcoeff(exp(sum(m=1, n+1, (4^m+1)^m*x^m/m)+x*O(x^n)), n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 04 2009
STATUS
approved