OFFSET
1,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..1000
FORMULA
a(n) ~ c * d^n / n^(3/2), where d = 4.4379395570069006103370017839662020486541... and c = 0.0779400283864079018838324912195942307... - Vaclav Kotesovec, Mar 16 2026
EXAMPLE
G.f.: A(x) = x + x^2 + 2*x^3 + 5*x^4 + 15*x^5 + 49*x^6 + 168*x^7 + 598*x^8 + 2190*x^9 + 8196*x^10 + ...
Series_Reversion(A(x)) = x - x*B(x) where
B(x) = x + x^4 + 4*x^7 + x^9 + 22*x^10 + 13*x^12 + 140*x^13 + 6*x^14 + 136*x^15 + 970*x^16 + ...
Series_Reversion(B(x)) = x - x^3*C(x) where
C(x) = x + x^6 + 6*x^11 + x^13 + 51*x^16 + 19*x^18 + 8*x^20 + 506*x^21 + x^22 + 300*x^23 + ...
Series_Reversion(C(x)) = x - x^5*D(x) where
D(x) = x + x^8 + 8*x^15 + x^17 + 92*x^22 + 25*x^24 + 10*x^26 + x^28 + 1240*x^29 + ...
Series_Reversion(D(x)) = x - x^7*E(x) where
E(x) = x + x^10 + 10*x^19 + x^21 + 145*x^28 + 31*x^30 + 12*x^32 + x^34 + 2470*x^37 + ...
Series_Reversion(E(x)) = x - x^9*F(x) where
F(x) = x + x^12 + 12*x^23 + x^25 + 210*x^34 + 37*x^36 + 14*x^38 + x^40 + 4324*x^45 + ...
Series_Reversion(F(x)) = x - x^11*G(x) where
G(x) = x + x^14 + 14*x^27 + x^29 + 287*x^40 + 43*x^42 + 16*x^44 + x^46 + 6930*x^53 + ...
PROG
(PARI) {a(n) = my(A=x, m=sqrtint(n+1)); for(k=0, m, A = serreverse(x - x^(2*(m-k)+1)*A + x*O(x^n))); polcoef(GF=A, n)}
{upto(n) = a(n); Vec(GF)}
upto(30)
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 11 2026
STATUS
approved