Euklids algoritme

Fra Wikipedia, den frie encyklopædi
Spring til navigation Spring til søgning
Euklid som den flamske maler Justus van Gent (ca. 1410 - ca. 1480) forestillede sig ham omkring 1474.

Euklids algoritme er en matematisk algoritme. Det er en effektiv metode til at beregne den største fælles divisor (forkortet SFD eller GCD efter greatest common divisor). For to tal er SFD det største tal, der går op i begge tal. Algoritmen er opkaldt efter den græske matematiker Euklid, der først beskrev det i sine Elementer (ca. 300 f.Kr.). Det er et eksempel på en algoritme, en trin-for-trin-procedure til udførelse af en beregning i henhold til veldefinerede regler, og er en af de ældste algoritmer i almindelig brug. Den kan bruges til at forkorte brøker så meget som muligt og er blot en af mange talteoretiske og kryptografiske beregningsmetoder. En vigtig fordel ved Euklids algoritme er, at den kan finde SFD effektivt uden at skulle beregne primtalsfaktorerne.[1] Faktorisering af store heltal antages at være et beregningsmæssigt meget vanskeligt problem, og sikkerheden i mange bredt anvendte kryptografiske protokoller er baseret på, at det er uoverkommeligt.[2] Euklids algoritme, der beregner SFD for to heltal, er tilstrækkelig til at beregne SFD for vilkårligt mange heltal.[kilde mangler]

Algoritmen[redigér | redigér wikikode]

Euklids algoritme er baseret på princippet om, at den største fælles divisor af to tal ikke ændres, hvis det større tal erstattes af forskellen mellem de to tal. For eksempel er 21 SFD for 252 og 105 (eftersom 252 = 21 × 12 og 105 = 21 × 5), og det samme tal, 21, er også SFD for 105 og 252 - 105 = 147. Da denne erstatning formindsker det største af de to tal, giver gentagelse af denne proces successivt mindre talpar, indtil de to tal bliver ens. Når dette sker, er de SFD for de oprindelige to tal. Ved at invertere processen, kan SFD udtrykkes som en sum af de to originale tal, der ganges med et positivt eller negativt heltal, f.eks. 21 = 5 × 105 + (−2) × 252. Det, at SFD altid kan udtrykkes på denne måde, er kendt som Bézouts identitet.[kilde mangler]

Den version af Euklids algoritme, der er beskrevet ovenfor (og af Euklid), kan tage mange trin for at finde SFD, når et af de givne tal er meget større end det andet. En mere effektiv version af algoritmen tager en genvej, hvor man i stedet erstatter det største af de to tal med dets rest, når det divideres med den mindste af de to (i denne version stopper algoritmen, når det største tal er deleligt med det mindste). Med denne forbedring kræver algoritmen aldrig flere trin end fem gange antallet af cifre (base 10) i det mindste heltal. Dette blev bevist af Gabriel Lamé i 1844 og markerer begyndelsen på kompleksitetsteori . Yderligere metoder til forbedring af algoritmens effektivitet blev udviklet i det 20. århundrede.[kilde mangler]

Procedure[redigér | redigér wikikode]

Euklids algoritme er en iterativ metode, således at outputtet fra hvert trin bruges som input til det næste. Lad være et helt tal, der tæller trinnene i algoritmen, og start med at tælle fra nul. Det første trin svarer således til , det næste trin svarer til osv.

Hvert trin begynder med to ikke-negative rester og . Da algoritmen sikrer, at resterne falder støt med hvert trin, er mindre end dens forgænger . Målet med 'ne skridt er at finde en kvotient og resten , der opfylder ligningen

og hvor der skal gælde at

Med andre ord fratrækkes det mindste tal gentagne gange fra det største tal , indtil resten er mindre end .

I det indledende trin () er resterne og lig med og , de tal, som SFD søges for. I det næste trin () er resterne og resten i det indledende trin, og så videre. Algoritmen kan således skrives som en sekvens af ligninger

Hvis er mindre end , vil det første trin i algoritmen bytte om på tallene. For eksempel, hvis vil den oprindelige kvotient være nul, og resten være . Således er mindre end sin forgænger for alle .

Da resterne falder i hvert trin, men aldrig kan være negative, vil en rest til sidst blive lig med nul og så stopper algoritmen.[3] Den sidste , som er forskellig fra nul, er den største fælles divisor af og . Tallet kan ikke være uendeligt, fordi der kun er et begrænset antal ikke-negative heltal mellem den oprindelige rest og 0.[kilde mangler]

Bevis for korrekthed[redigér | redigér wikikode]

