OFFSET
0,2
FORMULA
G.f. A(x) satisfies:
(1) A(x) = Sum_{n>=0} C(2^n,n)*x^n*A(x)^(2n),
(2) A(x) = sqrt((1/x)*Series_Reversion(x/B(x)^2)),
(3) A(x) = B(x*A(x)^2) and B(x) = A(x/B(x)^2),
where B(x) = Sum_{n>=0} C(2^n,n)*x^n is the g.f. of A014070.
(4) A(x) = F(x*A(x)) and F(x) = A(x/F(x)), where F(x) is the g.f. of A188193.
EXAMPLE
G.f.: A(x) = 1 + 2*x + 14*x^2 + 168*x^3 + 3756*x^4 + 261560*x^5 +...
which equals the series:
A(x) = 1 + log(1+2*x*A(x)^2) + log(1+4*x*A(x)^2)^2/2! + log(1+8*x*A(x)^2)^3/3! +...
Let B(x) equal the g.f. of A014070, which begins:
B(x) = 1 + 2*x + 6*x^2 + 56*x^3 + 1820*x^4 +...+ C(2^n,n)*x^n +...
then B(x) = A(x/B(x)^2) and A(x) = B(x*A(x)^2), so that:
A(x) = 1 + 2*x*A(x)^2 + 6*x^2*A(x)^4 + 56*x^3*A(x)^6 + 1820*x^4*A(x)^8 +...+ C(2^n,n)*x^n*A(x)^(2n) +...
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=sum(m=0, n, log(1+2^m*x*A^2+x*O(x^n))^m/m!)); polcoeff(A, n)}
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 23 2011
STATUS
approved