OFFSET
0,3
FORMULA
a(n) = n! * Sum_{k=0..n} (3*n+1)^(k-1) * binomial(k,n-k)/k!.
PROG
(PARI) a(n) = n!*sum(k=0, n, (3*n+1)^(k-1)*binomial(k, n-k)/k!);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 17 2023
STATUS
approved