poj 2891 Strange Way to Express Integers(扩展欧几里得，中国剩余定理)

题意：给出x%a1,x%a2……x%ak的余数r1,r2,r3……rk，求满足要求的最小的x的值。

思路：首先，满足第一个的最小解x1=r1,我们不妨先假设解为x=r1，对于第二个方程 x%a2=r2，为了让x满足第一个式子，那么x要满足x=x+a1*k，因此x+a1*k=a2*y+r2，然后用扩展欧几里得算法可以求得k的值，但此时k有可能为负，因此要取模……这样一直推下去，对于第i个式子，x要满足x=x+lcm(a1,a2…ai-1)，然后不断用扩展欧几里得求解即可，如果有一个式子无解，那么就没有答案。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}
ll ext_gcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    ll d=ext_gcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
ll lcm(ll x,ll y)
{
    return x/gcd(x,y)*y;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ll a,r,x,la,k,y,A,B,C,d;
    int n;
    while(cin>>n)
    {
        cin>>a>>r;
        x=r;
        la=a;
        bool flag=true;
        for(int i=1;i<n;++i)
        {
            cin>>a>>r;
            if(!flag) continue;
            A=la;B=a;C=r-x;
            d=ext_gcd(A,B,k,y);
            if(C%d!=0) {flag=false; continue;}
            C/=d;B/=d;
            k*=C;
            k=(k%B+B)%B;
            x=x+la*k;
            la=lcm(la,a);
        }
        if(!flag) x=-1;
        cout<<x<<endl;
    }
    return 0;
}