python - 在 Python 中处理大量数字时遇到问题

Question

遇到以下问题：

from pylab import *
x=arange(0,365,1)
y = []
for j in x:
    y.append(1-((math.factorial(365)/math.factorial(365-j))/(365**j)))

plot(x,y)
show()

有什么想法吗？我正在运行 python 2.7

Answer 1

将其更改为：

(math.factorial(365)/math.factorial(365-j))/long(365**j)

问题是numpy在进行除法之前尝试将 long 转换为 numpy 类型——但是 numpy 不能这样做，因为它会导致溢出，所以 numpy 会引发错误。如果您转换为long第一个，则除法全部在外部完成numpy。