login
A379996
Expansion of e.g.f. exp(-x)/(exp(-x) - x)^4.
3
1, 7, 61, 639, 7825, 109683, 1731645, 30403495, 587595649, 12395233539, 283385424829, 6979650164391, 184235963026833, 5188528445210035, 155284012799863453, 4921569063327156807, 164672737994453759617, 5800532536265417597571, 214559605389429001486557
OFFSET
0,2
FORMULA
a(n) = n! * Sum_{k=0..n} (k+3)^(n-k) * binomial(k+3,3)/(n-k)!.
MATHEMATICA
With[{nn=20}, CoefficientList[Series[Exp[-x]/(Exp[-x]-x)^4, {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, May 19 2025 *)
PROG
(PARI) a(n) = n!*sum(k=0, n, (k+3)^(n-k)*binomial(k+3, 3)/(n-k)!);
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Jan 07 2025
STATUS
approved