Python JIT素数判断代码处理超22位大数字时输出错误原因？

你的代码错误核心在于将输入的大整数转换为float类型。Python中的float是双精度浮点数，仅能精确表示53位以内的整数（约16位十进制数）。当输入的数字超过这个范围时，float会丢失精度，存储的是原数的近似值而非精确值。比如你输入的素数151978145606541879151，转成float后已经不是原来的精确整数，导致后续素数判断逻辑基于错误数值执行，从而得出错误结论。

直接将输入转换为Python的int类型即可——Python原生支持任意精度的整数，不会丢失精度。修改main函数中的调用代码：

def main():
    print(isprime(int(sys.argv[1])))

#!/home/fury/anaconda3/bin/python
from numba import jit
import sys
@jit(nopython=True)
def isprime(n):
    if n==2:return 1
    if n%2==0:return 0
    for i in range(3,int(n**0.5)+1,2):
        if n%i==0:return 0
    return 1
def main():
    print(isprime(int(sys.argv[1])))
if __name__=="__main__":main()

此时输入151978145606541879151，代码会正确返回1（表示是素数）。

• numba的nopython=True模式对大整数运算有支持，但Python原生int的精度保障才是解决该问题的核心。
• 对于超大素数的判断，试除法效率极低，若需处理更大的数，建议使用米勒-拉宾素性测试等更高效的算法。

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