OFFSET
0,2
EXAMPLE
A(x) = 1 + 3*x + 3*x^2 + x^3 + 6*x^4 + 6*x^5 + 9*x^6 +...
A(x)^3 = 1 + 9*x + 36*x^2 + 84*x^3 + 144*x^4 + 252*x^5 +...
A(x)^3 (mod 27) = 1 + 9*x + 9*x^2 + 3*x^3 + 9*x^4 + 9*x^5 +...
G(x) = 1 + 9*x + 9*x^2 + 3*x^3 + 9*x^4 + 9*x^5 + 3*x^6 +...
where G(x) is the g.f. of A083949.
PROG
(PARI) {a(n)=local(d=3, m=9, A=1+m*x); for(j=2, d*n, for(k=1, m, t=polcoeff((A+k*x^j+x*O(x^j))^(1/m), j); if(denominator(t)==1, A=A+k*x^j; break))); polcoeff(A, d*n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v and Paul D. Hanna, Aug 30 2005
STATUS
approved