此模板封装了中国剩余定理等函数,支持同余方程插入。
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define FOR(i,x,y) for(int i=(x);i<=(y);i++)
#define DOR(i,x,y) for(int i=(x);i>=(y);i--)
typedef long long LL;
using namespace std;
struct CRT //Chinese Remainder Theorem
{
LL A,P; //一个数模P为A
CRT(){A=0,P=1;}
LL gcd(LL a,LL b){return !b?a:gcd(b,a%b);} //欧几里得(最大公约数)
void exgcd(LL a,LL b,LL &x,LL &y) //扩展欧几里得(求ax+by=gcd(a,b)x的一个特解)
{
if(b==0){x=1,y=0;return;}
exgcd(b,a%b,y,x);y-=a/b*x;
return;
}
LL Exgcd(LL a,LL b,LL c) //求ax+by=c x的最小非负整数解
{
LL g=gcd(a,b);
if(c%g)return -1;
a/=g,b/=g,c/=g;
LL x,y;
exgcd(a,b,x,y);
return ((x*c%b)+b)%b;
}
void kase_insert(CRT x)
{
LL K=Exgcd(P,x.P,x.A-A);
//if(K==-1)...
A+=P*K;
P=(P*x.P)/gcd(P,x.P);
return;
}
void Read()
{
scanf("%lld%lld",&P,&A);
return;
}
LL ask(){return A?A:P;}
}crt;
int main()
{
int m;
scanf("%d",&m);
FOR(i,1,m)
{
CRT x;
x.Read();
crt.kase_insert(x);
}
printf("%lld\n",crt.ask());
return 0;
}