扩展中国剩余定理

废话

这玩意跟中国剩余定理没有半毛钱关系，所以所以没学过中国剩余定理的完全不用担心。

当然还是建议先学一学中国剩余定理。

介绍

这是一个用来求形如下面这样的同余方程组的解的算法：
$$\{\begin{array}{l}x\equiv a_1(\mathrm{mod}{m}_1)\\ x\equiv a_2(\mathrm{mod}{m}_2)\\ ⋮\\ x\equiv a_n(\mathrm{mod}{m}_n)\end{array}$$

其中， $m_1,m_2,...,m_n$ 之间不一定互质。

做法

很神奇，网上有两个版本，有两种不同的做法。（并不知道哪一个是正版……）

做法 $1$

假设我们已经求出了前 $p$ 个方程组的解 $x_1$。

不妨设 $M=\prod _{i=1}^p{m}_i$，那么前 $p$ 个方程组的通解为 $x_1+kM(k\in Z)$。

那么要满足第 $p+1$ 个方程，相当于要找到一个合法的 $k$，满足：
$$x_1+kM\equiv a_{p+1}(\mathrm{mod}{m}_{p+1})$$

相当于
$$x_1+kM+L{m}_{p+1}=a_{p+1}(L\in Z)(\mathrm{mod}{m}_{p+1})$$

移项得
$$kM+L{m}_{p+1}=a_{p+1}-x_1(\mathrm{mod}{m}_{p+1})$$

用扩展欧几里得算法求出 $k,L$ 即可，那么对于前 $p+1$ 个方程的解就是 $x_1+kM$。

但是这种做法有一个限制，就是 $M$ 要在 $longlong$ 范围内，而有些恶心题就不保证。

有一个小优化就是稍稍变一下 $M$ 的定义，让 $M=lc{m}_{i=1}^p{m}_i$，这样可以稍稍减小一点 $M$。

模板题传送门

代码如下：

#include <cstdio>
#include <cstring>
#define ll long long
#define maxn 100010

int n;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)return x=1,y=0,a;
    ll d=exgcd(b,a%b,y,x);
    return y-=a/b*x,d;
}
ll a[maxn],m[maxn];
ll ksc(ll x,ll y,ll mod)
{
	ll re=0,tot=x;
	while(y)
	{
		if(y&1ll)re=(re+tot)%mod;
		tot=(tot+tot)%mod;
		y>>=1;
	}
	return re;
}
ll excrt()
{
	ll M=m[1],x=a[1];
	for(int i=2;i<=n;i++)
	{
		ll k,L,c=(a[i]-x%m[i]+m[i])%m[i];
		ll gcd=exgcd(M,m[i],k,L),bg=m[i]/gcd;
		k=ksc(k,c/gcd,bg);
		x+=k*M;
		M*=bg;
		x=(x%M+M)%M;
	}
	return x;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	scanf("%lld %lld",&m[i],&a[i]);
	printf("%lld\n",excrt());
}

做法 $2$

这种做法就是大力将两个同余式和并为一个，然后大力合并所有同余式。

假设现在有两个同余式：
$$\{\begin{array}{l}x\equiv a_1(\mathrm{mod}{m}_1)\\ x\equiv a_2(\mathrm{mod}{m}_2)\end{array}$$

那么有
$$\{\begin{array}{l}x=a_1+k_1{m}_1\\ x=a_2+k_2{m}_2\end{array}$$

则有
$$\begin{array}{rl}{a}_1+k_1{m}_1& =a_2+k_2{m}_2\\ k_1{m}_1& =a_2-a_1+k_2{m}_2\end{array}$$

设 $g=\gcd (m_1,m_2)$，将柿子除以 $g$，则有：
$$\begin{array}{rl}\frac{k_1{m}_1}{g}& =\frac{a_2-a_1}{g}+\frac{k_2{m}_2}{g}\\ \frac{k_1{m}_1}{g}& \equiv \frac{a_2-a_1}{g}(\mathrm{mod}\frac{m_2}{g})\\ k_1& =\frac{a_2-a_1}{g}\times inv(\frac{m_1}{g})+y\times \frac{m_2}{g}\end{array}$$

其中，$inv$ 是指逆元。

将 $k_1$ 带回去，有
$$\begin{array}{rl}x& =a_1+m_1\times \left(\frac{a_2-a_1}{g}\times inv(\frac{m_1}{g})+y\times \frac{m_2}{g}\right)\\ x& =a_1+m_1\times \frac{a_2-a_1}{g}\times inv(\frac{m_1}{g})+y\times \frac{m_1{m}_2}{g}\\ x& \equiv a_1+m_1\times \frac{a_2-a_1}{g}\times inv(\frac{m_1}{g})(\mathrm{mod}\frac{m_1{m}_2}{g})\end{array}$$

于是两条同余式合并成了一条。

但是要注意，里面 $\frac{a_2-a_1}{g}\times inv(\frac{m_1}{g})$ 这一部分依然要在模 $\frac{m_2}{g}$ 意义下计算，但是前面的 $\times m_1$ 在模 $\frac{m_1{m}_2}{g}$ 意义下计算。

当 $g$ 不能整除 $a_2-a_1$ 时，该方程组无解。

代码如下：

#include <cstdio>
#include <cstring>
#define ll long long
#define maxn 100010

int n;
ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
void exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0){x=1,y=0;return;}
	exgcd(b,a%b,y,x);
	y-=a/b*x;
}
ll inv(ll a,ll mod)
{
	ll x,y;
	exgcd(a,mod,x,y);
	return (x%mod+mod)%mod;
}
ll ksc(ll x,ll y,ll mod)
{
	ll re=0,tot=x;
	y=(y%mod+mod)%mod;
	while(y)
	{
		if(y&1ll)re=(re+tot)%mod;
		tot=(tot+tot)%mod;
		y>>=1;
	}
	return re;
}
ll a[maxn],m[maxn];
ll excrt()
{
	ll a1=a[1],m1=m[1],a2,m2,g;
	for(int i=2;i<=n;i++)
	{
		a2=a[i],m2=m[i];
		g=gcd(m1,m2);
		if((a2-a1)%g!=0)return -1;
		a1+=ksc(inv(m1/g,m2/g),(a2-a1)/g,m2/g)%(m2/g)*m1;
		m1*=m2/g;
		a1=(a1%m1+m1)%m1;
	}
	return a1;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	scanf("%lld %lld",&m[i],&a[i]);
	printf("%lld",excrt());
}