import numpy n = 1600 m = 1500 def C(n,k): s = 0 B = numpy.array([[0]*(n+1)]*(n+1)) for i in range(n+1): B[i][0]=1; B[i][i]=1; for j in range(1,i): B[i][j]=long(B[i-1][j-1])+long(B[i-1][j]) return B[n][m] % 1000000 print C(n,m)
OverflowError: Python int too large to convert to C long