青蛙的约会

青蛙的约会

两只青蛙在网上相识了，它们聊得很开心，于是觉得很有必要见一面。
它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝西跳，直到碰面为止。
可是它们出发之前忘记了一件很重要的事情，既没有问清楚对方的特征，也没有约定见面的具体位置。
不过青蛙们都是很乐观的，它们觉得只要一直朝着某个方向跳下去，总能碰到对方的。但是除非这两只青蛙在同一时间跳到同一点上，不然是永远都不可能碰面的。
为了帮助这两只乐观的青蛙，你被要求写一个程序来判断这两只青蛙是否能够碰面，会在什么时候碰面。
我们把这两只青蛙分别叫做青蛙 A和青蛙 B，并且规定纬度线上东经 0 度处为原点，由东往西为正方向，单位长度 1米，这样我们就得到了一条首尾相接的数轴。设青蛙 A的出发点坐标是 x，青蛙 B 的出发点坐标是 y。
青蛙 A一次能跳 m 米，青蛙 B 一次能跳 n米，两只青蛙跳一次所花费的时间相同。纬度线总长 L米。
现在要你求出它们跳了几次以后才会碰面。
输入格式
输入只包括一行 5个整数 x，y，m，n，L。
输出格式
输出碰面所需要的跳跃次数，如果永远不可能碰面则输出一行 Impossible。
数据范围
x≠y<=2000000000,0<m,n<=2000000000,0<L<=2100000000
输入样例：
1 2 3 4 5
输出样例：
4

解题思路

扩展欧几里得算法

形如方程： ax + by =GCD(a，b)
根据数论中的相关定理，解一定存在。

LL exgcd(LL a,LL b,LL &x, LL &y)//扩展欧几里得算法  形如ax+by=GCD(a，b) 
{                                    //x,y返回一组解 
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    LL ans=exgcd(b,a%b,y,x);
    y-= a/b * x;
    return ans;//求出系数a，b的最大公约数 GCD(a,b)=GCD(b,a%b)
}

线性同余方程

形如方程： ax + by = c 有整数解的充要条件是 GCD(a,b)|c
线性同余方程的一个最小整数解x = (x%t+t)%t 其中t=b/GCD(a,b)

问题分析

两个青蛙跳到同一个点上才算是遇到了，所以得到如下式子：
（x+m * t） - （y+n * t） =p * L
t：跳的次数 , p：两青蛙跳的圈数差
路程差等于纬度线L的整数倍
变换方程： （n-m）* t +L * p =x-y
令： a=n-m ， b=L, c = GCD(a,b), d =x-y
有 a * t + b * p =d 此方程为系数为a，b的线性同余方程。
我们要求t 的最小正整数解,t不会是零。

求解ｔ的步骤

先用扩展欧几里得算法求出其中一组解 t0，p0，满足方程
a * t0 + b * p0 = GCD(a，b） ，c = GCD(a,b)
需要d%c=0 ，即必须 GCD(a,b)|(x-y) ， 否则无解。

两边都乘以（d/c），得到： a * t0 *（d/c) + b * p0 (d/c) =d
所以： t0*d/c是同余方程的一个特解，可能是负数。
又有： a * （t0 * （d/c）+ b／c*N ）+b * （p0（d/c）-a／c * N）=d，其中N是任意整数
所以，最小正整数解为 (t0*（d/c）%(b/c) + (b/c))%(b/c)

代码方案

#include<iostream>
#include<cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
LL exgcd(LL a,LL b,LL &x, LL &y) 
{                             
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    LL ans=exgcd(b,a%b,x,y);
    LL tmp = x;
    x = y;
    y = tmp - a/b *y;
    return ans;
}

int main()
{
    LL n,m,x,y,l,X,Y,gcd;
    cin>>x>>y>>m>>n>>l;
    
    if((x-y)%(gcd=exgcd(n-m,l,X,Y))!=0)
      cout<<"Impossible"<<endl;
    else{
        X *=((x-y)/gcd);
        LL M = abs(l /gcd);
        cout<<(X%M+M)%M<<endl;
        }
    return 0;
}