为什么numpy.prod（）对于一长串自然数会错误地返回负数或0？

Question

我只是在研究Euler项目问题12，所以我需要对数字进行测试，这些数字是500多个独特因子的倍数。

我认为数组[1、2、3 … 500]是一个很好的起点，因为该数组的乘积是该数字可能的最低值。但是，numpy.prod（）对于此数组返回 零
。我确定我缺少明显的东西，但是到底是什么呢？

>>> import numpy as np
>>> array = []
>>> for i in range(1,100):
...   array.append(i)
... 
>>> np.prod(array)
0
>>> array.append(501)
>>> np.prod(array)
0
>>> array.append(5320934)
>>> np.prod(array)
0

Answer 1

请注意，Python使用“无限”整数，但是在numpy中，所有类型都被键入，因此这里是“
C”风格（可能是64位）整数。您可能正在经历溢出。

如果查看的文档numpy.prod，则可以看到dtype参数：

返回数组的类型，以及与元素相乘的累加器的类型。

您可以执行以下操作：

1. 回到Python，并使用其“无限整数”相乘（有关方法，请参阅此问题）。

2. 考虑您是否真的需要找到如此庞大的产品。通常，当您使用非常小或非常大的乘积时，会切换到对数和。正如@WarrenWeckesser指出的那样，这显然是不精确的（这不像最后采用指数会为您提供确切的解决方案），而是用来衡量一种产品是否比另一种产品增长更快。