OFFSET
0,2
COMMENTS
Sum_{n>=0} a(n)/(3*n)! = exp( zeta(3) ) = 3.3269531100024997901915178...
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..150
EXAMPLE
E.g.f.: A(x) = 1 + 6*x^3/3! + 450*x^6/6! + 119280*x^9/9! + 78863400*x^12/12! +...
where
log(A(x)) = x^3 + x^6/2^3 + x^9/3^3 + x^12/4^3 + x^15/5^3 + x^18/6^3 +...
or,
log(A(x)) = 6*x^3/3! + 90*x^6/6! + 13440*x^9/9! + 7484400*x^12/12! + 10461394944*x^15/15! +...
MATHEMATICA
nmax=20; k=3; Table[(CoefficientList[Series[Exp[PolyLog[k, x^k]], {x, 0, k*nmax}], x] * Range[0, k*nmax]!)[[k*n-k+1]], {n, 1, nmax+1}] (* Vaclav Kotesovec, Jun 21 2015 *)
PROG
(PARI) {a(n) = (3*n)!*polcoeff( exp(sum(m=1, n, (x^m/m)^3)+x*O(x^(3*n))), 3*n)}
for(n=0, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 13 2015
STATUS
approved