Gyldigheden af Euklids algoritme kan bevises ved et to-trins argument.[4] I det første trin vises, at den sidste rest , der altså er forskellig fra nul, går op i både og . Da det er en fælles divisor, skal den være mindre end eller lig med den største fælles divisor . I det andet trin vises det, at enhver fælles divisor af og , inklusiv , skal gå op i , og derfor skal være mindre end eller lig med . Disse to konklusioner betyder sammen, at

For at vise at går op i både og (det første trin), ses først, at går op i sin forgænger

da den sidste rest er nul. går også op i den næste forgænger

fordi den går op i begge led på højre side af ligningen. Ved at gentage dette argument igen og igen, ses det at går op i alle de foregående rester, inklusive og . Ingen af de foregående rester , osv. går op i og , da de efterlader en rest. Da er en fælles divisor af og , må

I det andet trin vises det at ethvert naturligt tal c der går op i både og (med andre ord enhver fælles divisor af og ), også går op i alle resterne . Pr. definition kan og skrives som multipla af :

hvor og er naturlige tal. Derfor går op i den indledende rest , da

Et analogt argument viser, at også går op i de efterfølgende rester , osv. Derfor skal den største fælles divisor gå op i , hvilket indebærer, at

Da den første del af argumentet viste det omvendte (), følger det, at

Således er den største fælles divisor af alle de efterfølgende par:[5]

Gennemregnet eksempel[redigér | redigér wikikode]

Animation in which progressively smaller square tiles are added to cover a rectangle completely.
Animation af Euklids algoritme baseret på substraktion. Det indledende rektangel (grønt) har dimensionerne og . Kvadrater (orange) i størrelsen anbringes inden i, hvilket efterlader et rektangel. Dette rektangel belægges med kvadrater (blå), indtil der er et rektangel tilbage, som så igen belægges, denne gang med kvadrater (røde), hvilket ikke efterlader noget udækket område. Den mindste firkantede størrelse, 21, er dermed SFD for 1071 og 462.

Til illustration kan Euklids algoritme bruges til at finde den største fælles divisor af og . Til at begynde med trækkes 462 fra 1071, indtil resten er mindre end 462. Dettte kan gøres to gange (), hvilket efterlader en rest af 147:

Derefter trækkes 147 fra 462, indtil resten er mindre end 147. Det kan gøres tre gange ( q 1   =   3), hvilket efterlader en rest af 21:

462 = 3 × 147 + 21.

Derefter trækkes 21 fra 147, indtil resten er mindre end 21. Dette kan gøres syv gange ( q 2   =   7) uden at efterlade nogen rest:

147 = 7 × 21 + 0.

Da den sidste rest er nul, slutter algoritmen med 21 som den største fælles divisor for 1071 og 462. Dette stemmer overens med sfd(1071, 462), som beregnet overfor ved hjælp af primfaktorisering. I tabelform er trinnene:

Trin Ligning Kvotient og resten
0 1071 = q0 462 + r0 q0 = 2 og r0 = 147
1 462 = q1 147 + r1 q1 = 3 og r1 = 21
2 147 = q2 21 + r2 q2 = 7 og r2 = 0; algoritmen slutter

Visualisering[redigér | redigér wikikode]

Euklids algoritme kan visualiseres ved hjælp af den flisebelagte analogi, der er givet ovenfor for den største fælles divisor. Antag, at vi ønsker at dække et a x b rektangel med firkantede fliser nøjagtigt, hvor a er det største af de to tal. Vi forsøger først at fliselægge rektanglet ved hjælp af kvadratiske b x b fliser, men dette efterlader en resterende r0 x b rektangel ubelagt, hvor r 0   <   b . Vi forsøger derefter at fliselægge den resterende rektangel med kvadratiske r0 x r0 fliser. Dette efterlader et andet rektangel r1 x r0 ubelagt, som vi forsøger at fliselægge vha. kradratiske r1 x r 1 fliser, og så videre. Sekvensen slutter, når der ikke er noget resterende rektangel, dvs. når de kvadratiske fliser dækker hele det forrige rektangel. Sidelængden på den mindste flise er SFD for dimensionerne på det originale rektangel. For eksempel er den mindste firkantede flise i den tilstødende figur 21 x 21 (vist i rødt), og 21 er SFD for 1071 og 462, dimensionerne på det originale rektangel (vist i grønt).

Division med rest[redigér | redigér wikikode]

Ved hvert trin k beregner Euklids algoritme en kvotient og resten fra to tal og

hvor er ikke-negativ og strengt mindre end størrelsen. Dette kaldes division med rest. Den matematiske sætning der ligger bag division med rest, er at en sådan kvotient og resten altid findes og er unikke.[6]

