OFFSET
0,3
FORMULA
G.f.: exp((1/3) * Sum_{k>=1} A000670(2*k) * x^k/k).
a(n) ~ sqrt(Pi) * 2^(2*n) * n^(2*n - 1/2) / (3 * exp(2*n) * log(2)^(2*n+1)). - Vaclav Kotesovec, May 29 2025
PROG
(PARI) a000670(n) = sum(k=0, n, k!*stirling(n, k, 2));
my(N=20, x='x+O('x^N)); Vec(exp(sum(k=1, N, a000670(2*k)*x^k/k)/3))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 28 2025
STATUS
approved
