23 lines
356 B
Plaintext
23 lines
356 B
Plaintext
#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
|