欧几里德算法求解方程但有一个bug

时间:2018-03-24 20:19:53

标签: python discrete-mathematics

def kek(k,p):
    if k > p:
        return pew(k,p)
    elif k < p:
        return pew(p,k)
    else:
        return 1
def pew(j,k):
    if k == 0:
        return j,j,k
    q = j // k
    r = j % k
    R, h , o = pew(k,r)
    return R,o,h-q*o    
hmm = kek(231,1920)

'''
Math logic:
d = s x 231 + t x 1920
1920 = 8*231 + 72
231 = 3*72 + 15
72 = 4*15 + 12
15 = 1*12 + 3
12 = 4*3
d = gcd(231, 1920) = 3
3 = 15 – 12 
= 15 – (72 – 4*15) 
= 15 – 72 + 4*15 
= 5*15 – 72 
= 5*(231-3*72) – 72
= 5*231 – 15*72 – 72 
= 5*231 – 16*72 
= 5*231 – 16*(1920-8*231)
= 5*231 – 16*1920 + 128*231
= 133*231 – 16*1920 
= 133*231 + (-16)*1920
d = 3, s = 133 , t = -16
'''

然而,我得到-48为s,399为t,正确结果为3。哪里 是我的递归方法的问题?上面的数学逻辑是你如何手动计算它,我使用D = b R1 +(aQ b) R2和D = a R2 + b * R3来实现我的呼叫,但它给了我133 * 3和-16 * 3的结果,这是我不想得到的。

1 个答案:

答案 0 :(得分:2)

您的代码很难遵循,但我想出了一个不同的解决方案,它可以做同样的事情。

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

print(egcd(231, 1920)) # (3, 133, -16)