扩展中国剩余定理

在介绍扩展中国剩余定理之前，让我们先介绍一下中国剩余定理。

中国剩余定理

中国剩余定理是用来求解一堆线性同余方程组的通解的算法。
该算法的时间复杂度是$O(方程的个数)$。
$$\{\begin{array}{rc}x\equiv a_1& (mod{m}_1)\\ x\equiv a_2& (mod{m}_2)\\ x\equiv a_3& (mod{m}_3)\\ .........& \\ x\equiv a_n& (mod{m}_n)\end{array}$$
其中，$m_1,m_2,...,m_n$两两互质。
对于这个方程，我们很容易构造出其通解。
1.计算模数累积，令$M={\prod }_{i=1}^n$
2.枚举$i$从$1$到$n$，将$a_i\cdot \frac{M}{m_i}\cdot inv\left(\frac{M}{m_i}，m_i\right)$累加到答案中，注意边累加边$modM$。
其中，$inv\left(a,b\right)$表示$a\mathrm{mod}b$的乘法逆元，可以用$exgcd$求出。

#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define pii pair<int,int>
using namespace std;

const double eps = 1e-10;
const double pi = acos(-1.0);

int n;
ll a[30],m[30],M = 1;
/*
算法步骤:
计算模数积，然后遍历每个方程，让答案累加a[i]*M/m[i]*inv(M/m[i],m[i]);
其中inv(a,b)表示a在modb意义下的逆元。
注意，中国剩余定理必须要满足模数两两互质。
*/

template<class T> T gcd(T a, T b){
	return !b?a:gcd(b,a%b);
}

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

ll inv(ll a, ll b){
	ll x,y;
	exgcd(a,b,x,y);
	return (x%b+b)%b;
}

void solve(){
    scanf("%d",&n);
	for(int i = 1; i <= n; i++){
		scanf("%lld%lld",&m[i],&a[i]);
		M *= m[i];
	}
	ll ans = 0;
	for(int i = 1; i <= n; i++){
		ans = (ans + a[i]*M/m[i]*inv(M/m[i],m[i]))%M;
	}
	printf("%lld\n",(ans%M + M)%M);
}

int main()
{
	solve();
	return 0;
}

扩展中国剩余定理

$$\{\begin{array}{rc}x\equiv a_1& (mod{m}_1)\\ x\equiv a_2& (mod{m}_2)\\ x\equiv a_3& (mod{m}_3)\\ .........& \\ x\equiv a_n& (mod{m}_n)\end{array}$$
这时的方程模数不再满足两两互质，甚至还有可能无解。
还是老样子，我们先考虑简单情况。
先看前两个方程组:
$$\{\begin{array}{rc}x\equiv a_1& (mod{m}_1)\\ x\equiv a_2& (mod{m}_2)\end{array}$$
它等价于
$$\{\begin{array}{r}x=a_1+k_1{m}_1\\ x=a_2+k_2{m}_2\end{array}$$
于是消去$x$，我们就得到了一个方程:
$$k_1{m}_1-k_2{m}_2=a_2-a_1$$
显然，这是个简单的二元不定方程。
如果$gcd(m1,m2)\nmid a_2-a_1$，那么方程组无解。
否则我们可以用$exgcd$求出一组解。
设这组解为$k_1^{\ast },k_2^{\ast }$，那么我们可以得到这组方程组的通解：
$$\{\begin{array}{r}{k}_1=k_1^{\ast }+\frac{m_2}{gcd(m_1,m_2)}t\\ \\ k_2=k_2^{\ast }+\frac{m_1}{gcd(m_1,m_2)}t\end{array}$$
将第一个方程带入，得
$$\begin{array}{rl}x=& a_1+k_1^{\ast }{m}_1+\frac{m_1{m}_2}{d}t\\ =& a_1+k_1^{\ast }{m}_1+t[m_1,m_2]\end{array}$$
令$$x_0=a_1+k_1^{\ast }{m}_1$$
那么原来的两个方程可以转化为一个方程
$$x=x_0+t[m_1,m_2]$$
这样，我们可以通过n-1次循环将n个方程合并为1个方程，最后通过解最后一个方程得到整个方程组的解。
算法步骤:
1.初始化$a_1$，$m_1$;
2.对于第$i$个方程，令$a_2=a_i$，$m_2=m_i$，对其用拓展欧几里得算法:$$d=exgcd(m_1,m_2,k_1,k_2)$$
如果 $$d\nmid (a_2-a_1)$$那么无解；
否则$$k_1\ast =\frac{a_2-a_1}{d}$$
这里我们只需要$k_1$的值，故不用对$k_2$进行任何操作。
4.为了不让结果溢出且恒为正，我们需要做一些预处理:令$t=\frac{m_2}{d}$，然后$$k_1=(k_1\%t+t)\%t$$
令$$a_1=a_1+k_1{m}_1,m_1=\left|lcm(m_1,m_2)\right|$$。
重复循环$n-1$次，我们可以得到一个方程
$$x=a_1+t{m}_1$$
所得结果就是$\left(a_1\%m_1+m_1\right)\%m_1$。
各个变量关系非常复杂，建议在考场上实现代码的是由用手推一下式子。

#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define pii pair<int,int>
using namespace std;

const double eps = 1e-10;
const double pi = acos(-1.0);
const int maxn = 1e6 + 10;

int n;

inline ll mul(ll a,ll b,ll mod){//O(1)取模快速乘，不会爆long long
	return (a*b-(ll)((long double)a/mod*b)*mod+mod)%mod;
}

ll exgcd(ll a, ll b, ll& x, ll& y){
	if(!b){
		x = 1, y = 0;
		return a;
	}
	ll d = exgcd(b,a%b,y,x);
	y -= a/b*x;
	return d;
}

void solve(){
	ll a1,m1;
	scanf("%d",&n);
	scanf("%lld%lld",&m1,&a1);
	bool ok = 1;
	for(int i = 2; i <= n; i++){
		ll a2,m2,k1,k2;
		scanf("%lld%lld",&m2,&a2);
		ll d = exgcd(m1,m2,k1,k2);
		if((a2-a1)%d) ok = 0;
		else{
			k1 = mul(k1,(a2-a1)/d,m2/d);//这个地方必须要用取模快速乘
			a1 = a1+k1*m1;
			m1 = abs(m1/d*m2);
		}
	}
	if(ok)printf("%lld\n",(a1%m1+m1)%m1);
}

int main()
{
	solve();
	return 0;
}