numpy数组传入函数计算结果与单独传入相同数值结果不一致问题

这一问题的核心是NumPy固定精度整数类型的溢出与Python原生整数的任意精度特性差异，具体逻辑如下：

• Python原生int类型支持任意精度运算，数值不会溢出，你单独调用u(9)得到的3138105961是准确的计算结果。
• 你初始化X = np.array([8, 9])时，NumPy默认会为小整数分配np.int32类型（32位有符号整数），该类型的最大值为2^31 - 1 = 2147483647。计算u(9)时最高项是9^10 = 3486784401，已经远超int32的最大值上限，发生有符号整数溢出，符号位翻转后就会得到负数的错误结果，和原生Python计算结果不一致。
• 数值为8时没有出现问题，是因为8^10 = 1073741824，小于int32的最大值，运算全程没有溢出，因此结果和原生Python计算一致。

你可以通过显式指定数组的数据类型避免溢出：

1. 使用更高精度的整数类型：

X = np.array([8, 9], dtype=np.int64)

np.int64的最大值约为9e18，远大于9^10的数值，运算不会溢出。
2. 直接使用Python对象类型调用原生int运算：

X = np.array([8, 9], dtype=object)

该方案会完全复用Python原生整数的运算逻辑，和直接调用u(9)的计算规则完全一致。

内容的提问来源于stack exchange，提问作者Syrocco