OFFSET
1,3
COMMENTS
Not the same as A003237.
FORMULA
G.f. satisfies the following identities:
(1) A(x^2) = A(x)^2 / (1 + 2*A(x) + 2*A(x)^2),
(2) A(-x) = -A(x) / (1 + 2*A(x)),
(3) A(x) + A(-x) = -2*A(x)*A(-x),
(4) A(x)^2 / (1 + 2*A(x)) = A(x^2) / (1 - 2*A(x^2)).
EXAMPLE
A(x) = x + x^2 + 2x^3 + 3x^4 + 6x^5 + 10x^6 + 19x^7 + 33x^8 + ... = x/(1-x -(x^2 + x^4 + x^6 + 2x^8 + 3x^10 + 11x^12 + 23x^14 + ...)).
PROG
(PARI) {a(n) = my(A=x, u, v); for(k=2, n, u=A+x*O(x^k); v=subst(u, x, x^2); A-=x^k*polcoeff(u^2 -v*(1+2*u+2*u^2), k+1)/2); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 23 2004
EXTENSIONS
Changed offset to 1 and removed leading zero. - Paul D. Hanna, Aug 16 2016
STATUS
approved