HDU3579 Hello Kiki 解线性同余方程的应用

解析不写了，原理还是看 poj2891的题解 ：http://blog.csdn.net/u010682557/article/details/14146251

很典型的看出 x≡Ai(mod Mi)求解，解线性同余方程的应用，

写完一直WA，后来仔细一看，题目要求输出的是最小的正整数解，所以当同余方程组的解为0的时候，就要输出 Mi 的最小公倍数

估计没少有人被坑吧

#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<cmath>
#include<memory.h>
#include<set>

#define ll long long
#define LL __int64
#define eps 1e-8
#define e 2.718281828

//const ll INF=9999999999999;

#define M 400000100

#define inf 0xfffffff

using namespace std;

//vector<pair<int,int> > G;
//typedef pair<int,int> P;
//vector<pair<int,int>> ::iterator iter;
//
//map<ll,int>mp;
//map<ll,int>::iterator p;
//
//vector<int>G[30012];

ll aa[12],rr[12];

ll GCD(ll a,ll b)
{
	if(b==0)
		return a;
	return GCD(b,a%b);
}

ll exgcd(ll a,ll &x,ll b,ll &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	ll r=exgcd(b,x,a%b,y);
	ll tmp=x;
	x=y;
	y=tmp-a/b*y;
	return r;
}

int main(void)
{
	int t;
	int Case=0;
	ll n;
	cin>>t;
	while(t--)
	{
		cin>>n;
		ll lcm=1;
		for(int i=0;i<n;i++)
		{
			cin>>aa[i];
			lcm=lcm*aa[i]/GCD(aa[i],lcm);
		}
		for(int i=0;i<n;i++)
			cin>>rr[i];
		ll x0,y0;
		bool flag=false;
		for(int i=1;i<n;i++)
		{
			ll a=aa[0],b=aa[i],c=rr[i]-rr[0];
			ll gcd=exgcd(a,x0,b,y0);
			if(c%gcd!=0)
			{
				flag=true;
				break;
			}
			ll MOD=b/gcd;
			x0=(x0*c/gcd%MOD+MOD)%MOD;
			rr[0]=aa[0]*x0+rr[0];
			aa[0]*=(aa[i]/gcd);
		}
		printf("Case %d: ",++Case);
		if(flag)
			puts("-1");
		else
		{
			if(rr[0]!=0)
				cout<<rr[0]<<endl;
			else
				cout<<lcm<<endl;
		}
	}
}