题目描述
现有两组数字,每组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标识符。
∵ n-a[i]%b[i]==0
∴n同余a[i](mod b[i]);
所以这是一个求解线性方程组的,所以用中国剩余定理去写。
但题目的数据非常坑,存在负数,我们要把负数变成正数才可以。还有注意的一点就是,数据很有可能爆long long ,所以我们要用快速乘去做。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=110;
int n;
ll a[N],b[N];
ll sum,x,y;
ll exgcd(ll a,ll b,ll &x,ll &y) //扩展欧几里德
{
if(b==0) {x=1;y=0;return a;}
ll tmp=exgcd(b,a%b,x,y);
ll z=x;
x=y;
y=z-a/b*y;
return tmp;
}
ll ksc(ll a, ll b,ll p) //快速乘
{
ll res=0;
while(b)
{
if(b&1) res=(res+a)%p;
a=(a<<1)%p;b>>=1;
}
return res;
}
ll china()
{
sum=1;
ll ans=0;
for(int i=1;i<=n;i++) sum*=b[i];
for(int i=1;i<=n;i++)
{
ll w=sum/b[i];
exgcd(b[i],w,x,y);
y=(y%b[i]+b[i])%b[i];
ans=(ans+ksc(ksc(y,w,sum),a[i],sum))%sum;
}
return (ans+sum)%sum;
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++) cin>>a[i];
for(int i=1;i<=n;i++) cin>>b[i];
for(int i=1;i<=n;i++) a[i]=(a[i]%b[i]+b[i])%b[i]; //a[i]可能存在负值
cout<<china()<<endl;
}