OFFSET
0,3
MATHEMATICA
a[1] = 1; a[n_] := Sum[(n - k)! * (n - 1)^k, {k, 0, n}]; Array[a, 20, 0] (* Amiram Eldar, Dec 25 2021 *)
PROG
(PARI) a(n) = sum(k=0, n, (n-k)!*(n-1)^k); \\ Michel Marcus, Dec 25 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Dec 25 2021
STATUS
approved