中国剩余定理

中国剩余定理

令$M={\prod }_{i=1}^n{m}_i$, ，$M_i=\frac{M}{m_i}$ ,$t_i\ast M_i\equiv 1(mod{m}_i)$

则
$a_i\ast t_i\ast M_i\equiv a_i(mod{m}_i)$
$a_j\ast t_j\ast M_j\equiv 0(mod{m}_i)$

所以可以构造出答案$X$
$$\mathbf{X}={\prod }_{i=1}^n{\mathbf{M}}_{\mathbf{i}}{\mathbf{t}}_{\mathbf{i}}{\mathbf{a}}_{\mathbf{i}}$$

#include<bits/stdc++.h>
#define I inline
#define int long long
#define RI register int 
using namespace std;
I int read()
{
	RI res=0,f=1;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')f=-f;ch=getchar();}
	while(isdigit(ch)){res=(res<<1)+(res<<3)+(ch&15);ch=getchar();}
	return res*f;
}
int n,a,b,w[15];int x,y;
I void exgcd(int a,int b){
	if(!b){x=1;y=0;return;}
	exgcd(b,a%b);
	int k=x;x=y;y=k-a/b*y;
}
int M=1, m[150000];
int X;
signed main()
{
	n=read(); 
	for(RI i=1;i<=n;i++){m[i]=read();w[i]=read();M*=m[i];}
	for(RI i=1;i<=n;i++){
		exgcd(M/m[i],m[i]);
		X+=(x<0 ? x+m[i]:x)*(M/m[i])*w[i];
	}
	printf("%lld",X%M);
	return 0;
}