OFFSET
1,2
FORMULA
a(n) = [x^n] F_n(x) where F_n(x) = F_{n-1}(F_1(x)) with F_1(x) = x*(1+x)^2.
EXAMPLE
The successive iterations of F(x) = x*(1+x)^2 begin:
F(x) = (1)x + 2x^2 + x^3
F(F(x)) = x + (4)x^2 + 10x^3 + 18x^4 + 23x^5 + 22x^6 + 15x^7 + 6x^8 +...
F(F(F(x))) = x + 6x^2 + (27)x^3 + 102x^4 + 333x^5 + 960x^6 + 2472x^7 +...
F(F(F(F(x)))) = x + 8x^2 + 52x^3 + (300)x^4 + 1578x^5 + 7692x^6 +...
F(F(F(F(F(x))))) = x + 10x^2 + 85x^3 + 660x^4 + (4790)x^5 + 32920x^6+...
F(F(F(F(F(F(x)))))) = x + 12x^2 +126x^3 +1230x^4+11385x^5+(101010)x^6+...
PROG
(PARI) {a(n)=local(F=x*(1+x)^2, G=x+x*O(x^n)); if(n<1, 0, for(i=1, n, G=subst(F, x, G)); return(polcoeff(G, n, x)))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 01 2006
STATUS
approved