Python初学者如何处理OverflowError: (34, 'Result too large')错误？

问题原因分析 #

你用通项公式出现OverflowError是因为浮点数的表示范围有限，当n增大到1476时，((sqrt(5)+1)/2)**(n-1)的结果超出了Python浮点数能存储的最大值，导致溢出。下面给你两种简单易懂的解决方案：

方案一：迭代计算斐波那契数（入门友好） #

Python的整数类型支持无限精度，不会出现溢出问题。我们可以直接迭代生成斐波那契数，同时检查它的位数：

def find_fib_index_with_digits(target_digits):
    # 初始化斐波那契数列的前两项，索引从1开始：F(1)=1, F(2)=1
    a, b = 1, 1
    index = 2
    while len(str(b)) < target_digits:
        # 生成下一项：F(n) = F(n-1) + F(n-2)
        a, b = b, a + b
        index += 1
    return index

# 找第一个1000位的项的索引
print(find_fib_index_with_digits(1000))

代码说明 #

• 用len(str(b))直接判断数字的位数，简单直观
• 迭代过程只保存前两项，内存占用极低
• 运行后会输出结果：4782

方案二：利用对数公式计算（数学推导，效率更高） #

作为数学专业学生，你可以用斐波那契数列的渐近性质直接推导，不用生成实际的大数字：

斐波那契数列的通项公式近似为：

F(n) ≈ φⁿ / √5，其中φ=(1+√5)/2≈1.61803398875

一个数有1000位，等价于它 ≥ 10⁹⁹⁹（因为10⁹⁹⁹是第一个1000位的数）。结合近似公式可得不等式：
φⁿ / √5 ≥ 10⁹⁹⁹

两边取以10为底的对数：
n * log₁₀(φ) - log₁₀(√5) ≥ 999

整理后解n：
n ≥ (999 + log₁₀(√5)) / log₁₀(φ)

用Python计算这个值，取最小的整数即可：

import math

phi = (1 + math.sqrt(5)) / 2
log_phi = math.log10(phi)
log_sqrt5 = math.log10(math.sqrt(5))

n = (999 + log_sqrt5) / log_phi
# 向上取整得到最小的整数索引
print(math.ceil(n))

说明 #

• 因为当n足够大时，通项公式中的ψⁿ项（ψ=(1-√5)/2≈-0.618）绝对值极小，可以忽略，所以近似公式足够准确
• 计算结果同样是4782，而且几乎瞬间出结果，效率远高于迭代法

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