OFFSET
0,2
FORMULA
G.f.: exp(Sum_{k>=1} A000670(3*k) * x^k/k).
a(n) ~ sqrt(Pi) * 3^(3*n + 1/2) * n^(3*n - 1/2) / (sqrt(2) * exp(3*n) * log(2)^(3*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(3*k)*x^k/k)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 28 2025
STATUS
approved
