OFFSET
1,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..500 (terms 1..210 from Paul D. Hanna)
FORMULA
a(n) ~ c * d^n / n^(3/2), where d = 5.023750195318838381709... and c = 0.049554511600191... - Vaclav Kotesovec, Sep 02 2017
EXAMPLE
G.f.: A(x) = x + x^2 + 2*x^3 + 6*x^4 + 20*x^5 + 73*x^6 + 281*x^7 +...
Series_Reversion(A(x)) = x - x*B(x) where
B(x) = x + x^3 + 3*x^5 + x^6 + 12*x^7 + 9*x^8 + 59*x^9 + 67*x^10 +...
Series_Reversion(B(x)) = x - x^2*C(x) where
C(x) = x + x^4 + 4*x^7 + x^8 + 22*x^10 + 12*x^11 + 5*x^12 + 141*x^13 +...
Series_Reversion(C(x)) = x - x^3*D(x) where
D(x) = x + x^5 + 5*x^9 + x^10 + 35*x^13 + 15*x^14 + 6*x^15 + x^16 +...
Series_Reversion(D(x)) = x - x^4*E(x) where
E(x) = x + x^6 + 6*x^11 + x^12 + 51*x^16 + 18*x^17 + 7*x^18 + x^19 +...
Series_Reversion(E(x)) = x - x^5*F(x) where
F(x) = x + x^7 + 7*x^13 + x^14 + 70*x^19 + 21*x^20 + 8*x^21 + x^22 +...
Series_Reversion(F(x)) = x - x^6*G(x) where
G(x) = x + x^8 + 8*x^15 + x^16 + 92*x^22 + 24*x^23 + 9*x^24 + x^25 +...
PROG
(PARI) {a(n)=local(G=x); m=sqrtint(2*n+1); for(k=0, m, G=serreverse(x-x^(m-k+1)*G+x*O(x^n))); polcoeff(G, n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 07 2011
STATUS
approved