当
a
1
,
a
2
a_{1},a_{2}
a1,a2互质的时候
我们可以使用中国剩余定理合并
即
∑
i
b
i
∗
(
M
a
i
∗
I
n
v
(
M
a
i
)
)
\sum_{i}b_{i}*(\frac{M}{a_{i}}*Inv_{(\frac{M}{a_{i}})})
∑ibi∗(aiM∗Inv(aiM))
其中
M
=
∏
i
a
i
M=\prod_{i}a_{i}
M=∏iai
但是不同于别的定律
从CRT到EXCRT是完全不同的
观察下面
T
T
T
m
o
d
mod
mod
a
1
a_{1}
a1=
b
1
b_{1}
b1
T
T
T
m
o
d
mod
mod
a
2
a_{2}
a2=
b
2
b_{2}
b2
列出以下式子
a
1
∗
x
1
+
b
1
=
a
2
∗
x
2
+
b
2
a_{1}*x_{1}+b_{1}=a_{2}*x_{2}+b_{2}
a1∗x1+b1=a2∗x2+b2
转化为
a
1
∗
x
1
−
a
2
∗
x
2
=
b
2
−
b
1
a_{1}*x{1}-a_{2}*x_{2}=b_{2}-b_{1}
a1∗x1−a2∗x2=b2−b1
所以可以EXGCD
解出最小正整数解
x
1
x_{1}
x1后回带求出
T
T
T
这时有着一个新的同余式
T
T
T
m
o
d
mod
mod
L
C
M
(
a
1
,
a
2
)
LCM(a_{1},a_{2})
LCM(a1,a2) =
b
S
o
l
b_{Sol}
bSol
递推计算
n
−
1
n-1
n−1次即答案
#include<bits/stdc++.h>
using namespace std;
typedef int INT;
#define int long long
const int N=1e5+10;
inline void read(int &x){
x=0;
int f=1;
char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9'){
x=x*10+ch-'0';
ch=getchar();
}
x*=f;
}
struct Node{
int ret,mod;
}P[N];
int n;
int EXGCD(int &x,int &y,int a,int b){
if(!b){
x=1;
y=0;
return a;
}
int GCD=EXGCD(x,y,b,a%b);
int tmp=x;
x=y;
y=tmp-(a/b)*y;
return GCD;
}
int mul(int x,int y,int mod){
int ret=0;
while(y){
if(y&1)ret=(ret+x)%mod;
x=(x+x)%mod;
y/=2;
}
return ret;
}
void EXCRT(int tot){
int ret=P[1].ret;
int mod=P[1].mod;
for(int i=2;i<=tot;++i){
int x,y;
int GCD=EXGCD(x,y,mod,P[i].mod);
int now=((P[i].ret-ret)%P[i].mod+P[i].mod)%P[i].mod;
x=(x%P[i].mod+P[i].mod)%P[i].mod;
x=mul(x,now/GCD,P[i].mod)%P[i].mod;
ret=(ret+mul(x,mod,mod/GCD*P[i].mod));
mod=mod/GCD*P[i].mod;
ret=ret%mod;
}
cout<<ret;
}
INT main(){
// freopen("test.in","r",stdin);
read(n);
for(int i=1;i<=n;++i){
read(P[i].mod);
read(P[i].ret);
}
EXCRT(n);
}