(PARI) for(n=0, 50, print1(ceil(sum(i=0, 1000, i^(2*n)/(i)!)/exp(1)), ", "))
(SageMath) [bell_number(2*n) for n in range(0, 17)] #
Zerinvary Lajos, May 14 2009
(Python)
from itertools import accumulate, islice
def
A020557_gen(): # generator of terms
yield 1
blist, b = (1, ), 1
while True:
for _ in range(2):
blist = list(accumulate(blist, initial=(b:=blist[-1])))
yield b