扩展欧几里得（最小正整数解）

 

题意：

两只青蛙在一个循环坐标轴上，朝一个方向同时开始跳，问何时两只青蛙能落在同一个点上。

列出方程：(x + m * t) % L = (y + n* t) % L

根据同余定义：如果 a-b能被n整除，则a≡b(mod n), 即a%n==b%n，a和b同余n。

则上式可写成：

（x+m*t）-(y+n*t)==k*L

(x-y)+(m-n)*t==k*L

(m-n)*t - L*k== y-x

到这一步可以比较明显的看出来是一个 ax+by=k*(gcd(a,b)) ，可以用扩展欧几里得算法来求得x,y。

必须要清楚一点：exgcd解出来的那么一组x,y并不是我们所求的最小正整数解，它甚至可能是一个负的不知道哪里去的一个数。它们仅仅是一组解。

#include<iostream>
using namespace std;
typedef long long ll;
ll x,y,m,n,l;

ll exgcd(ll a,ll b,ll &x1,ll &y1)
{
    if(!b)
    {
        x1=1;y1=0;
        return a;
    }
    ll x2,y2;
    ll d=exgcd(b,a%b,x2,y2);
     x1=y2;
     y1=x2-(a/b)*y2;
    return d;
}

ll inv(ll a,ll b,ll c)
{
    ll xx,yy;
    ll d=exgcd(a,b,xx,yy);
    if(c%d)return -1;
    xx=xx*c/d;//按照比例扩大
    ll t=b/d;
    if(t<0)t=-t;
    return ((xx%t)+t)%t;
}
int main()
{
    cin>>x>>y>>m>>n>>l;
    ll res=inv(m-n,l,y-x);
  if(res==-1)cout<<"Impossible"<<endl;
  else cout<<res<<endl;
    return 0;
}

        inv函数中对x结果做出了处理：为求得x的最小正整数解，先是设了t = b/d(即b/gcd)，gcd不保证为正，所以如果t为负将它取为正。最后对x做出处理，先把x对t取模得出最小的解，再将这个值加模取模，将它变为一个正数，即可求得最小正整数的解。

为何用b/gcd来取模？
已知扩展欧几里得中的通解公式：

(x1,y1)是一组通解。由此可见，对于这个方程，x的解的周期为b/gcd，y的解的周期为a/gcd。用任意解 x 对 b/gcd 取模，能够得到的是最小解。再对这个解取正，即是最小正整数解。
求y最小正整数解同理，t = a/gcd。