题意:
两只青蛙在一个循环坐标轴上,朝一个方向同时开始跳,问何时两只青蛙能落在同一个点上。
列出方程:(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。