Rozšírený Euklidov algoritmus

Šablóna:Na úpravu Rozšírený Euklidov algoritmus je algoritmus, ktorým je možné nájsť Bézoutovú rovnosť, čiže vyjadriť najväčší spoločný deliteľ (NSD) dvoch kladných celých čísel ich lineárnou kombináciou.

Príklad

Nech sú dané dve kladné celé čísla m a n a nech d je ich najväčší spoločný deliteľ, t.j. d = NSD(m,n). Pomocou Euklidovho algoritmu je možné zistiť, že NSD(196,34) = 2:

i	ci	di	pi	zi
0	196	34	5	26
1	34	26	1	8
2	26	8	3	2
3	8	2	4	0

kde $i$ je iterácia, $c_0=m=196$, $d_0=n=34$, $p_i$ je celočíselný podiel $c_i/d_i$, $z_i$ je zvyšok po delení $c_i/d_i$, čiže platí $c_i=d_i{p}_i+z_i$, resp:

$z_i=c_i-d_i{p}_i(1)$

a pre každé $i=1,2,3,...$ je

$\begin{array}{cccc}{c}_i& =& d_{i-1}& (2)\\ d_i& =& z_{i-1}& (3)\end{array}$

Úlohou Rozšíreného Euklidovho algoritmu je nájsť dvojicu celých čísel $a$ a $b$, spĺňajúcu rovnosť $2=196a+34b$.

Nájdenie Bézoutovej rovnosti spätným dosadzovaním

V uvedenom príklade platí:

$\begin{array}{cccc}2& =& 26-8\cdot 3& (4)\\ 8& =& 34-26\cdot 1& (5)\\ 26& =& 196-34\cdot 5& (6)\end{array}$

Dosadením ľavej strany rovnosti (5) do rovnosti (4) dostaneme:

$2=34\cdot (-3)+26\cdot 4(7)$

a dosadením ľavej strany rovnosti (6) do rovnosti (7) dostaneme výsledok:

$2=196\cdot 4+34\cdot (-23)(8)$

t.j. jednu z možností, ako vyjadriť najväčší spoločný deliteľ dvoch čísel ich lineárnou kombináciou.

Odvodenie

Nech $k$ je iterácia v Euklidovom algoritme, v ktorej bol nájdený najväčší spoločný deliteľ dvoch kladných celých čísel $d=NSD(m,n)=d_k$, čiže pre ktoré platí $z_k=0$. Spätným dosadením vyjadríme najväčší spoločný deliteľ $d$ ako lineárnu kombináciu $c_i$ a $d_i$ postupne pre každé $i=k,k-1,...,0$, t.j:

$d=c_i{u}_i+d_i{v}_i(9)$

Dosadením (2) a (3) do (9) dostaneme:

$d=d_{i-1}{u}_i+z_{i-1}{v}_i(10)$

a dosadením (1) do (10):

$\begin{array}{cc}d& =d_{i-1}{u}_i+(c_{i-1}-d_{i-1}{p}_{i-1})v_i\\ & =c_{i-1}{v}_i+d_{i-1}(u_i-p_{i-1}{v}_i)(11)\end{array}$

Porovnaním (11) a (9) dostávame výsledok:

$\begin{array}{cccc}{u}_{i-1}& =& v_i& (12)\\ v_{i-1}& =& u_i-p_{i-1}{v}_i& (13)\end{array}$

Tiež platí:

$d=d_k=c_k\cdot 0+d_k\cdot 1$

preto:

$\begin{array}{cccc}{u}_k& =& 0& (14)\\ v_k& =& 1& (15)\end{array}$

Formálny zápis

1. [Inicializácia premenných.] Priraďte i ← k, u ← 0, v ← 1.

2. [Ukončovacia podmienka.] Ak i = 0 tak koniec, pričom NSD(m,n) = m·u + n·v

3. [Výpočet u a v.] Priraďte u_new ← v, v_new ← u - p[i-1]·v.

4. [Ďalšia iterácia.] Priraďte i ← i-1, u ← u_new, v ← v_new. Vráťte sa do bodu 2.

Poznámka: Výpočet prebieha opačným smerom než výpočet NSD a okrem naposledy vypočítanej dvojice hodnôt $u_i$ a $v_i$ je potrebné mať k dispozícii aj všetky hodnoty $p_i$.

Použitie na konkrétnom príklade

i	ci	di	pi	zi	ui	vi
0	196	34	5	26	4	-23
1	34	26	1	8	-3	4
2	26	8	3	2	1	-3
3	8	2	4	0	0	1

Rozšírený Euklidov algoritmus

