cp's OEIS Frontend

A213009 G.f. A(x) satisfies: A(A(x)) = x+x^2 + x*A(A(A(A(x)))).

%I A213009 #8 Aug 03 2012 14:51:06
%S A213009 1,1,1,5,21,125,825,6133,49925,439417,4142945,41544161,440710117,
%T A213009 4924691541,57766255689,709205703565,9090541134373,121389729560633,
%U A213009 1685431945085489,24289856880005441,362776874949660485,5606980244843123077,89560387072919814553
%N A213009 G.f. A(x) satisfies: A(A(x)) = x+x^2 + x*A(A(A(A(x)))).
%C A213009 Given g.f. A(x), A(A(x)) equals the g.f. of A213010.
%H A213009 Paul D. Hanna, <a href="/A213009/b213009.txt">Table of n, a(n) for n = 1..256</a>
%F A213009 a(n) == 1 (mod 4).
%e A213009 G.f.: A(x) = x + x^2 + x^3 + 5*x^4 + 21*x^5 + 125*x^6 + 825*x^7 +...
%e A213009 where
%e A213009 A(A(x)) = x + 2*x^2 + 4*x^3 + 16*x^4 + 80*x^5 + 480*x^6 + 3296*x^7 +...
%e A213009 A(A(A(A(x)))) = x + 4*x^2 + 16*x^3 + 80*x^4 + 480*x^5 + 3296*x^6 +...
%o A213009 (PARI) {a(n)=local(A=x+x^2,B=x+2*x^2);for(i=1,n,B=x+x^2+x*subst(B,x,B+x*O(x^n)));
%o A213009 for(i=1,n,A=(A+subst(B,x,serreverse(A+x*O(x^n))))/2);polcoeff(A,n)}
%o A213009 for(n=1,31,print1(a(n),", "))
%Y A213009 Cf. A213010, A215115, A215117, A215119.
%K A213009 nonn
%O A213009 1,4
%A A213009 _Paul D. Hanna_, Jun 01 2012