OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..500
FORMULA
G.f. A(x) satisfies [Paul D. Hanna, Nov 24 2012]:
(1) A(x) = (1/x)*series_reversion(x*eta(x)^2).
(2) A(x) = 1 / Product_{n>=1} (1 - x^n*A(x)^(2*n)).
(3) A(x) = Sum_{n>=0} x^n*A(x)^(2*n) / Product_{k=1..n} (1-x^k*A(x)^(2*k)).
(4) A(x) = Sum_{n>=0} (x*A(x)^2)^(n^2) / Product_{k=1..n} (1-x^k*A(x)^(2*k))^2.
(5) A(x) = exp( Sum_{n>=1} (x^n/n) * A(x)^(2*n)/(1 - x^n*A(x)^(2*n)) ).
a(n) ~ c * d^n / n^(3/2), where d = 8.42516721063251541777601555584151410936... and c = 0.2128745515668564974075326286129891378270... - Vaclav Kotesovec, May 13 2018
EXAMPLE
G.f.: A(x) = 1 + x + 4*x^2 + 20*x^3 + 115*x^4 + 714*x^5 +...
G.f. satisfies A(x/P(x)^2) = P(x) where:
P(x) = 1 + x + 2*x^2 + 3*x^3 + 5*x^4 + 7*x^5 + 11*x^6 + 15*x^7 +...
and x/P(x)^2 = x - 2*x^2 - x^3 + 2*x^4 + x^5 + 2*x^6 - 2*x^7 - 2*x^9 +...
MATHEMATICA
nmax = 20; A[_] = 0; Do[A[x_] = 1/Product[1 - x^k*A[x]^(2*k), {k, 1, nmax}] + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x] (* Vaclav Kotesovec, Sep 26 2023 *)
(* Calculation of constants {d, c}: *) eq = FindRoot[{s*QPochhammer[r*s^2] == 1, 1/s + 2*r*s^2*Derivative[0, 1][QPochhammer][r*s^2, r*s^2] == (2*(Log[1 - r*s^2] + QPolyGamma[0, 1, r*s^2]))/(s* Log[r*s^2])}, {r, 1/8}, {s, 1}, WorkingPrecision -> 1000]; {N[1/r /. eq, 100], val = -s* Log[r*s^2]*(Sqrt[1 - r*s^2]/ Sqrt[4*Pi*(16*r*s^2*ArcTanh[1 - 2*r*s^2] + (1 - r*s^2)*(Log[r*s^2] - 2*Log[1 - r*s^2])*(3*Log[r*s^2] - 2*Log[1 - r*s^2]) - 8*Log[1 - r*s^2] + 8*(1 - r*s^2)*(-1 + 2*ArcTanh[1 - 2*r*s^2]) * QPolyGamma[0, 1, r*s^2] + 4*(1 - r*s^2)*QPolyGamma[0, 1, r*s^2]^2 - 4*(1 - r*s^2)*(QPolyGamma[1, 1, r*s^2] + r*s^2*Log[r*s^2]*(r*s^3*Log[r*s^2]* Derivative[0, 2][QPochhammer][r*s^2, r*s^2] - 2* Derivative[0, 0, 1][QPolyGamma][0, 1, r*s^2])))]) /. eq; N[Chop[val], -Floor[Log[10, Abs[Im[val]]]] - 3]} (* Vaclav Kotesovec, Sep 26 2023 *)
PROG
(PARI) a(n)=polcoeff((1/x*serreverse(x*eta(x+x*O(x^n))^2))^(1/2), n)
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 19 2009
STATUS
approved