editing
approved
editing
approved
a(n) ~ LambertW(3*exp(4))^n * n^(n-1) / (sqrt(3*(1 + LambertW(3*exp(4)))) * exp(n) * (-3 + LambertW(3*exp(4)))^(4*n + 1)). - Vaclav Kotesovec, Nov 07 2023
Table[1/(3*n+1)! * Sum[(3*n+k)! * Abs[StirlingS1[n, k]], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Nov 07 2023 *)
approved
editing
proposed
approved
editing
proposed
a(n) = (1/(3*n+1)!) * Sum_{k=0..n} (3*n+k)! * |Stirling1(n,k)|.
a(n) = (1/(3*n+1)!) * Sum{k=0..n} (3*n+k)! * |Stirling1(n,k)|.
1, 1, 9, 167, 4780, 186004, 9173780, 548563140, 38573633016, 3119384230176, 285237426927552, 29102185296785160, 3277703460197645232, 403931173342682581296, 54066960915411480743520, 7811249803193620134996864, 1211525560869437165319590400
(PARI) a(n) = sum(k=0, n, (3*n+k)!*abs(stirling(n, k, 1)))/(3*n+1)!;