OFFSET
1,4
COMMENTS
Radius of convergence is r = sqrt(2)-1, where r = r^2*(1+r)/(1-r), with A(r) = 1.
Compare to a g.f. M(x) of Motzkin numbers: M(x)^2 = M(x^2/(1-2*x)) where M(x) = (1-x - sqrt(1-2*x-3*x^2))/(2*x).
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..520
EXAMPLE
G.f.: A(x) = x + x^2 + x^3 + 2*x^4 + 3*x^5 + 5*x^6 + 9*x^7 + 18*x^8 + 38*x^9 + 80*x^10 + 166*x^11 + 342*x^12 + 705*x^13 + 1463*x^14 + 3065*x^15 + 6492*x^16 +...
such that A(x^2*(1+x)/(1-x)) = A(x)^2, where
A(x)^2 = x^2 + 2*x^3 + 3*x^4 + 6*x^5 + 11*x^6 + 20*x^7 + 38*x^8 + 76*x^9 + 159*x^10 + 338*x^11 + 719*x^12 + 1526*x^13 + 3235*x^14 + 6868*x^15 + 14638*x^16 +...
RELATED SERIES.
Let B(x) be the series reversion of g.f. A(x), so that A(B(x)) = x, then
B(x) = x - x^2 + x^3 - 2*x^4 + 5*x^5 - 12*x^6 + 29*x^7 - 75*x^8 + 201*x^9 - 542*x^10 + 1480*x^11 - 4112*x^12 + 11562*x^13 - 32770*x^14 + 93593*x^15 +...
where
B(x^2) = B(x)^2*(1 + B(x))/(1 - B(x)).
PROG
(PARI) {a(n) = my(A=x); for(i=1, n, A = sqrt( subst(A, x, x^2*(1+x)/(1-x +x*O(x^n))) ) ); polcoeff(A, n)}
for(n=1, 50, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 08 2015
STATUS
approved