欧几里得算法及其拓展欧几里得算法

欧几里得算法（辗转相除法）的原理（点击这里）就可以解决。

代码实现思路 ：

对于任意的正整数a, b， 都有
$a=kb+r(k,r\in N)$
$r=a\%b$

假设 c 为a, b的最大公约数

$c=gcd(a,b)$
$\therefore c|a,c|b$
($c|a=a$可以整除以$c$)
$\because r=a-kb$
$\therefore c|r$
$\therefore c=gcd(b,r)$

结论：$gcd(a,b)=gcd(b,a\%b)$

最大公约数（欧几里得算法）

递归版本：

int gcd(int a,int b) {
    return b ? gcd(b, a % b) : a;
}

非递归（迭代）版本：

int gcd(int a, int b) {
    while (b) {
        a = a % b;
        swap(a, b);
    }
    return a;
}

顺便安利一波求最小公倍数的代码

最小公倍数

int lcm(int a, int b) {
    return a / gcd(a, b) * b;  //先除后乘可以防止溢出
}

拓展欧几里得算法

简单来讲，就是求解贝祖等式的解（注意求解出来的x，y不是最小的整数解）形如

$ax+by=gcd(a,b)$
$a,b\in N$

证明：

当 $b=0$ 时，等式为 $ax+0\ast y=gcd(a,0)$
$\therefore ax=a$
$\therefore x=1$
此时 $y$ 我们令它为 $0$

也就是 x = 1, y = 0

$a{x}_1+b{y}_1=gcd(a,b)$
$b{x}_2+(a\%b)y_2=gcd(a,b)$
根据欧拉几得算法又得出 ：$gcd(a,b)=gcd(b,a\%b)$
$\therefore a{x}_1+b{y}_1=b{x}_2+(a\%b)y_2$
$\because a\%b=a-a/b\ast b$
$\therefore a{x}_1+b{y}_1=b{x}_2+(a-a/b\ast b)y_2$
$\therefore a{x}_1+b{y}_1=b{x}_2+a{y}_2-(a/b\ast b)y_2$
$\therefore a{x}_1+b{y}_1=a{y}_2+b(x_2-a/b\ast y_2)$

$\therefore x_1=y_2$
$\therefore y_1=x_2-a/b\ast y_2$
接下来的步骤直接递归就好了
代码实现：

int exgcd(int a, int b, int &x, int &y) {
	if (!b) {
		x = 1;
		y = 0;
		return a;
	}
	int tx, ty;
	int res = exgcd(b, a % b, tx, ty);
	x = ty;
	y = tx - a / b * ty;
	return res;
}

简化版

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