扩展中国剩余定理模板

题目链接：点击这里

题目大意：
给定 $n$ 组非负整数 $a_i,b_i$ ，求解关于 $x$ 的方程组的最小非负整数解。
$$\{\begin{array}{c}x\equiv b_1\mathrm{m}\mathrm{o}\mathrm{d}{m}_1\\ x\equiv b_2\mathrm{m}\mathrm{o}\mathrm{d}{m}_2\\ ...\\ x\equiv b_n\mathrm{m}\mathrm{o}\mathrm{d}{m}_n\end{array}$$
与中国剩余定理不同的是，此时的 $m_1,m_2,...,m_n$ 两两不必互质

题目分析：
中国剩余定理求出的 $m_1,m_2,...,m_n$ 两两互质的同余方程组的解为：
$$x_0\equiv \sum _{i=1}^n{b}_i{M}_{i}inv(M_i,m_i)\mathrm{m}\mathrm{o}\mathrm{d}n$$
因为 $inv(M_i,m_i)$ 在 $M_i,m_i$ 不互质时是不存在的，所以中国剩余定理只能解决 $m_1,m_2,...,m_n$ 两两互质的情况，对于 $m_1,m_2,...,m_n$ 两两不一定互质的方程很难在 $crt$ 的基础上进行修改来使之完备，这是 $excrt$ 为了解决这种问题边应运而生了。
我们用数学归纳法，假设我们已经求出来了前 $k-1$ 个方程的解 $x$ ，令 $M={\sum }_{i=1}^{k-1}{m}_i$ ，则 $x+aM(a\in Z)$ 也是前 $k-1$ 个方程的解。
在计算第 $k$ 个方程时，我们求 $t$ 使之满足：
$$x+tM\equiv b_i\mathrm{m}\mathrm{o}\mathrm{d}n$$
等价于求：
$$tM\equiv b_i-x\mathrm{m}\mathrm{o}\mathrm{d}n$$
转化为求不定方程：
$$Mt+ny=b_i-x$$
此式可通过扩欧判断其是否有解，若有解， $x^{{}^{\prime }}=x+tM$ 就是前 $k$ 个方程的解，通过这种迭代方法从求出的前两个方程的解就可以求出对全部方程的解了

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<unordered_map>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
ll read()
{
	ll res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 1e5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
ll n,b[maxn],m[maxn],mul = 1;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b == 0)
	{
		x = 1;
		y = 0;
		return a;
	}
	ll d = exgcd(b,a%b,x,y);
	int tmp = x;
	x = y;
	y = tmp-y*(a/b);
	return d;
}
ll ksc(ll a,ll b,ll mod)
{
	return (a*b-(ll)((long double)a/mod*b)*mod+mod)%mod;
}
ll excrt(ll b[],ll m[],ll n)
{
	ll ans = b[1],M = m[1],d,x,y;
	for(int i = 2;i <= n;i++)
	{
		d = exgcd(M,m[i],x,y);
		ll c = ((b[i]-ans)%m[i]+m[i])%m[i];
		if(c%d) return -1; // 无解 
		x = ksc(c/d,x,m[i]/d);//x = c/d*x%(m[i]/d);
		ans += M*x;
		M = M/d*m[i];
		ans = (ans%M+M)%M;
	}
	return (ans%M+M)%M;
}
int main()
{
	n = read();
	for(int i = 1;i <= n;i++)
		m[i] = read(),b[i] = read();
	printf("%lld\n",excrt(b,m,n));
	return 0;
}