51Nod 1079 中国剩余定理（中国剩余定理模板题）

点击转到

1079 中国剩余定理 

基准时间限制：1 秒 空间限制：131072 KB 分值: 0 难度：基础题

 收藏

 关注

一个正整数K，给出K Mod 一些质数的结果，求符合条件的最小的K。例如，K % 2 = 1, K % 3 = 2, K % 5 = 3。符合条件的最小的K = 23。

Input

第1行：1个数N表示后面输入的质数及模的数量。（2 <= N <= 10)
第2 - N + 1行，每行2个数P和M，中间用空格分隔，P是质数，M是K % P的结果。（2 <= P <= 100, 0 <= K < P)

Output

输出符合条件的最小的K。数据中所有K均小于10^9。

Input示例

3
2 1
3 2
5 3

Output示例

23

就是刚接触中国剩余定理时的样例题，套用模板即可。一开始定义的变量都是整形，但是到第14样例以后就WA，第14个样例内容如下：

8
73 44
11 4
2 0
3 0
37 11
7 0
23 0
71 59

此时，输出sum的话，发现整形变量并不足以存储：

所以，还是定义成long long类型：

/*#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y)
{
	
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	int temp=exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-a/b*y;
	return temp;
}
int chinese_Remainder(int r[],int prime[],int n)
{
	int d,x,y,m;
	int lcm=1;
	int sum=0;
	for(int i=0;i<n;i++)
	{
		lcm=lcm*prime[i];
	}
	for(int i=0;i<n;i++)
	{
		m=lcm/prime[i];
		d=exgcd(prime[i],m,x,y);
		sum=(sum+y*m*r[i])%lcm;
	}
	return (lcm+sum%lcm)%lcm;
} 
int main()
{
	int n;
	int prime[16];
	int r[16];
	cin>>n;
	for(int i=0;i<n;i++)
	{
		scanf("%d%d",&prime[i],&r[i]);
	}
	printf("%d\n",chinese_Remainder(r,prime,n));
}*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	ll temp=exgcd(b,a%b,x,y);
	ll t=x;
	x=y;
	y=t-a/b*y;
	return temp;
}
ll chinese_Remainder(ll r[],ll prime[],ll n)
{
	ll d,x,y,m;
	ll lcm=1;
	ll sum=0;
	for(ll i=0;i<n;i++)
	{
		lcm=lcm*prime[i];
	}
	for(ll i=0;i<n;i++)
	{
		m=lcm/prime[i];
		d=exgcd(prime[i],m,x,y);
		sum=(sum+y*m*r[i])%lcm;
		//cout<<"sum="<<sum<<endl;
	}
	return (lcm+sum%lcm)%lcm;
} 
int main()
{
	ll n;
	ll prime[16];
	ll r[16];
	scanf("%lld",&n);
	for(int i=0;i<n;i++)
	{
		scanf("%lld%lld",&prime[i],&r[i]);
	}
	printf("%lld\n",chinese_Remainder(r,prime,n));
	return 0;
}