51nod 1079中国剩余定理

最新推荐文章于 2019-10-04 21:33:06 发布

陶鸿杰

最新推荐文章于 2019-10-04 21:33:06 发布

阅读量217

点赞数

分类专栏：数学

本文链接：https://blog.csdn.net/Hpuer_Random/article/details/80672834

版权

数学专栏收录该内容

45 篇文章 0 订阅

订阅专栏

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示例

Output示例

中国剩余定理是用来解决模线性方程组问题的

$(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 满足这个等式

wiki传送

第一种方法：

直接写

#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;
}