P3868 [TJOI2009]猜数字

题目描述
现有两组数字，每组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
说明/提示
所有数据中，第一组数字的绝对值不超过10^{9（可能为负数），第二组数字均为不超过6000的正整数，且第二组里所有数的乘积不超过10}18

每个测试点时限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;
}