题目
题目连接: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;
}