OFFSET
0,2
FORMULA
a(n) ~ exp(-gamma) * n * n!, where gamma is the Euler-Mascheroni constant A001620.
MATHEMATICA
nmax = 25; CoefficientList[Series[Product[(1+x^(2*k-1)/(2*k-1))/(1-x^(2*k-1)/(2*k-1)), {k, 1, Floor[nmax/2]+1}], {x, 0, nmax}], x] * Range[0, nmax]!
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Jul 27 2019
STATUS
approved