OFFSET
0,3
FORMULA
a(n) = n! * Sum_{k=0..n} (-3)^(n-k) * binomial(k+3,3)/(n-k)!.
a(0) = a(1) = 1; a(n) = n*a(n-1) + 3*(n-1)*a(n-2).
a(n) = A137775(n+2)/(3*(n+1)).
a(n) ~ sqrt(2*Pi) * n^(n + 7/2) / (6*exp(n+3)). - Vaclav Kotesovec, Apr 25 2025
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(exp(-3*x)/(1-x)^4))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Apr 24 2025
STATUS
approved
