POJ2891 Strange Way to Express Integers 模线性方程组

　　题目链接：http://poj.org/problem?id=2891

　　题意：求解模线性方程组，且任意的ai可能不互质。

　　直接用扩展欧几里得求解就可了，具体可以看：http://www.cnblogs.com/zhsl/archive/2013/04/12/3017109.htm

//STATUS:C++_AC_0MS_196KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<map>
using namespace std;
#define LL __int64
#define pii pair<int,int>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int N=110,INF=0x3f3f3f3f,MOD=10000,STA=8000010;
const double DNF=1e13;

LL a[N],m[N];
int n;

void exgcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b){d=a;x=1;y=0;}
    else {exgcd(b,a%b,d,y,x);y-=x*(a/b);}
}

LL Modline(int n)
{
    LL d,x,y,A,M,Mod;
    A=a[n-1],M=m[n-1];
    n--;
    // m1*x-m2*y=a2-a1
    while(n--){
        exgcd(M,m[n],d,x,y);
        if((A-a[n])%d!=0){
            return -1;
        }
        Mod=m[n]/d;
        x=(x*((a[n]-A)/d)%Mod+Mod)%Mod;
        A+=M*x;
        M=M/d*m[n];
    }
    return A;
}

int main()
{
 //   freopen("in.txt","r",stdin);
    int i,j;
    while(~scanf("%d",&n))
    {
        for(i=0;i<n;i++)
            scanf("%I64d%64d",&m[i],&a[i]);

        printf("%I64d\n",Modline(n));
    }
    return 0;
}