POJ1061 青蛙的约会 推广的欧几里得算法

Problem Address：http://poj.org/problem?id=1061

 

详细的推广欧几里得算法已经在http://blog.csdn.net/Human_CK/archive/2011/05/23/6439658.aspx写过了。

 

觉得找几道推广的欧几里得来做，最后找到了这道题。

 

很纯粹的题目。

 

很容易理解。

 

但是收获了点东西。

 

【代码】

 

#include <iostream> using namespace std; __int64 exgcd(__int64 a, __int64 b, __int64 &x, __int64 &y) { __int64 t, r; if (b==0) { x = 1; y = 0; r = a; } else { r = exgcd(b, a%b, x, y); t = x; x = y; y = t - a/b*y; } return r; } int main() { __int64 x,y,m,n,l; __int64 a, b, c, d, t, xx, yy; while(scanf("%I64d %I64d %I64d %I64d %I64d", &x, &y, &m, &n, &l)!=EOF) { a = ((m-n)%l+l)%l; b = l; c = ((y-x)%l+l)%l; d = exgcd(a, b, xx, yy); if (c%d!=0) { printf("Impossible/n"); } else { t = b/d; xx *= c/d; printf("%I64d/n", (xx%t+t)%t); } } return 0; }

 

【注意】

 

一般情况下，最后我们都是求出了x（这里用xx表示）。

 

但是当x为负时，我们必须把它化为正数。

 

那么怎么化为正数呢？

 

一个不定式 3x+4y=1。解得x=-1，y=1。

 

所以，所有解应该是：x=-1+4t，y=1-3t。这里的-1和1为上面的解。而4和3分别为y和x的系数。

 

令t=0，得到x=-1，y=1。

 

令t=1，得到x=3，y=-2。

 

把x=-1+4t，y=1-3t带进原式3x+4y=1可验证其是成立的。

 

而这里的4和3其实是4和3分别除以它们的最大公约数，即代码中 xx*=c/d 中的 c/d。

 

这样子，结果就出来了。

 

【P.S】

 

这样子推广的欧几里得的这部分简单应用算是学好了。

 

不过还有一个逆元的应用，接下来还要好好学习。