标签：type string ace quick void def 中国 -o mes
现有两组数字，每组k个，第一组中的数字分别为：a1，a2，...，ak表示，第二组中的数字分别用b1，b2，...，bk表示。其中第二组中的数字是两两互素的。求最小的非负整数n，满足对于任意的i，n - ai能被bi整除。
输入数据的第一行是一个整数k，（1 ≤ k ≤ 10）。接下来有两行，第一行是：a1，a2，...，ak，第二行是b1，b2，...，bk
输出格式：输出所求的整数n。
所有数据中，第一组数字的绝对值不超过10^9（可能为负数），第二组数字均为不超过6000的正整数，且第二组里所有数的乘积不超过10^18
每个测试点时限1秒
注意：对于C/C++语言，对64位整型数应声明为long long，如使用scanf, printf函数（以及fscanf, fprintf等），应采用%lld标识符。
#include<iostream> #include<cstring> #include<string> #include<cmath> #include<stdio.h> using namespace std; typedef long long ll; ll a[20],b[20],n; void exgcd(ll a,ll b,ll &x,ll &y) { if(b==0) x=1,y=0; else{ exgcd(b,a%b,y,x); y-=a/b*x; } } ll quickmul(ll a,ll b,ll mod) { ll ans=0; while(b){ if(b&1)ans=(ans+a)%mod; a=(a+a)%mod; b>>=1; } return ans; } ll solve() { ll M=1; for(int i=1;i<=n;i++){ M*=b[i]; } ll ans=0; for(int i=1;i<=n;i++){ ll x,y; exgcd(M/b[i],b[i],x,y); //M/b[i]*x==1(mod b[i])的解x x=(x%b[i]+b[i])%b[i];//化为最小正整数解 ans=(ans+quickmul(quickmul(M/b[i],x,M),a[i],M))%M;//运用快速乘，不然会爆 } return ans; } int main() { scanf("%lld",&n); for(int i=1;i<=n;i++)scanf("%lld",&a[i]); for(int i=1;i<=n;i++)scanf("%lld",&b[i]); for(int i=1;i<=n;i++)a[i]=(a[i]%b[i]+b[i])%b[i]; printf("%lld\n",solve()); return 0; }
标签：type string ace quick void def 中国 -o mes
原文地址：https://www.cnblogs.com/cherish-lin/p/11134599.html
