首页 > 代码库 > 欧几里得算法以及扩展欧几里得算法(过河noip2005提高组第二题)

欧几里得算法以及扩展欧几里得算法(过河noip2005提高组第二题)

欧几里得算法:也被称作辗转相除法

gcd(a,b)=gcd(b,a%b);

终止条件a=gcd b=0;

(gcd为a,b的最大公约数)

 

扩展欧几里得算法: a 和 b 的最大公约数是 gcd ,一定能够找到这样的 x 和 y ,使得: a*x + b*y = gcd 成立

我们只需要找到特殊解x0,y0;

则通解为 x = x0 + (b/gcd)*t    y = y0 – (a/gcd)*t

那如何求出下一组解呢

仿照欧几里得算法a=b,b=a%b代入。

a%b = a - (a/b)*b(这里的 “/” 指的是整除)可以得到:

        gcd = b*x1 + (a-(a/b)*b)*y1

            = b*x1 + a*y1 – (a/b)*b*y1

            = a*y1 + b*(x1 – a/b*y1)

对比易得      x = y1

                  y = x1 – a/b*y1

终止为    a*1 + b*0 = gcd即a的系数为1,b的系数为0或者其他值(*0=0)

我们再看过河这个问题 设a=s,b=s+1,ax+by=n a,b的最大公约数为1因为是相邻的自然数

ax+by=n是由ax+by=1(1为最大公约数)de系数*n得到的。

未完待续。

 

欧几里得算法以及扩展欧几里得算法(过河noip2005提高组第二题)