中国剩余定理——一分钟从入门到放弃

Description

给定$n$个限制条件，形如$x$膜$a_i$等于$b_i$，求$x$的最小整数解。

例如，三个限制条件为

3 1
5 1
7 2

(每一行第一个数是$a_i$，第二个是$b_i$)

答案为

16

Solution

这是中国剩余定理的板子题。

直接给出结论: 设$mul={\prod }_{i=1}^n{a}_i$, $m_i=\frac{mul}{a_i}$, $t_i$为在膜$a_i$意义下$m_i$的逆元。

则$x$的一个整数解为${\sum }_{i=1}^n{b}_i{t}_i{m}_i$。

Proving

首先，由于$m_i=\frac{mul}{a_i}$，所以对于任意$j\ne i$，有$m_i$为$a_j$的倍数。

所以，在某一项$b_i{t}_i{m}_i$中，由于$m_i$膜其他的$a_j$都是$0$，所以这一项与其他的限制无关。

又因为$t_i$为$m_i$在膜$a_i$意义下的逆元，所以$m_i{t}_i$膜$a_i$为$1$；所以$m_i{t}_i{b}_i$膜$a_i$为$b_i$。即，对于任意一项$b_i{t}_i{m}_i$，我们既满足了这个限制条件，也没有影响别的限制条件的成立。

$证毕$

Code

模板: 洛谷P1495

#include <bits/stdc++.h>
#define int long long
using namespace std;

int n,mul=1,ans=0,x,y;
int a[20005],b[20005],M[20005],t[20005];

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

signed main()
{
	cin>>n;
	for (int i=1;i<=n;i++)  cin>>a[i]>>b[i];
	for (int i=1;i<=n;i++)  mul*=a[i];
	for (int i=1;i<=n;i++)  M[i]=mul/a[i];
	for (int i=1;i<=n;i++)  exgcd(M[i],a[i],x,y),t[i]=(x+2*a[i])%a[i];
	for (int i=1;i<=n;i++)  ans=(ans+M[i]*t[i]*b[i])%mul;
	cout<<ans<<endl;
	
	return 0;
}

注意点

①一定要开$longlong$。并且用这个东西的时候，千万小心溢出$longlong$，要仔细观察数据范围；
②求逆元的时候，由于$a_i$不一定是质数，所以要用扩欧来求逆元，而不能用费吗小定理(费马小定理满足要求当且仅当模数为质数)。如果$mul$与$a_i$不互质咋办？这个就说明模数不两两互质，那么普通的中国剩余定理就不再适用，得用扩展中国剩余定理。
③在膜意义下的逆元可能是负数，所以要再加上一个$a_i$，详见代码倒数第$6$行；
④最终解出的$x$是一个特解而不是最小整数解，还要膜上$mul$。

时间复杂度$O(n)$，这个问题被中国人巧妙解决。但是，最后再强调一遍，如果模数不两两互质，普通中国剩余定理不再适用。