OFFSET
0,3
FORMULA
a(n) = n! * Sum_{k=0..floor(n/2)} (n-2*k)! * |Stirling1(n-k,n-2*k)|/(n-k)!.
a(n) ~ n! * (1 - r^2) / ((r^2 + 2*r - 1) * r^n), where r = 0.714556384743009681601449126434362887596497938663830826955917... (see A201750). - Vaclav Kotesovec, Sep 30 2025
MATHEMATICA
nmax = 20; CoefficientList[Series[1/(1 + Log[1 - x^2]/x), {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Sep 30 2025 *)
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(1/(1+log(1-x^2)/x)))
(PARI) a(n) = n!*sum(k=0, n\2, (n-2*k)!*abs(stirling(n-k, n-2*k, 1))/(n-k)!);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 29 2024
STATUS
approved
