Python 2.7中计算超大分数遇循环重复问题求助

我来帮你搞定这个迭代计算分数的麻烦！先给你拆解问题根源，再给你直接能用的解决方案。

问题原因分析 #

你遇到的两个问题本质上都是精度丢失或者不必要的近似截断导致的：

• 用limit_denominator后陷入循环：每次迭代就截断分数到分母小于1e6的近似值，当迭代到某个阶段，近似后的分数会和之前的某个值重复，导致迭代进入死循环，再也无法推进。
• 不用limit_denominator时重复输出：你一开始用了浮点数1/2.0来初始化，而Python 2.7的float是双精度类型，只有15-17位有效数字。迭代15次后，浮点数已经无法存储新的有效数值，后续计算出来的结果都是同一个值，转换成分数自然就重复了。

解决方案：全程用精确分数运算 #

核心思路是彻底抛弃浮点数，全程用Fraction对象进行精确计算，只在最后需要近似的时候再调用limit_denominator。这样即使迭代1000次，也能保持完全精确的分数值，不会出现循环或重复。

修改后的代码示例 #

假设你的迭代逻辑是类似牛顿法逼近某个值（比如√2，从你给出的重复分数来看很像），修改后的代码如下：

from fractions import Fraction

i = 0
# 用精确分数初始化，代替浮点数1/2.0
x = Fraction(1, 2)
max_iterations = 1000

while i < max_iterations:
    # 这里替换成你的迭代公式，比如牛顿法求√2的公式是 x = (x + 2/x)/2
    # 所有运算都用Fraction对象，保证精确性
    x = (x + Fraction(2, 1) / x) / Fraction(2, 1)
    i += 1
    # 可以每隔100次输出一次进度，避免刷屏
    if i % 100 == 0:
        print("Iteration {}: Fraction = {}".format(i, x))

# 最后如果需要近似到分母小于1,000,000的分数，只调用一次limit_denominator
final_approx = x.limit_denominator(1000000)
print("\nFinal approximate fraction (denominator < 1,000,000): {}".format(final_approx))

关键注意事项 #

1. 绝对不要用浮点数初始化或参与迭代：浮点数的精度限制会直接导致迭代到一定次数后无法得到新值，必须全程用Fraction来表示所有数值。
2. 不要在迭代过程中调用limit_denominator：这个方法是用来做近似的，中途调用会破坏迭代的精确性，只有当你最终需要一个分母较小的近似分数时，再调用一次即可。
3. Python 2.7的Fraction完全支持大分数运算：即使迭代1000次，得到的分子分母会很大，但Fraction可以完美处理，性能上也完全没问题。

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