OFFSET
1,3
COMMENTS
Compare to C(x)^2 = C( x^2 + 2*x*C(x)^2 ) where C(x) = x + C(x)^2 is a g.f. of the Catalan numbers (A000108).
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..500
FORMULA
G.f. A(x) = Sum_{n>=1} a(n)*x^n satisfies the following formulas in which G(x) = x + (1/2)*(G(x)^2 + G(x^2)) is the g.f. of A001190, the Wedderburn-Etherington numbers.
(1) A(x)^2 = A( x^2 + 2*x*A(x)^2 + 2*A(x)^4 ).
(2) x = A( x - x^2 - x*G(x^2) ).
(3) x = A( x + x^2 - 2*x*G(x) + x*G(x)^2 ).
(4) x = A( x*sqrt(1 - 2*x^2 - G(x^4)) - x^2 ).
(5) x^2 = A( x^2*((1 - G(x))^2 + 2*x)^2 + x^4 ).
(6) G(A(x)) = 1 - sqrt(x/A(x) - A(x)).
EXAMPLE
G.f.: A(x) = x + x^2 + 3*x^3 + 10*x^4 + 39*x^5 + 161*x^6 + 699*x^7 + 3135*x^8 + 14427*x^9 + 67716*x^10 + 322959*x^11 + 1560585*x^12 + ...
where A(x)^2 = A( x^2 + 2*x*A(x)^2 + 2*A(x)^4 ).
RELATED SERIES.
Let G(x) be the g.f. of the Wedderburn-Etherington numbers, then
A( x - x^2 - x*G(x^2) ) = x, where G(x) = x + (1/2)*(G(x)^2 + G(x^2)) begins
G(x) = x + x^2 + x^3 + 2*x^4 + 3*x^5 + 6*x^6 + 11*x^7 + 23*x^8 + 46*x^9 + 98*x^10 + 207*x^11 + 451*x^12 + 983*x^13 + 2179*x^14 + ... + A001190(n)*x^n + ...
Let B(x) be the series reversion of g.f. A(x) so that B(A(x)) = x, then
B(x) = x - x^2 - x^3 - x^5 - x^7 - 2*x^9 - 3*x^11 - 6*x^13 - 11*x^15 - 23*x^17 - 46*x^19 - 98*x^21 - 207*x^23 + ...
where B(x) = x - x^2 - x*G(x^2).
SPECIFIC VALUES.
A(t) = 1/3 at t = 0.1804894059505127351310871071614416167035910065610113327...
where 1/9 = A( t^2 + 2*t/9 + 2/81 ).
A(t) = 1/4 at t = 0.1708289565101545485579649480920097855916395263217351536...
where 1/16 = A( t^2 + t/8 + 1/128 ).
A(t) = 1/5 at t = 0.1516661092515691718015998101146470241027491658579501286...
where 1/25 = A( t^2 + 2*t/25 + 2/625 ).
A(t) = 1/6 at t = 0.1341268789797555579297424694390747929782019601987848246...
where 1/36 = A( t^2 + t/18 + 1/648 ).
A(1/6) = 0.2368314953172156547771056118501694080205525703518284958...
A(1/7) = 0.1824082884402163049324182135107985537409785918465705698...
A(1/8) = 0.1515179821748020682616541846638756124979071552818869937...
A(1/9) = 0.1303577455916869424988611259176631850931169441135101392...
A(1/10) = 0.1146797533131163787803333792504789207692884367435306666...
PROG
(PARI) {a(n) = my(A=[0, 1], Ax=x); for(i=1, n, A = concat(A, 0); Ax=Ser(A);
A[#A] = (1/2)*polcoeff( subst(Ax, x, x^2 + 2*x*Ax^2 + 2*Ax^4) - Ax^2, #A) ); A[n+1]}
for(n=1, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 20 2024
STATUS
approved