OFFSET
0,3
FORMULA
a(0) = 1; a(n) = 2 * n * Sum_{k=2..n} binomial(n-1,k-1) * a(n-k).
a(n) = n! * Sum_{k=0..floor(n/2)} 2^k * k! * Stirling2(n-k,k)/(n-k)!.
MATHEMATICA
terms = 21; CoefficientList[Series[1/(1 - 2 * x * (Exp[x] - 1)), {x, 0, terms}], x]*Range[0, terms]! (* Stefano Spezia, Dec 09 2025 *)
PROG
(PARI) a(n) = n!*sum(k=0, n\2, 2^k*k!*stirling(n-k, k, 2)/(n-k)!);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Dec 03 2023
EXTENSIONS
a(20)-a(21) from Stefano Spezia, Dec 09 2025
STATUS
approved
