OFFSET
0,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..160
FORMULA
E.g.f. satisfies: A(x) = exp( Series_Reversion( Integral 1/A(x) dx ) ).
E.g.f.: exp(G(x)) where G(x) = exp(G(G(x))) is the e.g.f. of A214645.
EXAMPLE
E.g.f.: A(x) = 1 + x + 2*x^2/2! + 7*x^3/3! + 38*x^4/4! + 292*x^5/5! + 2975*x^6/6! +...
Related expansions.
Integral 1/A(x) dx = x - x^2/2! - x^4/4! - 6*x^5/5! - 52*x^6/6! - 591*x^7/7! - 8404*x^8/8! +...
The series reversion of Integral 1/A(x) dx equals log(A(x)) and begins:
log(A(x)) = x + x^2/2! + 3*x^3/3! + 16*x^4/4! + 126*x^5/5! + 1333*x^6/6! + 17895*x^7/7! + 293461*x^8/8! +...+ A214645(n)*x^n/n! +...
and equals the e.g.f. of A214645.
PROG
(PARI) {a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=exp(serreverse(intformal(1/A+x*O(x^n))))); n!*polcoeff(A, n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 07 2013
STATUS
approved