OFFSET
0,3
FORMULA
a(0) = 1; a(n) = Sum_{k=1..n} (-1)^(n-k) * Stirling1(n,k) * Bell(k) * 2^(k-1) for n > 0.
a(n) = A136658(n) / 2 for n > 0.
MATHEMATICA
nmax = 21; CoefficientList[Series[(1/2) (Exp[1/(1 - x)^2 - 1] + 1), {x, 0, nmax}], x] Range[0, nmax]!
Join[{1}, Table[Sum[(-1)^(n - k) StirlingS1[n, k] BellB[k] 2^(k - 1), {k, 1, n}], {n, 1, 21}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 03 2025
STATUS
approved
