欧几里得及扩展欧几里得

欧几里得算法又名辗转相除法，用于计算两个正整数a, b的最大公约数
公式 gcd(a, b) = gcd(b, a%b)（gcd()表示求最大公约数）

代码

int Gcd(int a, int b)
{
    while(b != 0)
    {
        int r = a % b;
        a = b;
        b = r;
    }
    return a;
}

顺便可以求出a, b的最小公倍数
lcm(a, b) = a*b/gcd(a, b);

所以有a * b = gcd(a, b) * lcm(a, b) , 即lcm(a, b) = a * b/gcd(a, b);
但是先算乘法可能会溢出， 所以写成
lcm(a, b) = a/gcd(a, b)*b;

欧几里得算法扩展：

int exgcd(int a, int b, int &x, int &y)
{
    if(b==0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int r = exgcd(b, a, %b, x, y);
    int t = y;
    y = t - a/b * y;
    return r;
}

求得一组特解（x0， y0）后，可以得到通解（x+kb’, y-ka’）,其中a’ = a/gcd(a, b), b’ = b/gcd(a, b), k取任意整数。

详解见扩展欧几里德算法详解