poj 1061 青蛙的约会(扩展欧几里得解不定方程)

设两只青蛙跳了t步，A的坐标为x+mt,B的坐标为y+nt，他们相遇的时候满足x+mt-(y+nt) = pL(p表示两青蛙走过的路程相差p圈)
移项后：(n-m)t+Lp=x-y
设n-m为A，x-y为B，则At+Lp=B，，求满足这个式子的最小t，也就是求一次同余方程At=B(mod L)的最小整数解。求出GCD(A,L),如果B%GCD(A,L)!=0,则无解。解出的t也可能是负数，要转化为正数，我一直都错在这里。

从别人博客摘的：
欧几里得算法的拓展应用中有如下三条定理：
定理一：如果d = gcd(a, b)，则必能找到正的或负的整数k和l，使d = a*x+ b*y。
定理二：若gcd(a, b) = 1，则方程ax ≡ c (mod b)在[0, b-1]上有唯一解。
定理三：若gcd(a, b) = d，则方程ax ≡ c (mod b)在[0, b/d - 1]上有唯一解。
http://www.cnblogs.com/nwpuacmteams/articles/5545976.html

因为解的范围在[0, b/d - 1]，所以(t%(b/d)+b/d)%(b/d)，就可以得到最小的正整数解。因为在b/d之内，所以要模它。

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;

long long extendGcd(long long a, long long b, long long &x, long long &y)
{
    if(a == 0 && b == 0) return-1;
    if(b == 0)
    {
        y = 0;
        x = 1;
        return a;
    }
    long long d = extendGcd(b,a%b,y,x);
    y -= a/b*x;
    return d;
}

int main()
{
    long long x,y,m,n,l;
    while(cin >> x >> y >> m >> n >> l)
    {
        if(n == m)
        {
            printf("Impossible\n");
            continue;
        }
        long long a = n-m;
        long long c = x-y;
        long long b = l;
        long long rx = 0, ry = 0;
        long long d = extendGcd(a,b,rx,ry);
        if(c%d != 0)
        {
            printf("Impossible\n");
            continue;
        }
        int r = b/d;
        printf("%d\n",(((rx*c/d)%r)+r)%r);
    }
    return 0;
}