poj 1061 青蛙的约会（拓展欧几里得学习）

题意

青蛙A每天走$m$，青蛙B每天走$n$,它们在一个环形的轨道上，轨道长度为$l$，青蛙A的初始位置为$a$,青蛙B的初始位置为$b$.
求：青蛙走多少天才可以相遇？

题解

设青蛙走$x$天可以相遇，则可以列出方程
$mx+a-(nx+b)=kl$
整理可得：
$(n-m)x+kl=a-b$可以看出这是一个$ax+by=c$的式子，可以用拓展欧几里得，求出一个$(n-m)x+kl=gcd(n-m,l)$的解$x，y$
由于不确定$n,m$的大小，需要分类讨论一下
$(n-m)x+kl=a-b(n>m)$
$(m-n)x+(-k)l=b-a(n<m)$
因为我们不需要求$k$的大小，只需要求 x，所以我们只处理系数变换。

最后我们需要求出来的应该是一个最小正整数$x_{min}$

经过拓展欧几里得，我们可以求出：
$(n-m)x+kl=gcd(n-m,l)$的一个特解$(x,y)$
$x=x\frac{a-b}{gcd(n-m,l)}$
得到$(n-m)x+kl=a-b$的一个特解
这个方程式的通解为$X=x_{min}+k\frac{b}{gcd(n-m,l)}(k\in Z)$
那么$x_{min}=X\%(\frac{b}{gcd(n-m,l)})$;
因为有求出来$x<0$的情况
$x_{min}=[x\%(\frac{b}{gcd(n-m,l)})+\frac{b}{gcd(n-m,l)}]\%(\frac{b}{gcd(n-m,l)})$

代码

#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
ll m,n,a,b,l;
ll gcd;
void ex_gcd(ll a,ll b,ll &x,ll &y){
	if(b){
	    ex_gcd(b,a%b,x,y);
	    ll t=x;
	    x=y;
	    y=t-(a/b)*y;
	}
	else{
	    x=1,y=0;
	    gcd=a;
	}
	return;
}
int main()
{
	cin>>a>>b>>m>>n>>l;
	ll x=0,y=0;
    	ll c=a-b;
   	 if(n-m>0)
	     ex_gcd(n-m,l,x,y);
   	 else{
   	     c=-c;
   	     ex_gcd(m-n,l,x,y);
   	 }
   	 if(c%gcd!=0){
             cout<<"Impossible"<<endl;
             return 0;
         }
         ll t=l/gcd;
   	 x*=(c/gcd);
   	 x=(x%t+t)%t;
   	 cout<<x<<endl;
   	 
}