最大公约数与最小公倍数

1. 最大公约数

思路：最大公约数，是一个很典型的递归算法，每次递归使得原来求两个大数之间的公约数转变成求两个稍微小点的数之间的公约数。

解法一：辗转相除法（也叫欧几里得算法）：如果一个数能够同时整除x和y，则必能同时整除y和x%y，反之亦然，即f(x, y)=f(y, x%y)

(1) 对于已知两数x和y，x除以y记余数为r；

(2)若r==0，则y即为所求最大公约数，否则执行第(3)步；

(3)令x=y，y=r，重复(1)；

int gcd(int x, int y)
{
    return y == 0 ? x : gcd(y, x%y);
}

解法二：扩展欧几里德算法：如果一个数能够同时整除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);
    if(0 == y)
        return x;
    return gcd(x-y, y);
}

解法三：（摘自编程之美）

解法一的瓶颈是取模运算，对于大整数而言开销非常大，解法二问题在于减法的迭代次数太大。

分析如下：

对于x和y来说，令y=k*y1，x=k*x1，则有f(x, y)=k*f(x1, y1)

另外，如果x=p*x1，假设p是素数，且y%p!=0（即y不能被p整除），那么f(x, y)=f(p*x1, y)=f(x1, y)

根据以上思路，取p=2，可以将大整数的除法运算转化为移位运算，从而提高效率。

取p=2

若x，y均为偶数，则f(x, y)=2*f(x/2, y/2)=2*f(x>>1, y>>1)

若x为偶数，y为奇数，则f(x, y)=f(2*x/2, y)=f(x>>1, y)

若x为奇数，y为偶数，则f(x, y)=f(x, 2*y/2)=f(x, y>>1)

若x，y均为奇数，则f(x, y)=f(x, x-y)

bool isEven(int x)
{
    return !(x & 1);
}

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

2. 最小公倍数

对于x和y，假设其最大公约数为k，则最小公倍数为x*y/k

int lcm(int x, int y)
{
    return x*y / gcd(x,y);
}