剩余定理

#include<iostream>
using namespace std;
void egcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;
y=0;
return;
}
egcd(b,a%b,x,y);
int xx=x;
x=y;
y=xx-a/b*y;
return ;
}
int CRT(int a[],int m[],int n)
{
int i,M=1,x,y,Mi;
for(i=1;i<=n;i++) M*=m[i];
int res=0;
for(i=1;i<=n;i++)
{
Mi=M/m[i];
egcd(Mi,m[i],x,y);
x%=m[i];            //防止x为负数。当然了，把return res 改成return （res%M+M）%M更加简洁
x=(m[i]+x)%m[i];
res=(res+Mi*x*a[i])%M;
}
return res;
}
int main()
{
int sum,n,a[1000],m[1000];
cin>>n;
for(int i=1;i<=n;i++) cin>>a[i]>>m[i];
sum=CRT(a,m,n);
cout<<sum<<endl;
}