OFFSET
1,5
COMMENTS
Reversion of g.f. (with constant term omitted) for A003823.
LINKS
N. J. A. Sloane, Transforms
Eric Weisstein's World of Mathematics, Series Reversion
FORMULA
G.f. A(x) satisfies: A(x)/(1 + A(x)^2/(1 + A(x)^3/(1 + A(x)^4/(1 + A(x)^5/(1 + ...))))) = x.
a(n) ~ (-1)^(n+1) * c * d^n / n^(3/2), where d = 2.3794295463748306617... and c = 0.1900533719371157... - Vaclav Kotesovec, May 07 2024
MATHEMATICA
Rest[CoefficientList[InverseSeries[Series[ContinuedFractionK[x^i, 1, {i, 1, 36}], {x, 0, 36}], x], x]]
Rest[CoefficientList[InverseSeries[Series[-1 + QPochhammer[x^2, x^5] QPochhammer[x^3, x^5]/(QPochhammer[x, x^5] QPochhammer[x^4, x^5]), {x, 0, 36}], x], x]]
(* Calculation of constant d: *) -1/r /. FindRoot[{1 + r == QPochhammer[s^2, s^5]*QPochhammer[s^3, s^5] / (QPochhammer[s, s^5]*QPochhammer[s^4, s^5]), 5*s^4*QPochhammer[s^3, s^5] * Derivative[0, 1][QPochhammer][s^2, s^5] + (1/s)* QPochhammer[s^2, s^5]*((1/Log[s^5])*QPochhammer[s^3, s^5] * (QPolyGamma[0, Log[s]/Log[s^5], s^5] - 2*QPolyGamma[0, Log[s^2]/Log[s^5], s^5] - 3*QPolyGamma[0, Log[s^3]/Log[s^5], s^5] + 4*QPolyGamma[0, Log[s^4]/Log[s^5], s^5]) - (5*s^5*QPochhammer[s^3, s^5] * Derivative[0, 1][QPochhammer][s, s^5])/ QPochhammer[s, s^5] + 5*s^5*Derivative[0, 1][QPochhammer][s^3, s^5] - (5*s^5*QPochhammer[s^3, s^5] * Derivative[0, 1][QPochhammer][s^4, s^5])/ QPochhammer[s^4, s^5]) == 0}, {r, -2/5}, {s, -2/ 3}, WorkingPrecision -> 60] (* Vaclav Kotesovec, May 07 2024 *)
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Aug 23 2017
STATUS
approved