OFFSET
0,2
COMMENTS
In general, for m > 1, Sum_{k=0..n} (m*k)! * Stirling1(n,k) ~ c * (m*n)!, where c = exp(-1/8) if m = 2, and c = 1 if m > 2. - Vaclav Kotesovec, Jan 24 2026
FORMULA
a(n) = Sum_{k=0..n} (3*k)! * Stirling1(n,k).
a(n) ~ sqrt(2*Pi) * 3^(3*n + 1/2) * n^(3*n + 1/2) / exp(3*n). - Vaclav Kotesovec, Jan 24 2026
MATHEMATICA
Table[Sum[(3*k)!*StirlingS1[n, k], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jan 24 2026 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(sum(k=0, N, (3*k)!*log(1+x)^k/k!)))
(PARI) a(n) = sum(k=0, n, (3*k)!*stirling(n, k, 1));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 21 2022
STATUS
approved
