#define MAX 21

#macro ncp(n, p)
   (fact(n)/(fact(p))/(fact(n-p)))
#endmacro

dim as ulongint fact(0 to MAX), bell(0 to MAX)
dim as uinteger n=0, k

fact(0) = 1
for k=1 to MAX
    fact(k) = k*fact(k-1)
next k

bell(n) = 1
print n, bell(n)
for n=0 to MAX-1
    for k=0 to n
        bell(n+1) += ncp(n, k)*bell(k)
    next k
    print n+1, bell(n+1)
next n