pku 1061 欧几里德扩展定理

       早在几个月之前就在做北大上"青蛙的约会"那道题，当时想得简单，以
为这个题目像两个人在追赶一样，只要一个能追上另一个人就达到目的了。其实
不然，青蛙是跳着走的，一只青蛙超过了另一只青蛙，它们不一定能见着面，因
为只有当它们跳到同一点时，才能算它们见面。
       设两只青蛙分别为A,B，它们所在的坐标分别为x,y,每一步能跳m,n米,纬度总长为L米,如果x = y；那么它们不需要走就能见面，所走步数为0。否则
当m>n时，有(m-n)*p = -(x-y)+q*L； 同样，当m < n时，有(n-m)*p=-(y-x)
+q*L；其中q = 1,2,3....我们要求的是满足条件的第一个p的值。
把等式变形我们可以发现，它就可以转化为一个欧几时德的扩展公式，即对于
一个等式ax + by = c；如果c mod gcd(a,b)  == 0，那么它有无数组整解，其求解过程如下：
由扩展欧几里德定理可知：对于两个都不为0的非负整数a,b,令gcd(a,b)为a,b
的最大公约数,那么必存在整数x,y使得ax+by=gcd(a,b)成立。
其算法是这样的：

1. int exGcd(int a,int b,int &x,int &y)
2. {
3.         if(b == 0)
4.         {
5.                 x = 1;
6.                 y = 0;
7.                 return a;
8.          }
9.          int r = exGcd(b,a%b,x,y),t;
10.          t = x;
11.          x = y;
12.          y = t-a/b*y;
13.          return r;
14. }

其中x,y为满足条件的第一组解，函数返回的值是gcd(a,b)；
那么方程ax+by=c的第一组解为X=x*(c/gcd(a,b)), Y = y*(c/gcd(a,b))。
其它的整数解满足：p = X+b/gcd(a,b)*k；
                 q = Y-a/gcd(a,b)*k；（其中k为任意整数）