浪费我好多时间因为没找着一篇好的教程嘤嘤嘤
推荐一位大神的博客(是真的厉害):dalao
以下是c++版本的程序
#include <iostream>
#include <cstdio>
#define ll __int128
using namespace std;
ll C1,C2,M1,M2,C[1000002],M[1000002],x,y,n,t;
ll gcd(ll x,ll y){
return y==0?x:gcd(y,x%y);
}
ll exgcd(ll a,ll b,ll &x,ll &y){
if(b==0){
x=1,y=0;
return a;
}
int r=exgcd(b,a%b,x,y);
int tmp;
tmp=x,x=y,y=tmp-(a/b)*y;
return r;
}
ll inv(ll a,ll b){
ll r=exgcd(a,b,x,y);
return (x%b+b)%b;
}
int main(){
scanf("%lld",&n);
for(ll i=1;i<=n;i++)
scanf("%lld%lld",&M[i],&C[i]);
for(ll i=2;i<=n;i++){
C1=C[i-1],C2=C[i],M1=M[i-1],M2=M[i];
t=gcd(M1,M2);
M[i]=(M1*M2)/t;
C[i]=(inv(M1/t,M2/t)*(C2-C1)/t)%(M2/t)*M1+C1;
C[i]=(C[i]%M[i]+M[i])%M[i];
}
printf("%lld",C[n]);
return 0;
}