算法学习（十四）最大公约数问题

题目描述：
求两个正整数的最大公约数，如果两个正整数都很大，有什么简单的算法？
例如：给定1100100210001,120200021，求其最大公约数。
分析：
求最大公约数问题最经典的就是“辗转相除法”，其实就是一个数学问题了，来自欧几里得的著作《几何原本》，原来也就是知道这么个方法，但是具体的原理还真没考虑过，今天就证明一下：
假设f(x,y)是x,y的最大公约数，这里有k = x/y,b = x%y，有x= k*y + b。如果一个数能够同时整除x,y,那么必能同时整除y和b,所以x,y的最大公约数和b,y的最大公约数相同。则f(x,y) = f(y,x%y)（x>=y>0），将原问题转化为求两个更小的数的问题。
所以解法一，就是辗转相除法。

int gcd(int x,int y)
{
    if(y != 0)
    {
        return gcd(y,x%y);
    }
    else
        return x;
}

我们用到了取模运算，对于大整数而言，取模运算时非常昂贵的开销，将成为整个算法的瓶颈。
我们能不能不用除法呢？
我们仿照欧几里得的证明方法，如果一个数能同时整除x和y,必能整除x-y,y，而反过来也是可以的，所以f(x,y) = f(x-y,y)。我们就不需要进行除法运算，转换成简单的减法运算。
解法二：

int  gcd(int x,int y)
{
    if(x < y)
        return gcd(y,x-y);
    if(y == 0)
        return x;
    else
        return gcd(x,x-y);
}

虽然避免了大整数除法，但是使用减法，在大整数面前，迭代次数也是够够的。
能不能将解法一和解法二结合呢？
分析公约数特点：
如果x= k* x1;y = k*y1，那么f(x,y) = k*f(x1,y1)。
如果x = p* x1,假设p是素数，且y%p != 0,那么f(x,y) = f(p*x1,y) = f(x1,y)。
那么素数p如何选择呢，2是一个特殊的素数，因为它在二进制中如鱼得水，x/2 = x>>1 。
取p = 2;
若x,y都是偶数，f(x,y) = f(x/2,y/2) = 2*f(x>> 1,y>>1)
若x为偶，y为奇数，f(x,y) =f(x/2,y) = f(x>>1,y)
若x为奇数，y为偶，f(x,y) = f(x,y/2) = f(x,y>>1)
若x,y均为奇数，f(x,y) = f(y,x-y)