【模板】擴展中國剩餘定理

模板題

推薦博客

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
#define N 1005
inline LL exgcd(LL a,LL b,LL &x,LL &y)
{
	LL d=a;
	if(b>0)
	{
		d=exgcd(b,a%b,y,x);
		y-=x*(a/b);
	}
	else
		x=1,y=0;
	return d;
}
LL m[N],r[N],M,R,x,y;
int n;
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%lld%lld",&m[i],&r[i]);
	M=m[1],R=r[1];
	for(int i=2;i<=n;i++)
	{
		LL d=exgcd(M,m[i],x,y);
		LL c=r[i]-R;
		if(c%d)
		{
			printf("-1");
			return 0;
		}
		exgcd(M/d,m[i]/d,x,y);
		x=(x*c/d)%(m[i]/d);
		R+=x*M;
		M=M/d*m[i];
		R%=M;
	}
	printf("%lld",(R+M)%M);
}