扩展欧几里得算法

前置知识：

1.给予两个整数a，b。必存在整数x，y使得ax+by=gcd(a,b)。
2.0和任何一个数的最大公约数是这个数的本身。
3.辗转相除：
$$\begin{gathered} gcd(a,b)=gcd(b,amodb) \\ \because amodb=r\therefore a=kb\ast r; \\ 设gcd(a,b)=c\therefore a|c,b|c \\ \therefore bk|c \\ \because a|c\therefore (bk+r)|c,r|c \\ 由此得出：gcd(a,b)=gcd(b,amodb) \end{gathered}$$

推理：

1.当b=0，gcd(a,b) = a。所以x = 1, y = 0 --> ax=a;
2.当 a != 0, b != 0:
由辗转相除法可知：
$$\begin{gathered} gcd(a,b)=gcd(b,amodb) \\ \therefore ax+by=b{x}^{\prime }+(amodb)y^{\prime } \\ \therefore ax+by=b{x}^{\prime }+(a-a/b\ast b)y^{\prime } \\ \therefore ax+by=a{y}^{\prime }+b(x^{\prime }-a/b\ast y^{\prime }) \\ 对比系数：x=y^{\prime }y=x^{\prime }-a/b\ast y^{\prime } \end{gathered}$$
也就是说，如果我们知道了x’,y’就可以倒退得到x, y然而由辗转相除可知，经gcd(a,b) = gcd(b, a%b)如此反复b将变为0，然后利用b = 0时，x = 1,y = 0逐步倒推，最终得到ax+by=gcd(a, b)的一组解。

代码：

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

例题

Acwing 五指山

#include<bits/stdc++.h>

#define IOS ios::sync_with_stdio(false);cin.tie(nullptr)
#define int long long 
#define endl "\n"

using namespace std;

const int N = 2e5 + 10;

int exgcd(int a, int b, int& xx, int& yy)
{
	if (!b)
	{
		xx = 1, yy = 0;
		return a;
	}

	int g = exgcd(b, a % b, yy, xx);

	yy -= a / b * xx;
	return g;
}

void solve()
{
	int n, d, x, y;
	cin >> n >> d >> x >> y;

	int  xx, yy;
	int tmp = exgcd(n, d, xx, yy);
	if ((y - x) % tmp)
	{
		cout << "Impossible" << endl;
		return;
	}
	yy *= (y - x) / tmp;
	n /= tmp;
	cout << (yy % n + n) % n << endl;
}

signed main()
{
	IOS;
	int T;
	cin >> T;
	while (T--)
		solve();
	return 0;
}

如何根据已经求得的一组整数解求得其他的整数解：

x, y是我们根据ax+by=gcd(a,b)得一组通解，现在我们设新得解为x‘，y’;那么x’ = x+s1, y‘ = y-s2，则有
$$\begin{gathered} a(x+s_1)+b(y-s_2)=gcd(a,b) \\ \therefore ax+a{s}_1+by-b{s}_2=gcd(a,b) \\ \because ax+by=gcd(a,b) \\ \therefore a{s}_1=b{s}_2 \end{gathered}$$
让a，b同时除以一个较大的数（即a和b的最大公约数）
$$\begin{gathered} \frac{s_1}{s_2}=\frac{\frac{b}{gcd(a,b)}}{\frac{a}{gcd(a,b)}} \\ \therefore s_1=k\ast \frac{b}{gcd(a,b)},s_2=k\ast \frac{a}{gcd(a,b)} \end{gathered}$$
由此我们可以得到ax+by = gcd(a, b)
$$\begin{gathered} x=x_0+k\ast \frac{b}{gcd(a,b)} \\ y=y_0-k\ast \frac{a}{gcd(a,b)} \\ {x}_0,y_0是EXGCD求得的一组解 \end{gathered}$$

其他：

以上是求ax+by=gcd的通解的推导，如果我们想求ax+by=c的通解，首先我们要判断c是否能整除gcd(a, b)，如果不能整除代表此式子无解。

然后我们根据exgcd求出ax+by=gcd(a,b)一组解x’, y’，然后我们再求出原式的一组解
$$\begin{gathered} x=x^{\prime }\ast \frac{c}{gcd(a,b)}\ast k\ast \frac{b}{gcd(a,b)} \\ y=y^{\prime }\ast \frac{c}{gcd(a,b)}\ast k\ast \frac{a}{gcd(a,b)} \end{gathered}$$