歐幾里得算法與擴展歐幾里得算法

注意：歐幾里得算法和擴展歐幾里得算法是解決不一樣問題的兩種，其中擴展歐幾里得算法須要用到歐幾里得算法

update in 2018.8.22：補充歐幾里得算法及其相關證實

歐幾里得算法

歐幾里得算法是用來求$(a, b)$的最大公約數

若$d \mid a$且$d \mid b$，則稱$d$爲$(a, b)$的公約數

在$(a, b)$的公約數中，最大的$d$成爲$(a, b)$的最大公約數

設$gcd(a, b)$表示$(a, b)$的最大公約數

性質1：$gcd(a, b) = gcd(a - b, b)$

證實：

不妨設$a > b$，$d = gcd(a, b)$

那麼$a = dx$，$b = dy$。其中$x > y$

那麼$a - b = d(x - y)$

即$gcd(a, b ) = gcd(dx, dy) = gcd(d(x - y), dy) = d = gcd(a - b, b)$

性質2：$gcd(a , b) = gcd(a \% b, b)$

證實：由性質1及模運算的性質不可貴到

這樣直接遞歸計算便可，每次遞歸都有一個數減少一半

所以複雜度爲$log(a + b)$

 

擴展歐幾里得算法

用途

當咱們已知$a,b$

擴展歐幾里得算法能夠求出知足$a*x+b*y=gcd(a,b)$的$(x,y)$解集

$gcd(a,b)$表示$a,b$的最大公約數

 

前導知識

上面有證實。

$gcd(a,b)=gcd(b,a\%b)$

$gcd(a,0)=0$

$a\%b=a-a/b*b$

 

推導過程

其實擴展歐幾里得的推導過程挺天然的

$a*x+b*y$

$=gcd(a,b)$

$=gcd(b,a\%b)$

$=b*x+(a\%b)*y$

$=b*x+(a-a/b*b)*y$

$=b*x+a*y-a/b*b*y$

$=a*y+b*x-a/b*b*y$

$=a*y+(x-y*a/b)*b$

這樣不斷的遞歸下去

當$b=0$時

$x=1,y=0$

 

代碼

注意：

咱們在求$(x-y*a/b)$的時候須要用到上一層的$x$

但此時上一層$x$已經被賦值成了$y$

因此咱們須要開一箇中間變量來記錄上一層的$x$

 

int exgcd(int a, int b, int &x, int &y) {
    if(!b) {x = 1; y = 0; return a;}
    int r = exgcd(b, a % b, x, y);
    int tmp = x; x = y, y = tmp - a / b * y;
    return r;
}

 

應用

1

擴展歐幾里得最重要的應用就是求形如$a*x+b*y=c$的解

那麼如何求呢？

首先，這個方程可以能力的條件是$c\%gcd(a,b)=0$，這個應該比較顯然

根據前面將的擴展歐幾里得算法

咱們能夠先求出$a*x_0+b*y_0=gcd(a,b)$的解$x_0,y_0$

而後方程兩邊同時除以$gcd(a,b)$

就獲得$a*x_0/gcd(a,b)+b*y_0/gcd(a,b)=1$的解

再在方程兩邊同乘$c$

就獲得了方程

$a*x_0/gcd(a,b)*c+b*y_0/gcd(a,b)*c=c$

是否是很簡單？

2

若$gcd(a,b)=1$，且$x0,y0$爲$a*x+b*y=c$的一組解，則該方程的任一一解能夠表示爲

$x=x_0+b*t,y=y_0-a*t$

證實:

$a*x+b*y$

$=a*(x_0+b*t)+b*(y_0-a*t)$

$=a*x_0+a*b*t+b*y_0-a*b*t$

$=a*x_0+b*y_0$

 

例題

洛谷P1516 青蛙的約會

根據題目要求列出等式，化簡便可

題解

洛谷P2421 [NOI2002]荒島野人

題解