Algorisme d'Euclides ampliat

L'algorisme d'Euclides ampliat o algorisme d'Euclides estès és una millora de l'algorisme d'Euclides de càlcul del màxim comú divisor de dos nombres enters, que dona, a més del màxim comú divisor dels dos nombres, els coeficients de cadascun d'aquests dos nombres a la identitat de Bézout.

Descripció

Siguin a {\displaystyle a} i b 0 {\displaystyle b\neq 0} dos nombres enters. L'algorisme d'Euclides consisteix a construir la recurrència finita

{ d 1 = | a | d 2 = | b | d i = d i 2 d i 1 q i , i > 2 , 0 < d i < d i 1 {\displaystyle {\begin{cases}d_{1}=|a|\\d_{2}=|b|\\d_{i}=d_{i-2}-d_{i-1}q_{i}\,,\quad i>2\,,\quad 0<d_{i}<d_{i-1}\end{cases}}}

en la qual d i = d i 2 d i 1 q i {\displaystyle d_{i}=d_{i-2}-d_{i-1}q_{i}} no és més que el residu de la divisió entera de d i 2 {\displaystyle d_{i-2}} i d i 1 {\displaystyle d_{i-1}} amb quocient q i {\displaystyle q_{i}} . La successió és estrictament decreixent i la condició d i > 0 {\displaystyle d_{i}>0} obliga a que sigui finita. L'últim terme, posem d k {\displaystyle d_{k}} arriba quan hi ha q {\displaystyle q} que fa d k 1 = d k q {\displaystyle d_{k-1}=d_{k}q} . La successió té, doncs, k {\displaystyle k} termes i d k = m . c . d . ( a , b ) {\displaystyle d_{k}=m.c.d.(a,b)} .

Però si ara considerem aquestes altres dues recurrències finites:

{ x 1 = 1 x 2 = 0 x i = x i 2 x i 1 q i , k i > 2 { y 1 = 0 y 2 = 1 y i = y i 2 y i 1 q i , k i > 2 {\displaystyle {\begin{cases}x_{1}=1\\x_{2}=0\\x_{i}=x_{i-2}-x_{i-1}q_{i}\,,\quad k\geq i>2\end{cases}}\qquad {\begin{cases}y_{1}=0\\y_{2}=1\\y_{i}=y_{i-2}-y_{i-1}q_{i}\,,\quad k\geq i>2\end{cases}}}

amb els valors q i {\displaystyle q_{i}} de la successió de l'algorisme d'Euclides, resulta que, per i > 2 {\displaystyle i>2} amb 0 < d i < d i 1 {\displaystyle 0<d_{i}<d_{i-1}} , es té

d i = d 1 x i + d 2 y i {\displaystyle d_{i}=d_{1}x_{i}+d_{2}y_{i}}

com es comprova fàcilment per inducció.

Per tant, si d k = m . c . d . ( a , b ) {\displaystyle d_{k}=m.c.d.(a,b)} , resulta

m . c . d . ( a , b ) = | a | x k + | b | y k {\displaystyle m.c.d.(a,b)=|a|x_{k}+|b|y_{k}}

i x k {\displaystyle x_{k}} i y k {\displaystyle y_{k}} , amb els signes adequats, són els coeficients de a {\displaystyle a} i b {\displaystyle b} a la identitat de Bézout.

Càlcul pràctic

Hom sol disposar els càlculs en una graella com aquesta

1 {\displaystyle 1} 0 {\displaystyle 0} x 3 {\displaystyle x_{3}} x 4 {\displaystyle x_{4}} x 5 x i 2 {\displaystyle x_{5}\ldots \ldots x_{i-2}} x i 1 {\displaystyle x_{i-1}} x i x k 2 {\displaystyle x_{i}\ldots \ldots x_{k-2}} x k 1 {\displaystyle x_{k-1}} x k {\displaystyle x_{k}}
0 {\displaystyle 0} 1 {\displaystyle 1} y 3 {\displaystyle y_{3}} y 4 {\displaystyle y_{4}} y 5 y i 2 {\displaystyle y_{5}\ldots \ldots y_{i-2}} y i 1 {\displaystyle y_{i-1}} y i y k 2 {\displaystyle y_{i}\ldots \ldots y_{k-2}} y k 1 {\displaystyle y_{k-1}} y k {\displaystyle y_{k}}
q 3 {\displaystyle q_{3}} q 4 {\displaystyle q_{4}} q 5 {\displaystyle q_{5}} q 6 q i 1 {\displaystyle q_{6}\ldots \ldots q_{i-1}} q i {\displaystyle q_{i}} q i + 1 q k 1 {\displaystyle q_{i+1}\ldots \ldots q_{k-1}} q k {\displaystyle q_{k}} q {\displaystyle q}
| a | {\displaystyle |a|} | b | {\displaystyle |b|} d 3 {\displaystyle d_{3}} d 4 {\displaystyle d_{4}} d 5 d i 2 {\displaystyle d_{5}\ldots \ldots d_{i-2}} d i 1 {\displaystyle d_{i-1}} d i d k 2 {\displaystyle d_{i}\ldots \ldots d_{k-2}} d k 1 {\displaystyle d_{k-1}} d k {\displaystyle d_{k}} 0 {\displaystyle 0}

Hom comença obtenint q 3 {\displaystyle q_{3}} com a quocient de la divisió entera de d 1 {\displaystyle d_{1}} entre d 2 {\displaystyle d_{2}} , és a dir, | a | {\displaystyle |a|} entre | b | {\displaystyle |b|} i d 3 {\displaystyle d_{3}} a partir de d 3 = d 1 d 2 q 3 = | a | | b | q 3 {\displaystyle d_{3}=d_{1}-d_{2}q_{3}=|a|-|b|q_{3}} . Els termes x 3 {\displaystyle x_{3}} i y 3 {\displaystyle y_{3}} resulten de x 3 = x 1 x 2 q 3 = 1 {\displaystyle x_{3}=x_{1}-x_{2}q_{3}=1} i y 3 = y 1 y 2 q 3 = q 3 {\displaystyle y_{3}=y_{1}-y_{2}q_{3}=-q_{3}} . Els termes següents, q i {\displaystyle q_{i}} , d i {\displaystyle d_{i}} , x i {\displaystyle x_{i}} i y i {\displaystyle y_{i}} s'obtenen de la mateixa manera i en el mateix ordre:

{ q i   e ´ s el quocient de la divisi o ´  entera de  d i 2  entre  d i 1 d i = d i 2 d i 1 q i   ( e ´ s el residu de la divisi o ´  entera de  d i 2  entre  d i 1 ) x i = x i 2 x i 1 q i y i = y i 2 y i 1 q i {\displaystyle {\begin{cases}q_{i}{\mbox{ }}{\acute {\mbox{e}}}{\mbox{s el quocient de la divisi}}{\acute {\mbox{o}}}{\mbox{ entera de }}d_{i-2}{\mbox{ entre }}d_{i-1}\\d_{i}=d_{i-2}-d_{i-1}q_{i}{\mbox{ }}({\acute {\mbox{e}}}{\mbox{s el residu de la divisi}}{\acute {\mbox{o}}}{\mbox{ entera de }}d_{i-2}{\mbox{ entre }}d_{i-1})\\x_{i}=x_{i-2}-x_{i-1}q_{i}\\y_{i}=y_{i-2}-y_{i-1}q_{i}\end{cases}}}

i el procés acaba quan trobem d k 1 = d k q {\displaystyle d_{k-1}=d_{k}q} . Aleshores, m . c . d . ( a , b ) = d k = | a | x k + | b | y k {\displaystyle m.c.d.(a,b)=d_{k}=|a|x_{k}+|b|y_{k}}

Exemple

Il·lustrem aquest procés amb un exemple: es tracta de calcular m . c . d . ( 763 , 175 ) {\displaystyle m.c.d.(763,175)} :

1 {\displaystyle 1} 0 {\displaystyle 0} 1 {\displaystyle 1} 2 {\displaystyle -2} 3 {\displaystyle 3} 11 {\displaystyle -11}
0 {\displaystyle 0} 1 {\displaystyle 1} 4 {\displaystyle -4} 9 {\displaystyle 9} 13 {\displaystyle -13} 48 {\displaystyle 48}
4 {\displaystyle 4} 2 {\displaystyle 2} 1 {\displaystyle 1} 3 {\displaystyle 3} 2 {\displaystyle 2}
763 {\displaystyle 763} 175 {\displaystyle 175} 63 {\displaystyle 63} 49 {\displaystyle 49} 14 {\displaystyle 14} 7 {\displaystyle 7} 0 {\displaystyle 0}

que prové de

1 {\displaystyle 1} 0 {\displaystyle 0} 1 = 1 4 0 {\displaystyle 1=1-4\cdot 0} 2 = 0 2 1 {\displaystyle -2=0-2\cdot 1} 3 = 1 1 ( 2 ) {\displaystyle 3=1-1\cdot (-2)} 11 = 2 3 3 {\displaystyle -11=-2-3\cdot 3}
0 {\displaystyle 0} 1 {\displaystyle 1} 4 = 0 4 1 {\displaystyle -4=0-4\cdot 1} 9 = 1 2 ( 4 ) {\displaystyle 9=1-2\cdot (-4)} 13 = 4 1 9 {\displaystyle -13=-4-1\cdot 9} 48 = 9 3 ( 13 ) {\displaystyle 48=9-3\cdot (-13)}
4 = 763 ÷ 175 {\displaystyle 4=763\div 175} 2 = 175 ÷ 63 {\displaystyle 2=175\div 63} 1 = 63 ÷ 49 {\displaystyle 1=63\div 49} 3 = 49 ÷ 14 {\displaystyle 3=49\div 14} 2 = 14 ÷ 7 {\displaystyle 2=14\div 7}
763 {\displaystyle 763} 175 {\displaystyle 175} 63 = 763 175 4 {\displaystyle 63=763-175\cdot 4} 49 = 175 63 2 {\displaystyle 49=175-63\cdot 2} 14 = 63 1 49 {\displaystyle 14=63-1\cdot 49} 7 = 49 3 14 {\displaystyle 7=49-3\cdot 14} 0 = 14 2 7 {\displaystyle 0=14-2\cdot 7}

(Les divisions a ÷ b {\displaystyle a\div b} se sobreentenen enteres) Aleshores, m . c . d . ( 763 , 175 ) = 7 = 763 ( 11 ) + 175 48 {\displaystyle m.c.d.(763,175)=7=763\cdot (-11)+175\cdot 48} .

Referències

  • PlanetMath: Euclid's algorithm (anglès)