PE-05 欧几里得算法

欧几里得算法：$$gcd(a,b)=gcd(b,a\%b)$$

适用场景：用于求最大公约数。

#include <stdio.h>
int gcd(int a, int b)
{
    return b ? gcd(b, a % b) : a;
}
int lcm(int a, int b)
{
    return a / gcd(a, b) * b;
}
int main()
{
    int a, b;
    while (~scanf("%d%d", &a, &b))
    {
        printf("%d\n", gcd(a, b));
        printf("%d\n", lcm(a, b));
    }
}

备注：通过最大公约数（gcd）可求最小公倍数（lcm）

证明

求证 $gcd(a,b)=gcd(b,a\%b)$; – gcd 代表最大公约数

1. 假设$a>b，t=gcd(a,b)；$则存在p1,p2使下列等式成立：
   $$\{\begin{array}{c}a=t\ast p1;\\ b=t\ast p2;\end{array}=>gcd(a,b)=t=t\ast gcd(p1,p2)$$
   由上式得到$gcd(p1,p2)=1$，即p1和p2互素（互为质数）
2. 假设$a\%b=r$，设$a=k\ast b+r$
3. 将1代入2得：$t\ast p1=k\ast t\ast p2+r，所以r=(p1-k\ast p2)\ast t$，整理得：
   $$\{\begin{array}{c}b=p2\ast t;\\ a\%b=r=(p1-k\ast p2)\ast t;\end{array}$$
4. 设$gcd(p2,p1-k\ast p2)=m$，则：
   $$\{\begin{array}{c}p2=m\ast n1;\\ p1-k\ast p2=m\ast n2;\end{array}$$
   化简得：
   $$\{\begin{array}{c}p2=m\ast n1;\\ p1=m\ast (n2+k\ast n1)\end{array}$$
5. 因为步骤1得到了$gcd(p1,p2)=1$，结合步骤4知，m=1;即$gcd(p2,p1-k\ast p2)=1$；代入3即可得：$gcd(b,a\%b)=t$。所以：$$gcd(a,b)=gcd(b,a\%b)$$