[中国剩余定理+快速乘]

题目

题目连接：https://www.luogu.com.cn/problem/solution/P3868

思路

模板题
注意快速乘不是提速的 而是减小数据范围的，因为这道题最大的时候会达到10^32 所以要用快速乘

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
int a[15],b[15];
void exgcd(int a,int b,int &x,int &y){
	if(b==0){
		x=1;
		y=0;
		return;
	}
	exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-a/b*y;
}
int quick_mul(int a,int b,int mod){//快速乘 
	int res=0;
	while(b){
		res=(res+b%2*a)%mod;
		a=(a<<1)%mod;
		b>>=1;
	}
	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>>a[i];
	}
	for(int i=1;i<=n;i++){
		cin>>b[i];
	} 
	int M=1;
	for(int i=1;i<=n;i++){
		M=M*b[i];
	}
	int res=0;
	for(int i=1;i<=n;i++){
		int m=M/b[i];
		int x,y;
		exgcd(b[i],m,x,y);//y为逆元
		y=(y%b[i]+b[i])%b[i];//取模
		res=(res+quick_mul(y,quick_mul(m,(a[i]+M)%M,M),M)+M)%M;
	}
	cout<<res<<endl;
    return 0;
}