最大公约数问题 (python)

最大公约数，在python的math库中有直接调用的函数gcd()。

解法 1

最简单的方法是使用辗转相除法。利用递归能够很轻松的解决这个问题。

def gcd(x,y):
    return x if y<=0 else gcd(y,x%y)
print(gcd(10,5))

解法 2

在解法1中用到了取模运算。但是对于大整数而言，取模运算是非常昂贵的开销。在解法2中我们用到类似辗转相除法的大整数减法，如果一个数能够同时除以x和y，则必能同时整除x-y和y。及F(x,y)=F(x-y,y)。

def gcd2(x, y):
    if x < y:
        return gcd2(y, x)
    if y == 0:
        return x
    else:
        return gcd2(x-y, y)
n,m = map(int,input().split())
print(gcd2(n,m))

解法 3

对于公约数有如下特点：

1. 对于y和x来说，如果y=k*y1,x=k*x1,那么有f(y,x) = k*(y1,x1)。
2. 如果x=p*x1,假如p是素数，并且y%p！=0，那么有f(x,y)=f(p*x1,y)=f(x1,y)

2是一个素数，对于二进制表示的大整数而已，可以很容易地将除以2和乘以2的运算转换为移位运算，从而避免大整数除法。

1. 取p=2
2. 若x，y均为偶数，f(x,y) = 2*f(x/2,y/2)=2*f(x>>1,y>>1)
3. 若x为奇数，y为偶数，f(x,y) = f(x,y/2)=f(x,y>>1)
4. 若x为偶数，y为奇数，f(x,y) = f(x/2,y)=f(x>>1,y)
5. 若x，y均为奇数，f(x,y) =f(y,x-y)

def gcd3(x,y):
    if x < y:
        return gcd3(y,x)
    if y == 0:
        return x
    else:
        if x%2 == 0:
            if y%2 == 0:
                return gcd3(x>>1,y>>1)<<1
            else:
                return gcd3(x>>1,y)
        else:
            if y%2 == 0:
                return  gcd3(x,y>>1)
            else:
                return  gcd3(y,x-y)
print(gcd3(10000,60))