对于拓展欧几里德算法的理解

对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，必然存在整

数对 x，y ，使得 gcd（a，b）=ax+by。

c++语言的实现：

int exgcd(ll a,ll b,ll &x,ll &y)
{
    if(a==0)
    {
        x=0;y=1;
        return b;
    }
    else
    {
        ll tx,ty;
        ll d=exgcd(b%a,a,tx,ty);
        x=ty-(b/a)*tx;
        y=tx;
        return d;
    }
}

求解 x，y的方法的理解

设 a>b。

1，显然当 b=0，gcd（a，b）=a。此时 x=1，y=0；

2，ab<>0 时

设 ax ₁+ by ₁= gcd(a,b);

bx ₂+ (a mod b)y ₂= gcd(b,a mod b);

根据朴素的 欧几里德原理有 gcd(a,b) = gcd(b,a mod b);

则:ax ₁+ by ₁= bx ₂+ (a mod b)y ₂;

即:ax ₁+ by ₁= bx ₂+ (a - [a / b] * b)y ₂=ay ₂+ bx ₂- [a / b] * by ₂;

也就是ax ₁+ by1 == ay ₂+ b(x ₂- [a / b] *y ₂);

根据恒等定理得：x ₁=y ₂; y ₁=x ₂- [a / b] *y ₂;

这样我们就得到了求解 x ₁,y ₁ 的方法：x ₁，y ₁ 的值基于 x ₂，y ₂.

上面的思想是以递归定义的，因为 gcd 不断的递归求解一定会有个时候 b=0，所以递归可以结束。

对于上述说明，可以给道题目帮助理解：
http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=10755

先说一下大概题意：有两只青蛙，一只在坐标x，另一直在坐标y，青蛙x一次跳跃可以前进m单位距离，青蛙y一次跳跃可以前进n单位的距离，两青蛙都在同一纬度，该纬度长度为L。两只青蛙同方向同时跳啊跳，问你最少跳多少次，它们才可以相遇，如果不能相遇，输出impossble

 

分析一下：假设跳了T次以后，青蛙1的坐标是x+m*T，青蛙2的坐标是y+n*T.他们能相遇的情况是（x+m*T）-(y+n*T)==P*L，其中P为某一个整数，变形一下

得到(n-m)*T+P*L==x-y   我们设a=(n-m),b=L,c=x-y,T=x,P=y.于是便得到ax+by==c。激动啊，这不就是上面一样的式子吗~

直接套用扩展欧几里得函数，得到一组解x,y。由于问题是问最少跳多少次，于是只有x是我们需要的信息。那么再想，x是最小的吗？

答案是可能不是！那么如何得到最小解呢？  我们考虑x的所有解的式子： x=x0+b/d*t。x0是我们刚刚求到的，很显然右边是有个单调函数，当t为某一个与x正负性质相反的数时，可以得到最小的x。 令x的正负性质为正，那么x=x0-b/d*t1 (t1==-t)。令x==0，那么t=x0*d/b，最小的x等于x0减去t*b/d。这里得到的x可能是负数，如果是负数，我们再为它加上一个b/d即是所求答案了！

类似的题目还有poj2142,poj2891.