最大公约数及求解x*a+y*b=1

#-*- coding:utf-8 -*-

def gcd(x,y):
    #保证x、y为整数
    if type(x) is not int or type(y) is not int :
        return 'x and y must be integer'
    #保证x、y为正数
    if x <= 0 or y <= 0:
        return 'x and y must be positive number'
    #将x、y按大到小排序,保证(x > y)
    x,y = sorted((x,y),reverse=True)
    A, B = x, y
    rt = []
    '''
    x1 = a1 * x2 + x3    x3 = -a1 * x2 + x1   x3 = p1 * x2 + q1 * x1
                                                p1 = -a1  q1 = 1
    x2 = a2 * x3 + x4    x4 = -a2 * x3 + x2   x4 = p2 * x2 + q2 * x1
                                                p2 = -a2*p1+1 q2 = -a2*q1
    x3 = a3 * x4 + x5    x5 = -a3 * x4 + x3   x5 = -a3*(p2 * x2 + q2 * x1) + (p1 * x2 + q1 * x1)
                                                x5 = (-a3*p2 + p1)*x2 + (-a3*q2 + q1)*x1
    ...
    x = temp * y + z
    x = y , y = z (z=x%y)
                                        
    '''
    #a1
    temp = A/B
    p1 = -temp 
    q1 = 1
    rt.append('%4d = %4d*%-4d + %4d'%(A,temp,B,A%B))
    #说明不互质,即有最大公约数,输出最大公约数
    if A%B == 0:
        rt.append('gcd(%d,%d)=%d'%(B,A,B))
        return rt
    #最大公约数为1,即互质,输出等式  1 = x*a+y*b
    if A%B == 1:
        rt.append(' 1 = %d*%d + %d*%d'%(p1,y,q1,x))
        return rt
    
    A, B = B, A%B
    #a2
    temp = A/B
    p2 = -temp*p1+1
    q2 = -temp*q1

    
    
    while True:
        rt.append('%4d = %4d*%-4d + %4d'%(A,A/B,B,A%B))
        #print '%4d = %4d*%4d + %4d*%4d'%(A%B,p2,y,q2,x)
        #print '*'*30
        #同上
        if A%B == 0:
            rt.append('gcd(%d,%d)=%d'%(y,x,B))
            return rt
        if A%B == 1:
            rt.append(' 1 = %d*%d + %d*%d'%(p2,y,q2,x))
            return rt
            
        temp1 ,temp2 = p2,q2
        
        A, B = B, A%B
        temp = A/B
        
        p2 = -temp*p2 + p1
        q2 = -temp*q2 + q1

        p1, q1 = temp1,temp2

if __name__ == '__main__':       
    rt = gcd(71,192)
    for v in rt:
        print v
运行结果
 192 =    2*71   +   50
  71 =    1*50   +   21
  50 =    2*21   +    8
  21 =    2*8    +    5
   8 =    1*5    +    3
   5 =    1*3    +    2
   3 =    1*2    +    1
 1 = -73*71 + 27*192


转载于:https://my.oschina.net/cheerfullchen/blog/282164

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值