exgcd

题目传送门

对于不等式：$ax+by=gcd(a,b)$，给定$a,b$，求出符合条件的一组$x,y$。

推导： (假定a>=b)

1. $a{x}_1+b{y}_1=gcd(a,b)$
2. $b{x}_2+a\%b{y}_2=gcd(b,a\%b)$
3. $gcd(a,b)=gcd(b,a\%b)$
4. $a{x}_1+b{y}_1=b{x}_2+a\%b{y}_2$
5. $a\%b=a-\lfloor \frac{a}{b}\rfloor \times b$
6. $a{x}_1+b{y}_1=b{x}_2+(a-\lfloor \frac{a}{b}\rfloor \times b)y_2$
7. $a{x}_1+b{y}_1=a{y}_2+b{x}_2-\lfloor \frac{a}{b}\rfloor \times b{y}_2$
8. $a{x}_1+b{y}_1=a{y}_2+b(x_2-\lfloor \frac{a}{b}\rfloor \times y_2)$
9. $x_1=y_2,y_1=x_2-\lfloor \frac{a}{b}\rfloor \times y_2$

向下迭代，直到$b=0$，此时将$x=1，y=0$返回，向上迭代，即可得到$x_1,y_1$。

通解： 我们刚才求的只是$ax+by=gcd(a,b)$的一组特解，但是题目往往让我们求最小整数解。通过求x的最小整数解为例，假设存在另外的一组解$x_0,y_0$：

1. $ax+by=gcd(a,b)$
2. $y=\frac{gcd(a,b)-ax}{b}$
3. $y_0=\frac{gcd(a,b)-a{x}_0}{b}$
4. 假设$x_0=x+i$
5. $y_0=\frac{gcd(a,b)-ax}{b}-\frac{ai}{b}$
6. $y_0=y-\frac{ai}{b}$
7. 因为$y\in Z,y_0\in Z$
8. $a\ast i/b\in Z$
9. 即$a/b\ast i\in Z$
10. 若$a,b$互质，则$i_{min}=b$
11. 若$a,b$不互质，那么存在$a^{\prime }/b^{\prime }\ast i\in Z$
12. $a^{\prime }=a/gcd(a,b),b^{\prime }=b/gcd(a,b)$

由此，可知$x_i=x+k\ast i,i=b^{\prime }=b/gcd(a,b)$
同理可知，$y_i=y+k\ast j,j=a^{\prime }=a/gcd(a,b)$

代码(递归版）：

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;
	return ;
}

代码（非递归版）：

void exgcd(int a,int b,int &x,int &y)
{
	x = 1, y = 0;
	int x1 = 0, y1 = 1, a1 = a, b1 = b;
	while (b1)
	{
		int q = a1 / b1;
		tie(x, x1) = make_tuple(x1, x - q * x1);
		tie(y, y1) = make_tuple(y1, y - q * y1);
		tie(a1, b1) = make_tuple(b1, a1 - q * b1);
	}
	return ;
}