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 外部完成。