比赛 EYOI暨SBOI暑假快乐赛2nd 评测结果 AAAAAAAAAA
题目名称 曹冲养猪 最终得分 100
用户昵称 张恒畅 运行时间 0.000 s
代码语言 C++ 内存使用 0.00 MiB
提交时间 2022-06-26 11:55:40
显示代码纯文本
    #include<iostream>
    #include<cstring>
    #include<cstdio>
    using namespace std;
    typedef long long L;
    int m;
    L a[21],mod[21];
    inline L gcd(L a,L b){
    	return a%b?gcd(b,a%b):b;
    }
    inline L ext_gcd(L a,L b,L &x,L &y){
    	if(b==0){
    		x=1;
    		y=0;
    		return a; 
    	}
    	L gcd(ext_gcd(b,a%b,x,y));
    	L tmp(x);
    	x=y;
    	y=tmp-(a/b)*y;
    	return gcd;
    }
    inline L ny(L a,L b){
    	L x,y;
    	L gcd(ext_gcd(a,b,x,y));
    	if(gcd!=1)
    		return -1;
    	return (x%b+b)%b;
    }
    inline bool merge(L a1,L m1,L a2,L m2,L &a3,L &m3){
    	L d(gcd(m1,m2));
    	L c=a2-a1;
    	if(c%d)
    		return false;
    	c=(c%m2+m2)%m2;
    	m1/=d;
    	m2/=d;
    	c/=d;
    	c*=ny(m1,m2);
    	c%=m2;
    	c*=m1*d;
    	c+=a1;
    	m3=m1*m2*d;
    	a3=(c%m3+m3)%m3;
    	return true;
    }
    L CRT(L a[],L m[],int n){
    	L a1(a[1]),m1(m[1]);
    	for(int i=2;i<=n;i++){
    		L a2(a[i]),m2(m[i]),a3,m3;
    		if(!merge(a1,m1,a2,m2,a3,m3))
    			return -1;
    		a1=a3;
    		m1=m3;
    	}
    	return (a1%m1+m1)%m1;
    }
    inline int gg(){
    	freopen("ccyz.in","r",stdin);
    	freopen("ccyz.out","w",stdout);
    	scanf("%d",&m);
    	for(int i=1;i<=m;i++)
    		scanf("%lld%lld",&mod[i],&a[i]);
    	printf("%lld",CRT(a,mod,m));
    }
    int k(gg());
    int main(){;}