poj-1061青蛙的约会（扩展欧几里德算法）

欧几里德算法（辗转相除法，求两个数的最小公约数）
扩展欧几里德算法是用来在已知a, b求解一组x，y，使它们满足贝祖等式等式： ax+by = gcd(a, b) =d（解一定存在，根据数论中的相关定理）。扩展欧几里德常用在求解模线性方程及方程组中。
由欧几里德算法知 ：gcd（a，b）= gcd（b，a%b）；前提b不为0；
第一次看到这个题用暴力去解(x+mk)modL=(y+nk)modL;(k为跳的次数)，结果wrong；
从网上看解题报告才知道是用扩展欧几里德算法；
将上式转变得：l*p=(x+mk)-(y-nk) （p为青蛙A,B的相差圈数）——>k*(n-m)+l*p=x-y;
此时是不是特像ax+by=c;
所以令 a=（n-m）,b=l, c=x-y;
因为已知ax+by=gcd(a,b)有解可求；（x是k，y是p）；gcd(b,a%b)=bx1+a%by1;
由欧几里德算法知 gcd(a,b)=gcd(b,a%b);
所以：ax+by=bx1+a%by1; a%b=a-a/b*b(这里是程序算法中的乘除法)
代入化简 得：ax+by=ay1+b(x1-a/b*y1);
所以 x=x1; y=x1-a/b*y1;
所以代码int exGcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;y=0;
return a;
}
int r=exGcd(b,a%b,x,y);
int t=x;x=y;y=t-a/b*y;
return r;
}
exgcd中得出最小公约数，而且有一组x,y;但不是题目中那个，题目是求满足ax+by=c的最小正整数x；
这里求得了ax+by=gcd的得数；
只有C整除gcd时才有解
所以两边乘以c/gcd; 即x=x*(c/gcd);
当x>=0时，x=xmod（b/gcd）;(因为求最小的正整数)
当x<0时，x=x%（b/gcd）+(b/gcd);

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

#include<iostream>
#include<stdio.h>
#define ll long long int
using namespace std;
long long exgcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0){
        x=1;
        y=0;
        return a;
    }
    else{
        int r=exgcd(b,a%b,y,x);
        y-=a/b*x;
        return r;
    }
}
int main()
{
    ll x,y,m,n,l;
    ll a,b,c,t,r,gcd;
    while(scanf("%lld %lld %lld %lld %lld",&x,&y,&m,&n,&l)!=EOF){
        a=n-m;
        b=l;
        c=x-y;
        if(a<0){
            a=-a;
            c=-c;
        }
        gcd=exgcd(a,b,x,y);
        if(c%gcd!=0) printf("Impossible\n");
        else{
            t=c/gcd;
            x=x*t;
            r=b/gcd;
            if(x>=0)
            {
                x=x%r;
            }
            else{
                x=x%r+r;
            }
            printf("%lld\n",x);
        }
    }
    return 0;
}