社区讨论

0tps求条 玄关 在线等 急!!

P4777【模板】扩展中国剩余定理(EXCRT)参与者 2已保存回复 1

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
1 条
当前快照
1 份
快照标识符
@mjv1bpp1
此快照首次捕获于
2026/01/01 13:59
2 个月前
此快照最后确认于
2026/01/03 18:55
2 个月前
查看原帖
CPP
#include <bits/stdc++.h>
using namespace std;
#define int long long
int n;
const int N=1e5+10;
int a[N],m[N];
int exgcd(int a,int b,int &x,int &y){
	if(b==0){
		x=1;
		y=0;
		return a;
	}
	int d=exgcd(b,a%b,x,y);
	y-=a/b*x;
	return d;
}
int mul(int a,int b,int m){
	int res=0;
	while(b>0){
		if(b&1) res=(res+a)%m;
		a=(a+a)%m;
		b>>=1;
	}
	return res;
}
int EXcrt(){
	int x,y;
	int m1=m[1],a1=a[1];
	int ans=(a1%m1+m1)%m1;
	for(int i=2;i<=n;i++){
		int m2=m[i],a2=a[i];
		int a=m1,b=m2,c=(a2-a1%m2+m2)%m2;
		int d=exgcd(a,b,x,y);
		if(c%d!=0) return -1;
		x=mul(x,c/d,b/d);
		ans=a1+x*m1;
		m1=m2/d*m1;
		ans=(ans%m1+m1)%m1;
		a1=ans;
	}
	return ans;
}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>m[i]>>a[i];
	}
	cout<<EXcrt();
	return 0;
}

回复

1 条回复,欢迎继续交流。

正在加载回复...