扩展欧几里得

扩展欧几里得可以求出
$ax+by=gcd(a,b)$的一组特解

证明：
当$b=0$时，显然有$x=1,y=0$这一组解，但其实$y$并非只能取零，取任何值都可以。

若$b\ne 0$，因为$ax+by=gcd(a,b)=gcd(b,a\%b)=bx+(a\%b)y$
$=bx+(a-\lfloor a/b\rfloor \times b)\times y$
$=bx+ay-\lfloor a/b\rfloor by$
$=ay+b(x-\lfloor a/b\rfloor y)$
即存在另一组解
$y$和$x-\lfloor a/b\rfloor y$
令$x^{{}^{\prime }}=y,y^{{}^{\prime }}=x-\lfloor a/b\rfloor y$，则有$a{x}^{{}^{\prime }}+b{y}^{{}^{\prime }}=gcd(a,b)$
即$x^{{}^{\prime }}$和 $y^{{}^{\prime }}$是一组解

证毕

因此当$b\ne 0$时，我们可以递归调用欧几里得算法，每次回溯的时候令$x=y,y=x-\lfloor a/b\rfloor y$，最后就会得到一组可行解，这个过程有两种表示方法。

// 1
ll exgcd(ll a, ll b, ll& x, ll& y){
	if (b == 0){
		x = 1;
		y = 0;//等于多少无所谓
		return a;
	}
	ll t = exgcd(b, a % b, y, x);
	y -= a / b * x;
	return t;
}

// 2
ll exgcd(ll a, ll b, ll& x, ll& y){
	if (b == 0){
		x = 1;
		y = -534;
		return a;
	}
	ll t = exgcd(b, a % b, x, y);
	ll z = x;
	x = y;
	y = z - (a / b) * y;
	return t;
}