Odvodenie

Nech $k$ je iterácia v Euklidovom algoritme, v ktorej bol nájdený najväčší spoločný deliteľ dvoch kladných celých čísel $d=NSD(m,n)=d_k$, čiže pre ktoré platí $z_k=0$. Pre každé $i=0,1,2,...,k-1$ vyjadrime $z_i$ ako lineárnu kombináciu čísel $m$ a $n$, t.j. hľadáme dvojicu celých čísel $a_i$ a $b_i$, pre ktoré platí $z_i=m{a}_i+n{b}_i(16)$

• Pre $i=0$ z rovnosti (1) vyplýva:

$\begin{array}{ccc}{z}_0& =& c_0-d_0{p}_0\\ & =& m-n{p}_0& (17)\end{array}$

preto

$\begin{array}{cccc}{a}_0& =& 1& (18)\\ b_0& =& -p_0& (19)\end{array}$

• Pre $i=1$ z (1) vyplýva:

$\begin{array}{cccc}{z}_1& =& c_1-d_1{p}_1& (20)\end{array}$

a podľa (2) a (3):

$\begin{array}{ccccc}{c}_1& =& d_0& =& n\\ d_1& =& z_0\end{array}$

Preto

$\begin{array}{ccc}{z}_1& =& n-z_0{p}_1\end{array}$

Po dosadení za $z_0$ zo (17) dostaneme:

$\begin{array}{ccc}{z}_1& =& n-m{p}_1+n{p}_0{p}_1\\ & =& -m{p}_1+n(1+p_0{p}_1)& (21)\end{array}$

A teda

$\begin{array}{cccc}{a}_1& =& -p_1& (22)\\ b_1& =& 1+p_0{p}_1& (23)\end{array}$

• Pre ľubovoľné $i>1$ je podľa (2) a (3):

$\begin{array}{ccccc}{c}_i& =& d_{i-1}& =& z_{i-2}\\ d_i& =& z_{i-1}\end{array}$

Dosadením do (1) dostávame:

$\begin{array}{ccccc}{z}_i& =& c_i-d_i{p}_i& =& z_{i-2}-z_{i-1}{p}_i\end{array}$

Ak poznáme (16) pre $i-1$ a $i-2$, dostávame:

$\begin{array}{ccc}{z}_i& =& m{a}_{i-2}+n{b}_{i-2}-(m{a}_{i-1}+n{b}_{i-1})p_i\\ & =& m(a_{i-2}-a_{i-1}{p}_i)+n(b_{i-2}-b_{i-1}{p}_i)\end{array}$

Preto pre ľubovoľné $i>1$ je:

$\begin{array}{cccc}{a}_i& =& a_{i-2}-a_{i-1}{p}_i& (24)\\ b_i& =& b_{i-2}-b_{i-1}{p}_i& (25)\end{array}$

• Ak sa nasledujúcim spôsobom zadefinujú hodnoty $z_i$, $a_i$ a $b_i$ aj pre $i=-1$ a $i=-2$:

$\begin{array}{ccccccccc}{z}_{-2}& =& m,& a_{-2}& =& 1,& b_{-2}& =& 0\\ z_{-1}& =& n,& a_{-1}& =& 0,& b_{-1}& =& 1\end{array}$

bude rovnosť (16) platiť aj pre ne a vzťahy (24), (25) je možné použiť aj na zistenie $a_0,a_1,b_0,b_1$.

Formálny zápis

1. [Inicializácia premenných.] Priraďte c ← m, d ← n, a_old ← 1, b_old ← 0, a ← 0, b ← 1.

2. [Výpočet podielu a zvyšku.] Vypočítajte celočíselný podiel p a zvyšok z po delení c / d.

3. [Ukončovacia podmienka.] Ak z = 0 tak koniec, pričom NSD(m,n) = d = m·a + n·b

4. [Výpočet a a b.] Priraďte a_new ← a_old - a·p, b_new ← b_old - b·p.

5. [Redukcia.] Priraďte c ← d, d ← z, a_old ← a, a ← a_new, b_old ← b, b ← b_new. Vráťte sa do bodu 2.

Poznámka: Výpočet prebieha súčasne s výpočtom NSD pričom je potrebné mať k dispozícii aj dve posledné dvojice hodnôt $a_i,b_i$.

Použitie na konkrétnom príklade

i	ci	di	pi	zi	ai	bi
-2				196	1	0
-1				34	0	1
0	196	34	5	26	1	-5
1	34	26	1	8	-1	6
2	26	8	3	2	4	-23
3	8	2	4	0

Zdroj

• Šablóna:Citácia knihy