OFFSET
1,3
EXAMPLE
G.f.: A(x) = x + x^2 + 2*x^3 + 5*x^4 + 13*x^5 + 36*x^6 + 104*x^7 +...
A(x*A(x)/(1-x)) = x^2 + 2*x^3 + 5*x^4 + 13*x^5 + 36*x^6 + 104*x^7 +...
Let G(x) = x*A(x)/(1-x) then
A(x) = x + G(x) + G(G(x)) + G(G(G(x))) + G(G(G(G(x)))) + ... where
G(x) = x^2 + 2*x^3 + 4*x^4 + 9*x^5 + 22*x^6 + 58*x^7 + 162*x^8 +...;
G(G(x)) = x^4 + 4*x^5 + 14*x^6 + 46*x^7 + 148*x^8 + 474*x^9 +...;
G(G(G(x))) = x^8 + 8*x^9 + 44*x^10 + 204*x^11 + 862*x^12 +...;
G(G(G(G(x)))) = x^16 + 16*x^17 + 152*x^18 + 1112*x^19 +...; ...
PROG
(PARI) {a(n)=local(A=x+x*O(x)); for(i=0, n, A=x+subst(A, x, x/(1-x)*A)); polcoeff(A, n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 16 2009
STATUS
approved