题目描述
现有两组数字,每组k个,第一组中的数字分别为:a1,a2,...,ak表示,第二组中的数字分别用b1,b2,...,bk表示。其中第二组中的数字是两两互素的。求最小的非负整数n,满足对于任意的i,n - ai能被bi整除。
输入输出格式
输入格式:
输入数据的第一行是一个整数k,(1 ≤ k ≤ 10)。接下来有两行,第一行是:a1,a2,...,ak,第二行是b1,b2,...,bk
输出格式:
输出所求的整数n。
输入输出样例
输入样例#1: 复制3 1 2 3 2 3 5输出样例#1: 复制
23
说明
所有数据中,第一组数字的绝对值不超过109(可能为负数),第二组数字均为不超过6000的正整数,且第二组里所有数的乘积不超过1018
每个测试点时限1秒
注意:对于C/C++语言,对64位整型数应声明为long long,如使用scanf, printf函数(以及fscanf, fprintf等),应采用%lld标识符。
1 #include"bits/stdc++.h" 2 using namespace std; 3 typedef long long ll; 4 5 ll k,a[11],b[11],s=1,ans=0; 6 7 8 void exgcd(ll a,ll b,ll &x, ll &y) 9 { 10 if (!b){ 11 x=1,y=0 ;return ; 12 } 13 14 exgcd(b,a%b,x,y); ll t; 15 t=x; 16 x=y; 17 y=t-(a/b)*y; 18 } 19 20 inline ll mul(ll a,ll b) 21 { 22 ll r=0; 23 while (b) 24 { 25 if (b&1)r+=a%=s; 26 b>>=1; 27 a+=a%=s; 28 }return r; 29 } 30 int main() 31 { 32 cin>>k; 33 for (int i=1;i<=k;i++)cin>>a[i]; 34 for (int i=1;i<=k;i++)cin>>b[i],s*=b[i]; 35 for(int i=1;i<=k;i++) 36 { 37 ll x,y; 38 exgcd(s/b[i],b[i],x,y); 39 x=(x%b[i]+b[i])%b[i];(1) 40 // cout<<i<<" "<<x<<endl; 41 ans+=mul(s/b[i]*x,((a[i]%b[i]+b[i])%b[i]));//快速乘,否则爆long long(2)
// 一和二式都取模意义下的最小正整数
42 ans%=s; 43 } 44 cout<<ans; 45 }