题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=3579
这是一道基本的孙子定理的数论题目。
代码如下:
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
using namespace std;
/*
freopen("input.txt", "r", stdin); //读数据
freopen("output.txt", "w", stdout); //注释掉此句则输出到控制台
*/
int w[10],b[10],n;
//扩展Euclid求解gcd(a,b)=ax+by
int ext_gcd(int a,int b,int& x,int& y)
{
int t,ret;
if (!b)
{
x=1,y=0;
return a;
}
ret=ext_gcd(b,a%b,x,y);
t=x,x=y,y=t-a/b*y;
return ret;
}
int China_left2()
{
int flag,w1,w2,b1,b2,gcd,x,y,i,t;
flag=0;w1=w[0];b1=b[0];
for(i=1;i<n;i++)
{
w2=w[i];b2=b[i];
gcd=ext_gcd(w1,w2,x,y);
if((b2-b1)%gcd)
{
flag=1;break;
}
t=w2/gcd;
x=(x*(b2-b1))/gcd;
x=(x%t+t)%t;
b1=w1*x+b1;
w1=(w1*w2)/gcd;
b1=(b1%w1+w1)%w1;
}
if(flag==1)
b1=-1;
if(b1==0&&n>1)
b1=w1;
if(b1==0&&n==1)
b1=w[0];
return b1;
}
int main()
{
int i,T,casei=0;
cin>>T;
while(T--)
{
scanf("%d",&n);
for(i=0;i<n;i++)
scanf("%d",&w[i]);
for(i=0;i<n;i++)
scanf("%d",&b[i]);
printf("Case %d: %d\n",++casei,China_left2());
}
return 0;
}