Pollards p–1-methode

Uit testwiki
Versie door imported>Jcvg16 op 18 apr 2023 om 14:56 (growthexperiments-addlink-summary-summary:3|0|0)
(wijz) ← Oudere versie | Huidige versie (wijz) | Nieuwere versie → (wijz)
Naar navigatie springen Naar zoeken springen

Pollards p–1-methode is een methode voor het ontbinden van een geheel getal in priemfactoren. In 1974 publiceerde John Pollard zijn algoritme voor redelijk grote getallen. Deze getallen moeten zodanig zijn dat van elke priemfactor p de voorganger p1 een glad getal is. Als een getal aan deze voorwaarde voldoet, kan met Pollards p–1-methode een priemfactor van dit getal worden gevonden.

Gladheid

De eis voor gladheid van de voorganger p1 van een priemfactor p houdt niet alleen in dat p1 B-glad is, maar ook dat de voorkomende machten van de priemfactoren van p1 niet groter zijn dan B.

Basisgedachten bij Pollards algoritme

Er liggen twee basisgedachten ten grondslag aan Pollards p–1-methode.

  • Als p een priemfactor is van het samengestelde gehele getal n, dan geldt volgens de kleine stelling van Fermat voor alle gehele getallen a relatief priem met p en voor alle positieve gehele getallen k:
ak(p1)1(modp)
  • Als een getal x congruent is aan 1 modulo een factor van n, dan is de grootste gemene deler van x1 en n deelbaar door deze factor.

Dat leidt ertoe dat

ak1(modp),

als p een deler is van n, en p1 een deler van k.

Hieruit volgt weer dat als ak1 deelbaar is door p en ook n deelbaar is door p, dat dan ook de grootste gemene deler van ak1 en n deelbaar is door p.

Het idee is om voor de exponent een groot veelvoud van p1 te nemen, door een getal met zeer veel priemfactoren te kiezen. Over het algemeen kiest men het product van machten van alle priemgetallen kleiner dan een bepaalde grens B. Begin met een willekeurige x, en bepaal iteratief een nieuwe x als xw(modn), waarbij w loopt door de machten van deze priemgetallen. Controleer in elk stadium, of eventueel aan het eind, of ggd(x1,n)1.

Het algoritme van Pollard

Het algoritme van Pollard om het getal n te factoriseren werkt als volgt:

  1. Kies een niet te groot en ook niet te klein getal B als grens voor de exponent.
  2. Bereken het kleinste gemene veelvoud van alle getallen B
  3. Kies een willeleurig getal a<n1
  4. Bepaal ak(modn)
  5. Bepaal de grootste gemene deler d van ak1 en n
  6. Als 1<d<n, dan is d een deler van n
  7. Ontbind n

Als het algoritme niet meteen een priemfactor vindt, is het mogelijk te variëren met de keuze van B en a.

Voorbeeld

Hoe kunnen we het getal 540143 ontbinden in priemfactoren met behulp van het algoritme van Pollard?

  1. Kies B=8
  2. k=kgv(2,3,4,5,6,7,8)=840
  3. Kies a=2
  4. 284053047(mod540143) (zie hieronder uitwerking 1)
  5. d=ggd(53046,540143)=421 (zie hieronder uitwerking 2)
  6. 1<421<540143, dus 421|540143
  7. 540143 = 421 × 1283

Uitwerking 1

Bepaal 2840(mod540143)

2840=2512×2256×264×28
28=256
232=4294967296290303(mod540143)
2642903032=8427583180920234(mod540143)
2128202342=409414756526505(mod540143)
22565265052=277207515025185852(mod540143)
25121858522=34540965904441483(mod540143)
2840=2512×2256×264×28=42501370546502246453047(mod540143)

Uitwerking 2

Bepaal de ggd van 540143 en 53046 met het algoritme van Euclides

540143 = 10 × 53046 + 9683
53046 = 5 × 9683 + 4631
9683 = 2 × 4631 + 421
4631 = 11 × 421

Dus ggd( 540143 , 53046 ) = 421

Bronnen