OFFSET
0,3
PROG
(Python)
from mpmath import mp, gammainc, exp, nint
def a(n: int) -> int:
if n == 0: return 1
mp.dps = 50
val = sum((-1)**(k + 1) * k**(n - k) * gammainc(k, -k) / exp(k) for k in range(1, n + 1))
return int(nint(val))
print([a(n) for n in range(21)])
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Mar 09 2026
STATUS
approved
