实现扩展欧几里得算法

Question

为什么扩展Euclid算法的以下实现失败？

def extended_euclid(a,b):
    if b == 0:
        return {a, 1, 0}

    d1,x1,y1 = extended_euclid(b, a % b)
    d = d1
    x = y1
    y = x1 - math.floor(a/b) * y1
    return {d, x, y} 

Answer 1

def extended_euclid(a,b):
    if b == 0:
        return a, 1, 0

    d1,x1,y1 = extended_euclid(b, a % b)
    d = d1
    x = y1
    y = x1 - math.floor(a/b) * y1
    return d, x, y

从您的退货中删除{} 。 看一下d1,x1,y1 = extended_euclid(b, a % b) ，如果保留{}作为return值，则没有足够的值可用于解包。

Answer 2

这是在https://en.wikibooks.org/wiki/Algorithm_Implementation/Mathematics/Extended_Euclidean_algorithm中找到的实现，与您的实现相似。

def egcd(a, b):
    if a == 0:
        return (b, 0, 1)
    else:
        g, x, y = egcd(b % a, a)
        return (g, y - (b // a) * x, x)