I Euklids originale version af algoritmen findes kvotienten og resten ved gentagen subtraktion; dvs. trækkes fra , indtil resten er mindre end . Derefter byttes der om på og (så næste gang trækkes fra ), og processen gentages. Euklidisk opdeling reducerer alle trin mellem to opbytninger til et enkelt trin, hvilket er mere effektivt. Desuden er kvotienterne ikke nødvendige, så man kan erstatte division med rest med modulooperationen, der kun giver resten. Således bliver hver iteration af Euklids algoritme simpelthen

Implementeringer[redigér | redigér wikikode]

Implementeringer af algoritmen kan udtrykkes i pseudocode . F.eks. kan den divisionsbaserede version programmeres som [7]

function sfd (a, b)
    while b ≠ 0
        t := b; 
        b := a mod b; 
        a := t; 
    return a; 

I begyndelsen af den k'te iteration indeholder variablen b den nyeste rest rk-1, hvorimod variablen a holder dens forgænger, rk-2. Trinnet b := a mod b svarer til ovennævnte rekursionsformel r kr k −2 mod r k −1 . Den midlertidige variabel t holder værdien af rk-1 mens den næste rest r k beregnes. Ved slutningen af løkke-iterationen indeholder variablen b resten rk, mens variablen a holder sin forgænger, rk−1 .

I den subtraktionsbaserede version, der var Euklids originale version, er beregningen b = a mod b erstattet af gentagen subtraktion.[8] I modsætning til den divisionsbaserede version, der fungerer med vilkårlige heltal som input, antager den subtraktionsbaserede version, at input består af positive heltal og stopper, når a = b :

function sfd(a,b)
    while a ≠ b
        if a > b
            a := a - b; 
        else
            b := b - a; 
    return a; 

Variablerne a og b veksler med de foregående rester r k −1 og r k −2 . Hvis a er større end b i begyndelsen af en iteration, så er a = r k -2, eftersom r k -2 > r k -1. I løbet af løkkens iteration, bliver a reduceret med den tidligere rest b, indtil a er mindre end b. Så er a den næste rest r k . Derefter reduceres b med a, indtil det igen er mindre end a, hvilket giver den næste rest r k +1, og så videre.

Den rekursive version [9] er baseret på at to rester efter hinanden alle har samme GDC, samt stoptilstanden sfd ( r N −1,   0)   =   r N −1 .

function gcd (a, b)
    if b = 0
        return a; 
    else
        return gcd(b, a mod b);

Anvendelser[redigér | redigér wikikode]

Euklids algoritme har mange teoretiske og praktiske anvendelser. Den bruges til at gøre brøker uforkortelige og til at foretage division i kongruensregning. Beregninger, der bruger denne algoritme, indgår i kryptografiske protokoller, der bruges til at sikre internetkommunikation, og i metoder til at bryde disse kryptosystemer ved at primtalsfaktorisere store tal. Euklids algoritme kan bruges til at løse Diofantiske ligninger såsom at finde tal, der tilfredsstiller flere kongruenser samtidig som i den kinesiske restklassesætning til at konstruere kædebrøker og til at finde gode rationelle tilnærmelser til reelle tal. Endelig kan den bruges som et grundlæggende værktøj til at bevise sætninger i talteori såsom Lagranges fire-kvadratsætning, og at primtalsfaktoriseringer er unikke. Den originale algoritme blev kun beskrevet for naturlige tal og geometriske længder (reelle tal), men algoritmen blev generaliseret i det 19. århundrede til andre typer tal såsom Gaussiske heltal og polynomier i en variabel. Dette førte til udviklingen af koncepter i moderne abstrakt algebra såsom euklidiske ringe.[kilde mangler]

Kildehenvisninger[redigér | redigér wikikode]

  1. ^ Schroeder 2005, s. 19
  2. ^ Schroeder 2005, s. 216–219
  3. ^ Stark 1978, s. 18
  4. ^ Stark 1978, s. 16–20
  5. ^ Knuth 1997, s. 320
  6. ^ Dummit, David S.; Foote, Richard M. (2004). Abstract Algebra. John Wiley & Sons, Inc. s. 270-271. ISBN 978-0-471-43334-7. 
  7. ^ Knuth 1997, s. 319–320
  8. ^ Knuth 1997, s. 318–319
  9. ^ Stillwell 1997, s. 14

Litteratur[redigér | redigér wikikode]

Oversættelse
Denne artikel eller en tidligere version er helt eller delvist oversat fra den engelsksprogede Wikipedia, der er tilgængelig under Creative Commons Kreditering-Deling på samme vilkår 3.0. Se versionshistorik for oplysninger om oprindelig(e) bidragyder(e).