Strange Way to Express Integers （模线性方程组（扩展欧几里得））

题目来源：https://vjudge.net/problem/POJ-2891
【题意】
给出两个a，r数组，求出一个尽可能小的数s保证对于每一组a[i],r[i]，都符合s%a==r。
【思路】
模线性方程组模板题。
列式解释：
S%a[1]==r[1]. 1
S%a[2]==r[2]. 2
式子1中S的通解是S=r[1]+k*a[1]（k属于0,1,2…）
式子2中S的通解是S=r[2]+k*a[2]（k属于0,1,2…）
因为S是一样的。所以可以把两个式子合成一个：
T%lcm(a[1],a[2])==S （此处的S是符合上面两式的最小的解）
为什么是lcm（最小公倍数）呢？
因为T的通解是T=S+k*lcm（a[1],a[2]）;
【代码】

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<limits.h>
#include<algorithm>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int mod=1e9+7;
typedef unsigned long long ll;
typedef long long LL;
LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    LL d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
int main()
{
    LL a1,r1,a2,r2,x,y,k;
    while(~scanf("%lld",&k))
    {
        bool flag=1;
        scanf("%lld%lld",&a1,&r1);
        for(int i=1; i<k; i++)
        {
            scanf("%lld%lld",&a2,&r2);
            LL a=a1,b=a2,c=r2-r1;
            LL d=exgcd(a,b,x,y);
            if(c%d!=0)
            {
                flag=0;
            }
            else
            {
                x=((c/d)*x%(b/d)+(b/d))%(b/d);//求出一个特解，有一个解系：X=x+lcm(a1,a2);
                r1=x*a1+r1;//计算答案，必须与上一行代码相对应
                a1=a/d*b;//最小公倍数
            }
        }
        if(!flag)
            printf("-1\n");
        else
            printf("%lld\n",r1);
    }
}

参考：http://blog.csdn.net/qq_21120027/article/details/51340218