POJ 2891 中国剩余定理

除数不会两两互质,不能套用中国余数定理~

只能不断利用线性模方程把两个式子合并


#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <string>
#include <cmath>
#include <stack>
#include <queue>
using namespace std;
#define ll long long
ll ext_gcd(ll a,ll n,ll &x,ll &y)
{
	if(n==0)
	{
		x=1;y=0;
		return a;
	}
	ll d=ext_gcd(n,a%n,x,y);
	ll x1=x,y1=y;
	x=y1;
	y=x1-a/n*y1;
	return d;
}
int main ()
{
	int k;ll a1,r1,a2,r2,x,y,d,e;
	while(scanf("%d",&k)!=EOF)
	{
		scanf("%lld%lld",&a1,&r1);
		bool f=true;
		for(int i=2;i<=k;++i)
		{
			scanf("%lld%lld",&a2,&r2);
			if(!f)
			{
				continue;
			}
			d=ext_gcd(a1,a2,x,y);
			if((r2-r1)%d)
				f=false;  // no solution
			else 
			{
				e=((r2-r1)/d*x%a2+a2)%(a2/d);
				r1=e*a1+r1;
				a1=a1/d*a2;
			}
		}
		if(f)
			printf("%lld\n",r1);
		else printf("-1\n");
	}
	return 0;
}


你可能感兴趣的:(POJ 2891 中国剩余定理)