OFFSET
0,2
COMMENTS
More generally, for m integer, exp( Sum_{n>=1} (m^n - 1)^n/(m-1)^(n-1) * x^n/n ) is a power series in x with integer coefficients.
Note that g.f. exp( Sum_{n>=1} (3^n - 1)^n/2^n * x^n/n ) has fractional coefficients as a power series in x.
EXAMPLE
G.f.: A(x) = 1 + 2*x + 18*x^2 + 1498*x^3 + 1283090*x^4 + 10377556482*x^5 +...
log(A(x)) = 2*x + 8^2/2*x^2/2 + 26^3/2^2*x^3/3 + 80^4/2^3*x^4/4 + 242^5/2^4*x^5/5 +...
PROG
(PARI) {a(n)=polcoeff(exp(sum(m=1, n+1, (3^m-1)^m/2^(m-1)*x^m/m)+x*O(x^n)), n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 04 2009
STATUS
approved