[excrt] luogu 4777

题目

题目链接：https://www.luogu.com.cn/problem/P4777

代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
int a[100010],b[100010];
int exgcd(int a,int b,int &x,int &y){
	if(b==0){
		x=1;
		y=0;
		return a;
	}
	int gcd=exgcd(b,a%b,x,y);
	int tmp = x;
	x=y;
	y= tmp-a/b*y;
	return gcd;
}
int mul(int a,int b,int mod){
	int res=0;
	while(b){
		if(b&1){
			res=(res+a)%mod;
		}
		b>>=1;
		a=(a+a)%mod;
	}
	return res;
}
signed main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	int n;
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>b[i]>>a[i];	
	}
	int x,y;
	int ans=a[1],M=b[1];
	bool flag=0;
	for(int i=2;i<=n;i++){
		int aa=M,bb=b[i];
		int c=(a[i]-ans%bb+bb)%bb;
		int gcd=exgcd(aa,bb,x,y);
		int gd=bb/gcd;
//		if(c%gcd!=0){//判断存不存在答案 
//			flag=1;
//			break;
//		}
		x=mul(x,c/gcd,gd);
		ans+=x*M;
		M*=gd;
		ans=(ans%M+M)%M;
	}
	cout<<ans<<endl; 
    return 0;
}