poj2891 Strange Way to Express Integers【解一元线性同余方程组模板】

题目大意：

求解一元线性同余方程组:
$x\equiv b_1(mod$ $m_1)$
$x\equiv b_2(mod$ $m_2)$
……
$x\equiv b_n(mod$ $m_n)$
若无解输出-1，否则输出最小正整数解。

解题思路：

考虑只有两个方程：
$x\equiv b_1(mod$ $m_1)$
$x\equiv b_2(mod$ $m_2)$
设$(m_1,m_2)=d$

可以转化为：
$x=m_1y_1+b_1$
$x=m_2y_2+b_2$

联立移项得
$m_2y_2+m_1y_1(y_1正负无所谓)=b_1-b_2$

这时我们就可以用拓展欧几里得把最小的$y_2$求出，而我们知道$y_2+\frac{m_1}{d}t,t\in Z$都是该方程的解，所以得到：
$x=m_2(y_2+\frac{m_1}{d}t)+b_2=m_2y_2+b_2+\frac{m_1m_2}{d}t$。

我们取特解$x_1=m_2y_2+b_2$，所以$x=x_1+\frac{m_1m_2}{d}t$

这时假设有第三个方程：
$x\equiv b_3(mod$ $m_3)$
那么同理：$x=m_3y_3+b_3=x_1+\frac{m_1m_2}{d}t$
移项得：$m_3y_3+\frac{m_1m_2}{d}t=x_1-b_3$，
可以发现这与$m_2y_2+m_1y_1=b_1-b_2$同构，所以迭代求解即可。

由$x=x_1+\frac{m_1m_2}{d}t$可知最小正整数解在$[1,\frac{m_1m_2}{d}]$之间，取模即可。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<ctime>
#include<vector>
#include<set>
#define ll long long
using namespace std;

int n;
ll b1,b2,m1,m2,d;

ll gcd(ll a,ll b)
{
    return b?gcd(b,a%b):a;
}

void exgcd(ll a,ll b,ll &x,ll &y)
{
    b?(exgcd(b,a%b,y,x),y-=a/b*x):(x=1,y=0);
}

ll solve()
{
    ll x,y,a=m2,b=m1,c=b1-b2;
    d=gcd(a,b);
    if(c%d!=0)return -1;
    a/=d,b/=d,c/=d;
    exgcd(a,b,x,y);
    x*=c;x=(x%b+b)%b;
    return x;
}

int main()
{
    //freopen("lx.in","r",stdin);
    //freopen("lx.out","w",stdout);
    ll y2;
    while(scanf("%d",&n)!=EOF)
    {
        y2=0;
        scanf("%lld%lld",&m1,&b1);
        for(int i=2;i<=n;i++)
        {
            scanf("%lld%lld",&m2,&b2);
            if(y2==-1)continue;
            y2=solve();
            m1=m1*(m2/d);b1=(b2+m2*y2)%m1;
        }
        if(y2==-1)puts("-1");
        else cout<<b1<<'\n';
    }
    return 0;
}