扩展欧几里得算法（exgcd）

前置知识

数论基本概念，欧几里得算法

裴蜀定理

在学习exgcd之前，我们需要先学习一下裴蜀定理。

裴蜀定理有两条：

1. 对于任意的整数 $a,b$，存在一组整数 $x,y$，使得 $ax+by=gcd(a,b)$
2. 若 $a,b$ 是整数，那么对于任意的整数 $x,y$，都有 $gcd(a,b)|ax+by$

第一条比第二条难证，我们先来证明第一条，下面是证明过程，我尽量讲的通俗易懂：

我们考虑辗转相除的最后一步，此时 $b=0,a=gcd(a,b)$，所以 $ax+by=gcd(a,b)$ 一定有解（ 解为 $x=1,y=0$ ）。

而辗转相除这一步的 $a$ 为上一步的 $b$，这一步的 $b$ 为上一步的 $a\mathrm{mod}b$，所以我们代进去，得： $$bx+(a\mathrm{mod}b)y=gcd(b,a\mathrm{mod}b)$$

因为最后一步的 $ax+by=gcd(a,b)$ 有解，所以 $bx+(a\mathrm{mod}b)y=gcd(b,a\mathrm{mod}b)$ 必定有解。我们设它的解为 $\{\begin{array}{l}x=x^{\prime }\\ y=y^{\prime }\end{array}$，则有：
$$b{x}^{\prime }+(a\mathrm{mod}b)y^{\prime }=gcd(b,a\mathrm{mod}b)$$

而我们知道：
$$gcd(b,a\mathrm{mod}b)=gcd(a,b)$$
$$a\mathrm{mod}b=a-\lfloor \frac{a}{b}\rfloor b$$

所以可得：
$$b{x}^{\prime }+(a-\lfloor \frac{a}{b}\rfloor b)y^{\prime }=gcd(a,b)$$

转化一下：
$$a{y}^{\prime }+b{x}^{\prime }-y^{\prime }\lfloor \frac{a}{b}\rfloor b=gcd(a,b)$$
$$a{y}^{\prime }+b(x^{\prime }-y^{\prime }\lfloor \frac{a}{b}\rfloor )=gcd(a,b)$$

我们发现，式子的左边成为了 $ax+by$ 的形式，所以我们得到，方程 $ax+by=gcd(a,b)$ 的解为：$$\{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-y^{\prime }\lfloor \frac{a}{b}\rfloor \end{array}$$

这样，对于任何的 $ax+by=gcd(a,b)$，我们都可以通过这种方法证得有解（数学归纳法）。

下面我们来证明第二条：

我们设 $d=gcd(a,b)$，则有
$$d|a,d|b$$

因为 $x,y$ 均为整数，所以
$$d|ax,d|by$$

所以得
$$d|ax+by$$

即
$$gcd(a,b)|ax+by$$

得证。

这样，我们就证明了裴蜀定理。我们来看一道例题：洛谷 P4549 【模板】裴蜀定理。

这个题目主要运用了裴蜀定理的第二条：若 $a,b$ 是整数，那么对于任意的整数 $x,y$，都有 $gcd(a,b)|ax+by$。那么，我们得到，$gcd(a,b)$ 是 $ax+by$ 的最小正整数值。

推广到多个形如 $A_i{X}_i$ 的单项式上，$\sum _{i=1}^n{A}_i{X}_i$ 的最小正整数值为 g c d { A i } gcd\{A_i\} gcd{Ai​}，也就是题目中的 $S$。

注意输入中有负数，取个绝对值就行了。

$Code$

#include<iostream>
#include<cstdio>
using namespace std;
int n;
int a[30];
int ans;
int gcd(int x,int y)
{
	if(y==0)
		return x;
	return gcd(y,x%y);
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
	{
		if(a[i]<0)
			a[i]=-a[i];
		ans=gcd(ans,a[i]);
	}
	printf("%d",ans);
	return 0;
}

扩展欧几里得算法

扩展欧几里得算法（exgcd），正如名字一样，建立在欧几里得算法的基础上。事实上，它还依赖于上述的裴蜀定理。

exgcd是干什么用的呢？我们看到裴蜀定理的第一条：

对于任意的整数 $a,b$，存在一组整数 $x,y$，使得 $ax+by=gcd(a,b)$

exgcd 就是求 $ax+by=gcd(a,b)$ 的解。

仔细阅读上述裴蜀定理第一条的证明过程，$ax+by=gcd(a,b)$ 的解是什么？

$\{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-y^{\prime }\lfloor \frac{a}{b}\rfloor \end{array}$

那么 $x^{\prime }$ 和 $y^{\prime }$ 又是怎么求出来的呢？它们是方程 $bx+(a\mathrm{mod}b)y=gcd(b,a\mathrm{mod}b)$ 的解。

可以发现，这是一个递归的过程，不断的递归求出 $x^{\prime }$ 和 $y^{\prime }$，从而得到 $x$ 和 $y$。

那么递归的边界又是什么呢？当 $b=0$ 时，解为 $\{\begin{array}{l}x=1\\ y=0\end{array}$ 。

int x,y;
void exgcd(int a,int b)
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	exgcd(b,a%b);
	int t=x;//暂存x
	x=y;
	y=t-a/b*y;//注意不能写成y*a/b，必须先计算a/b
}

来道例题：洛谷 P1082 [NOIP2012 提高组] 同余方程

我们转换一个形式：
$$ax\mathrm{mod}b=1$$

因为我们知道，$a\mathrm{mod}b=a-\lfloor \frac{a}{b}\rfloor \times b$，所以得：
$$ax-\lfloor \frac{ax}{b}\rfloor \times b=1$$

稍微改动，它成了这个形式：
$$ax+by=1$$

因为题目告诉我们 $ax\mathrm{mod}b=1$ 有解，所以 $ax+by=1$ 有解。

裴蜀定理第二条告诉我们，$gcd(a,b)|ax+by$，而 $ax+by=1$，则 $gcd(a,b)|1$，那么 $gcd(a,b)$ 只能等于 $1$。

把 $1$ 代换一下，方程又成了：
$$ax+by=gcd(a,b)$$

可以使用exgcd求解。

题目要求的是最小正整数解，而exgcd求出的可能更大，也可能是负数。那我们怎么把它转化为最小正整数解呢？

我们可以证明，所有的解之间的差均为 $b$，证明过程如下：

我们知道$ax+by=gcd(a,b)$，即 $\frac{gcd(a,b)-ax}{b}=y$，因为 $y$ 是整数，所以 $b|gcd(a,b)-ax$。

现在我们设 $x$ 增加了 $\Delta x$，则原式变为 $gcd(a,b)-ax-a\Delta x$。

而因为 $gcd(a,b)=1$ 即 $a,b$ 互质，则 $b|\Delta x$，得证。

那么，我们就可以通过加减 $b$ 的方式找到答案。具体做法么，直接模 $b$ 就行了……

$Code$

#include<iostream>
#include<cstdio>
using namespace std;
long long a,b;
long long x,y;
void exgcd(long long a,long long b)
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	exgcd(b,a%b);
	long long t=x;
	x=y;
	y=t-a/b*y;
}
int main()
{
	scanf("%lld%lld",&a,&b);
	exgcd(a,b);
	printf("%lld",(x%b+b)%b);
	return 0;
}