最大公因数进阶

辗转相除法和相减法：https://blog.csdn.net/xinwenhuayu/article/details/80367966

当我们求两个数的最大公约数，可以用辗转相除法和相减法来求，但是辗转相除法的取模运算的代价很大，而相减法遇到类似（100000，1）这样的数则很尴尬，从分析公约数的特点入手，对于 gcd(a，b)

当 a = k * a1，b = k * b1 时，gcd(a，b) = k * gcd(a1，b1)

当 a = p * a1，p 是素数时，gcd(a，b) = gcd(a1，b)

通过以上两点，对算法进行改造，取 p = 2，方便左移右移

当 x，y 均为偶数，gcd(a，b) = (gcd(a>>1，b>>1) << 1)

当 x为偶数，y 为奇数，gcd(a，b) = gcd(a>>1，b)

当 x为奇数，y 为偶数，gcd(a，b) = gcd(a，b>>1)

当 x，y 均为奇数，gcd(a，b) = gcd(b，a-b)

代码如下： 

#include <iostream>
using namespace std;

template <class T>
T gcd(T a, T b)
{
    if (a < b)
        return gcd(b, a);
    if (b == 0)
        return a;
    else 
    {
        if ((a & 0x01) == 0)
        {
            if ((b & 0x01) == 0)
                return (gcd(a>>1, b>>1) << 1);
            else 
                return gcd(a>>1, b);
        }
        else 
        {
            if ((b & 0x01) == 0)
                return gcd(a, b>>1);
            else 
                return gcd(b, a-b);
        }
    }
}

int main()
{
    int a, b;
    while (cin >> a >> b)
    {
        cout << gcd(a, b) << endl;
    }
    return 0;
}