【算法】扩展欧几里得

扩展欧几里得是用来求出类似于

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

的方程，满足 $x$ 和 $y$ 为正整数。本博客会讲解此方程解法以及此算法的运用方法。

---

方程解法

根据辗转相除法（欧几里得法）可以得出

$$\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$$

设 $x_0$ 和 $y_0$ 是此方程中的一组特解，那么

$$ax+by=b{x}_0+(a\mathrm{mod}b)y_0$$

$$=b{x}_0+(a-b\lfloor \frac{a}{b}\rfloor )y_0$$

$$=b{x}_0+a{y}_0-b\lfloor \frac{a}{b}\rfloor y_0$$

$$=a{y}_0+b(x_0-\lfloor \frac{a}{b}\rfloor y_0)=\gcd (a,b)$$

可得出答案

$$x=y_0,y=x_0-\lfloor \frac{a}{b}\rfloor y_0$$

其中，$x$ 和 $y$ 初始时无法得出，所以可任意赋初值。

因为任何正整数与 $0$ 的最大公因数都为自己，所以当 $b=0$ ，且需满足原方程，$x=1,y=0$。

代码

#include<bits/stdc++.h>
using namespace std;
int a,b,x,y;
void exgcd(int a,int b,int &x,int &y)
{
	if(!b)
	{
		x=1;y=0;
		return ;
	}
	exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-(a/b)*y;
}
int main()
{
	scanf("%d %d",&a,&b);
	exgcd(a,b,x,y);
	printf("%d %d",&x,&y);
	return 0;
}

二元一次不定方程

扩展欧几里得的一个主要作用就是来求二元一次不定方程。

二元一次方程是一个类似于以下形式的一个方程

$$ax+by=c$$

我们发现，这个方程的形式类似于扩展欧几里得的形式，又根据裴蜀定理，则可以得出

$$\gcd (a,b)\mid c$$

那么如果不满足当前性质，则不能有解

方程特解

根据上述性质，我们可以得出 $\frac{c}{\gcd (a,b)}$ 一定为整数，则如果将扩展欧几里得原式两边同时乘上一个，一定满足扩展欧几里得性质（$x$ 和 $y$ 为非负整数）

$$a\frac{c{x}_0}{\gcd (a,b)}+b\frac{c{y}_0}{\gcd (a,b)}=c$$

那么特解即为：

$$x_1=\frac{c{x}_0}{\gcd (a,b)},y_1=\frac{c{y}_0}{\gcd (a,b)}$$

方程通解

我们再考虑构造出一个等式

$$a{x}_1+b{y}_1+dab-dab=c$$

$$a(x_1+db)+b(y_1-da)=c$$

因为我们需要满足扩展欧几里得的性质，所以需要保证 $da$ 和 $db$ 均为非负整数。通解即为：

$$x=x_1+kdb,y=y_1-kda$$

其中，$k$ 为任意正整数

方程最值问题

首先，因为通解的和是一定的，所以 $y$ 一定随着 $x$ 的增大而减小。$k$ 越大，$x$ 越大，$y$ 越小，反之亦然。

所以，我们可以先求出 $x$ 和 $y$ 的最小值再带入求出它们的最大值。

如果我们要使 $x$ 最小，那么我们就要将 $x$ 减去 $k$ 个 $db$。但是我们并不知道 $k$，所以便要一直向下减一直到不能再减了（再减就不为正整数了），那我们就考虑有没有一种方法可以直接求出而不用一直减。

考虑与 $db$ 取模，取模后的答案，也就为 $x$ 的最小非负整数解。

那我们就要使 $db$ 尽可能的小，也就是使 $d$ 尽可能的小（$b$ 已经给你了），且还满足 $db$ 为非负整数。显然，$d=\frac{1}{gcd(a,b)}$

那么 $x$ 的最小非负整数解即为

$$x_2=(\frac{c{x}_0}{\gcd (a,b)}\mathrm{mod}\frac{b}{gcd(a,b)}+\frac{b}{gcd(a,b)})\mathrm{mod}\frac{b}{gcd(a,b)}$$

PS：注意这里一定要加一下再模一下，不然会为负数（我也不知道为什么会为负数）

将其设为 $x_2$ 随后带入求 $y$ 的最大非负整数解 $y_3$，即：

$$y_3=\frac{c-a{x}_2}{b}$$

反之亦然

代码

#include<bits/stdc++.h>
using namespace std;
int a,b,x,y,c;
int exgcd(int a,int b,int &x,int &y)
{
	if(!b)
	{
		x=1;y=0;
		return a;
	}
	int gcd=exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-(a/b)*y;
	return gcd;
}
int main()
{
	scanf("%d %d %d",&a,&b,&c);
	int gcd=exgcd(a,b,x,y);
	if(c%gcd)
		puts("impossible!");
	else
	{
		int minx=(((c*x)/gcd)%(b/gcd)+(b/gcd))%(b/gcd);
		int miny=(((c*y)/gcd)%(a/gcd)+(a/gcd))%(a/gcd);
		int maxx=(c-b*miny)/a;
		int maxy=(c-a*minx)/b;
		printf("%d %d %d %d",minx,miny,maxx,maxy);
	}
	return 0;
}