OFFSET
0,3
FORMULA
a(0) = 0; a(n) = n*a(n-1) + binomial(n+3,4).
E.g.f.: x * (1+3*x/2+x^2/2+x^3/24) * exp(x) / (1-x).
PROG
(PARI) my(N=30, x='x+O('x^N)); concat(0, Vec(serlaplace(x*sum(k=0, 3, binomial(3, k)*x^k/(k+1)!)*exp(x)/(1-x))))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Dec 31 2023
STATUS
approved