扩展欧几里得算法

一、ax+by=gcd(a,b)的求解

问题：给定a,b两个非零整数，求满足ax+by=gcd(a,b)的整数解x、y，其中gcd(a,b)是a和b的最大公约数。

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

在求最大公约数的算法中，当到达递归边界b=0, a=gcd时，有(x1, y1)=(1,0)使ax1+by1=gcd(a,b)成立，假设bx2+(a%b)y2=gcd(b,a%b)，那么如果能找到(x1,y1)与(x2,y2)的关系，就可以用(x1, y1)=(1,0)为边界，在不断退出递归的过程中最终得到a=a0, b=b0时ax+by=gcd(a,b)的解（设a0, b0为a, b的初始值）。
由
$$\{\begin{array}{rl}& a{x}_1+b{y}_1=gcd(a,b)\\ & b{x}_2+(a\%b)y_2=gcd(b,a\%b)\\ & gcd(a,b)=gcd(b,a\%b)\end{array}$$
得：ax1+by1=bx2+(a%b)y2，又因为a%b=a-(a/b)*b（其中a/b取整数），带入整理可以得到：ax₁+by₁=ay₂+b[x₂-(a/b)y₂]
因此有：
$$\{\begin{array}{rl}{x}_1& =y_2\\ y_1& =x_2-(a/b)y_2\end{array}$$
因此，求一个特解(x0,y0)的算法如下：

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

显然，ax+by=gcd(a,b)可能有多个解，假设通解为
$$\{\begin{array}{rl}& x=x_0+s_1\\ & y=y_0+s_2\end{array}$$
由
$$\{\begin{array}{rl}& a{x}_0+b{y}_0=gcd(a,b)\\ & ax+by=gcd(a,b)\end{array}$$
得：as₁+bs₂=0，所以s₁/s₂=-b/a=-(b/gcd)/(a/gcd)，因为b/gcd与a/gcd互质，所以s₁与s₂最小分别为b/gcd、-a/gcd（正负反过来也可以），这就得到了ax+by=gcd(a,b)的通解：
$$\{\begin{array}{rl}& x=x_0+K\ast b/gcd\\ & y=y_0-K\ast a/gcd\end{array}$$
**PS:**由任意一组解得到最小非负整数解，以x为例，因为x%gcd可能为负，所以最小非负整数解为(x%gcd+gcd)%gcd.

二、ax+by=c的求解

ax+by=c有解的充要条件是：c%gcd=0（我只知道ax+by=c有解一定可以得到c%gcd=0，因为a和b都是gcd的倍数，所以c必须也是）。
在ax+by=gcd(a,b)两边同乘c/gcd得到a(cx/gcd)+b(cy/gcd)=c. 因此，将ax+by=gcd(a,b)的特解(x0, y0)乘上c/gcd就得到了ax+by=c的特解，显然，ax+by=gcd(a,b)的通解乘上c/gcd也是ax+by=c的解（Kb/gcd变为Kb/gcd²，Ka/gcd变为Ka/gcd²）；但是，按照上面求通解的步骤：
$$\{\begin{array}{rl}& x=\frac{c{x}_0}{gcd}+s_1\\ & y=\frac{c{y}_0}{gcd}+s_2\end{array}$$
$$\{\begin{array}{rl}& a\frac{c{x}_0}{gcd}+b\frac{c{y}_0}{gcd}=c\\ & ax+by=c\end{array}$$
依然得到s₁/s₂=-b/a=-(b/gcd)/(a/gcd)，所以通解为：
$$\{\begin{array}{rl}& x=\frac{c{x}_0}{gcd}+K\ast b/gcd\\ & y=\frac{c{y}_0}{gcd}-K\ast a/gcd\end{array}$$

三、ax≡c(modm)求解

如果(a-b)%m=0，那么a与b模m同余，写为a≡b(modm)。
因为ax≡c(modm)，所以(ax-c)%m=0，即ax-c=my，即ax-my=c，这就变成了ax+by=c求解，其中b=-m，只需要求x的解就可以了。

参考：《算法笔记》5.7