Pi spigot算法Python实现偶现错误 进位相关digit偏差问题求助

我按照相关文档实现了Pi spigot算法的Python脚本，绝大多数数位计算正确，但偶尔会出现固定模式的错误：某一位digit比实际值小1，下一个digit输出10而非0，之后恢复正确，此类错误会重复出现。错误实例如下：

实现代码如下：

def denom(n): return 2*n+1
def nume(n): return n

count = 10000

array = [2]*count
carry = [0]*count
carried = [0]*count
remainder = [0]*count
quotient = [0]*count

pi = "<1000 digits of pi have been excluded>"

idx = 0

while True:
    array = list(map(lambda x:10*x, array))
    for i in range(count):
        i = count-1-i
        
        carried[i] = array[i] + carry[i]
        
        if i == 0:
            remainder[i] = carried[i] % 10
            break
        else:
            remainder[i] = carried[i] % denom(i)
        quotient[i] = carried[i] // denom(i)
        carry[i-1] = quotient[i] * nume(i)
    if not int(pi[idx]) == carried[i] // 10:
    # if True:
        print(carried[i] // 10, end="  ")
        print(pi[idx], end="  ")
        print(idx+1)
    if idx == 610:
        break
    array = remainder
    carry = [0]*count
    carried = [0]*count
    remainder = [0]*count
    quotient = [0]*count
    
    idx += 1

看起来“10”是代码尝试向前一位进位，但我原以为spigot算法的核心优势是可独立计算任意数位，为何此处需要计算后续数位来处理进位？

问题原因与解决思路 #

你遇到的是Spigot算法的伪进位陷阱，这是这类算法的常见问题，并非你对算法核心优势的理解有误：

1. "独立计算数位"的前提
   所谓"独立计算任意数位"，是指不需要从头计算所有前置数位就能定位到目标位，但算法本身依赖余数链的传递——它用整数近似生成数位，当后续余数累积到阈值时，当前位的初始计算值会存在偏差，需要回溯修正。

2. 错误的本质
   你看到的"某一位少1、下一位输出10"，本质是当前位的余数累积值本该触发进位，但代码没有处理回溯逻辑：当前位计算值少1，后续运算试图通过输出10来补进位，但这不符合数位规范，必须调整前一位的值并将当前位设为0。

3. 代码修复方向

   • 在生成数位后，增加检查步骤：如果后续余数链的累积值≥5（不同实现阈值略有差异），则将当前输出的数位加1；
   • 禁止输出超过9的数位，一旦出现这种情况，立即回溯将前一位加1，当前位设为0。

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