OFFSET
0,3
FORMULA
a(n) = (2*n-1) * a(n-1) - (n-1)*(n-3) * a(n-2) - binomial(n-1,2) * a(n-3) for n > 2.
a(n) ~ 2^(-3/4) * 3^(1/4) * exp(-5/4 + sqrt(6*n) - n) * n^(n - 1/4) * (1 + sqrt(3)/(2*sqrt(2*n))). - Vaclav Kotesovec, Mar 20 2023
MATHEMATICA
With[{nn=20}, CoefficientList[Series[Exp[(x+x^2/2)/(1-x)], {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Jun 08 2023 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(exp((x+x^2/2)/(1-x))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Mar 15 2023
STATUS
approved