最大公约数&最小公倍数(欧几里德算法相关)

//欧几里德算法即为辗转相除法，用于求两个整数的最大公约数
int main()   /*  辗转相除法求最小公倍数，最大公约数 */
{
    int m, n, a, b, c;
    scanf("%d%d", &a,&b);
    m=a;
    n=b;
    while(b!=0)
    { c=a%b; a=b;  b=c;}
    printf("%d\n", a);
    printf("%d\n", m*n/a);
    return 0;
}
int main ()  /* 相减法求最小公倍数，最大公约数 */
{
    int m, n, a, b;
    scanf ("%d%d", &a, &b);
    m=a;n=b;
    while(a!=b){
        if(a>b)
            a=a-b;
        else
            b=b-a;
    }
    printf("%d\n", a);
    printf("%d\n", m*n/a);
}

补充

由辗转相除法也可以推出，两数的最大公约数可以用两数的整数倍相加来表示，如21 = 5 × 105 + (−2) × 252。这个重要的结论叫做裴蜀定理。
欧几里得算法定理：gcd(a,b) = gcd(b,a mod b)

扩展欧几里得算法
已知整数a、b，扩展欧几里得算法可以在求得a、b的最大公约数的同时，能找到整数x、y（其中一个很可能是负数），使它们满足裴蜀等式。即gcd（a，b）=a*x+b*y 。
扩展欧几里德常用在求解模线性方程及方程组中。

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

Stein算法
这个方法也是计算两个数的最大公约数。和欧几里德算法 算法不同的是，Stein算法只有整数的移位和加减法。

int gcd(int a, int b)
{
    if(a == 0) return b;
    if(b == 0) return a;
    if(a % 2 == 0 && b % 2 == 0) return 2 * gcd(a >> 1, b >> 1);
    else if(a % 2 == 0)  return gcd(a >> 1, b);
    else if(b % 2 == 0) return gcd(a, b >> 1);
    else return gcd(abs(a - b), min(a, b));
}