OFFSET
0,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..200
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 12*x^3 + 27*x^4 + 105*x^5 + 420*x^6 +...
Related expansions:
A(x)^2 = 1 + 2*x + 7*x^2 + 30*x^3 + 87*x^4 + 336*x^5 + 1356*x^6 +...
A(x)^3 = 1 + 3*x + 12*x^2 + 55*x^3 + 189*x^4 + 756*x^5 + 3132*x^6 +...
Let w = exp(2*Pi*I/3), then A(x) = 1 + x*A(x)^3/(A(x)*A(w*x)*A(w^2*x)) where
A(x)*A(w*x)*A(w^2*x) = 1 + 28*x^3 + 1134*x^6 + 61857*x^9 + 3929121*x^12 + 272388420*x^15 + 19981576476*x^18 + 1524888581787*x^21 +...
PROG
(PARI) {a(n)=local(A=1+x*O(x^n)); for(i=1, n+1, A=1+x*A^3*exp(-3*sum(m=1, n\3, x^(3*m)*polcoeff(log(A), 3*m))+x*O(x^n))); polcoeff(A, n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 14 2012
STATUS
approved