1. 程式人生 > >51nod 1079中國剩餘定理

51nod 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

中國剩餘定理是用來解決模線性方程組問題的

(S) : \quad \left\{ \begin{matrix} x \equiv a_1 \pmod {m_1} \\ x \equiv a_2 \pmod {m_2} \\ \vdots \qquad\qquad\qquad \\ x \equiv a_n \pmod {m_n} \end{matrix} \right.

求出 x 滿足這個等式




第一種方法:

直接寫

#include<cstdio>
#include<cstring> 
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn=11;

int main()
{
	int n;
	while(~scanf("%d",&n))	
	{
		ll a[maxn],m[maxn];
		ll lcm=1;
		ll ans=0;
		for(int i=0;i<n;i++)
		{
			scanf("%lld %lld",&a[i],&m[i]);
			lcm*=a[i];
		}
		for(int i=0;i<n;i++)
		{
			int j=1;
			ll M=lcm/a[i];
			if(M%a[i]!=0)
				for(j=1;j*M%a[i]!=1;j++);
			ans+=M*m[i]*j;
		}
		printf("%lld\n",ans%lcm);
	}
	return 0;
}

第二種,利用exgcd

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn=105;
int n;
ll m[maxn],p[maxn];
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(!b)
	{
		x=1;
		y=0;
		return a;
	}
	int r=exgcd(b,a%b,y,x);
	y-=x*(a/b);
	return r;
}
ll CRT(ll *m,ll *a)
{
	ll M=1,ans=0;
	for(int i=0;i<n;i++)
		M*=m[i];
	for(int i=0;i<n;i++)
	{
		ll mi=M/m[i],x,y;
			exgcd(mi,m[i],x,y);	//要注意為止不能放反。。。 
		ans=(ans+x*mi*a[i])%M;	
	}
	return (ans+M)%M;
}
int main()
{

	while(~scanf("%d",&n))
	{
		for(int i=0;i<n;i++)
			scanf("%lld %lld",m+i,p+i);
		printf("%lld\n",CRT(m,p));
	}
	return 0;
}