浅谈扩展欧几里得算法

P1 拓展欧几里得

已知二元不定方程，其中 $a,b$ 为常数

$$ax+by=c(a,b\in Z)$$

若令 $k=\gcd (a,b)$ ，则有 $k|(ax+by)$

显然，只有 $k|c$ 时，不定方程才有整数解。

只需要求出
$$ax+by=k$$
的一组特解，我们设其为 $x_0$ 与 $y_0$，则有：
$$a{x}_0+b{y}_0=k$$
$$a\frac{x_{0}c}{k}+b\frac{y_{0}c}{k}=k$$

那么 $ax+by=c$ 的一组解则为：$x=\frac{x_{0}c}{k},y=\frac{y_{0}c}{k}$

那么怎么求 $ax+by=k$ 的解呢，我们把 $k=\gcd (a,b)$ 换回：
$$ax+by=\gcd (a,b)$$
若 $a\le b$ ，由欧几里得定理得：

$$\gcd (a,b)=\gcd (b,a\%b)$$
其中 $\%$ 为取模运算，设 $x_1$，$y_1$ 满足：
$$x_{1}b+y_1(a\%b)=\gcd (b,a\%b)$$
则：
$$ax+by=x_{1}b+y_1(a\%b)$$
由取模的意义得：

$$ax+by=x_{1}b+y_1(a-b\ast \lfloor \frac{a}{b}\rfloor )$$

化简得：
$$ax+by=y_{1}a+b(x_1-\lfloor \frac{a}{b}\rfloor y_1)$$

则一定有一组解为
$$x=y_1,y=x_1-\lfloor \frac{a}{b}\rfloor y_1$$

则问题已经转化为求
$$bx+(a\%b)y=\gcd (b,a\%b)$$
的解

很明显，这是一个递归，那么跳出条件是什么呢？
当 $b=0$ 时，一定有一组解
$x=1,y=0$ 满足
$$ax=gcd(a,0)$$
注：
$$gcd(0,k)=k(k\in Z)$$

代码：

void exgcd(ll a,ll b,ll &x,ll &y)
{
	if(!b)
	{
		x=1,y=0;
      		return;
	}
	else
	{
		int tmp=x;
      		x=y;
            	y=tmp+(a/b)*y
	}
}