OFFSET
1,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..100
FORMULA
E.g.f.: A(x) = Series_Reversion( exp(x) - exp(x^2) ).
a(n) = -sum(k=1..n-1, (k^n/n!+(k^(n/2)*(-1)^k*((-1)^n+1))/(2*(n/2)!)+sum(j=1..k-1, (-1)^(k-j)*binomial(k,j)*sum(m=1..n, (j^(2*m-n)*(k-j)^(n-m)*binomial(m,n-m))/m!)))*a(k)), n>1, a(1)=1. - Vladimir Kruchinin, Jun 25 2011
a(n) ~ n^(n-1) * sqrt((2*s-1)/(2-2*s+4*s^2)) / (exp(1+s)*(1-1/(2*s)))^n, where s = 0.6310764773894916166238... is the root of the equation exp(s) = 2*s*exp(s^2). - Vaclav Kotesovec, Jan 06 2014
MATHEMATICA
Rest[CoefficientList[InverseSeries[Series[E^x-E^(x^2), {x, 0, 20}], x], x] * Range[0, 20]!] (* Vaclav Kotesovec, Jan 06 2014 *)
PROG
(PARI) {a(n)=if(n<1, 0, n!*polcoeff(serreverse(exp(x+x*O(x^n))-exp(x^2+x*O(x^n))), n))}
(Maxima)
a(n):=if n=1 then 1 else -sum((k^n/n!+(k^(n/2)*(-1)^k*((-1)^n+1))/(2*(n/2)!)+sum((-1)^(k-j)*binomial(k, j)*sum((j^(2*m-n)*(k-j)^(n-m)*binomial(m, n-m))/m!, m, 1, n), j, 1, k-1))*a(k), k, 1, n-1); /* Vladimir Kruchinin, Jun 25 2011 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 27 2008
EXTENSIONS
a(19)-a(20) from Vincenzo Librandi, Feb 20 2018
STATUS
approved