Faktorenzerlegung grober Zahlen

Faktorenzerlegung großer Zahlen

Faktorenzerlegung à la Monte Carlo :

Bei der Methode von J. M. Pollard ist es nicht immer möglich bei einer teilbaren Zahl einen Faktor zu finden. Die Suche wird zu einem Glücksspiel. Man untersucht ob es von einer natürlichen Zahl N und einer errechneten Zahl P einen größten gemeinsamen Teiler gibt. P bekommt man, indem zwei neue Variablen (x y ) einführt, für welche gilt : f(x)=x²+c (c ) x x²+c ; y (y²+c)²+c »»»» P =P ·( y - x )

Als Startwert für x , y , und P wird verwendet.

1.Schritt:

Setze x=1 ;y=1 ;P=1

2.Schritt:

Setze x = x²+c ; y = (y²+c)²+c und P =P ·( y - x )

3.Schritt:

Berechne den ggT von P und N. Fals das Ergebnis = 1 ist , Wiederhole den Vorgang ab dem Schritt 2 .

Bei allen anderen Ergebnissen hat man mit dem ggT den gesuchten Teiler gefunden .

Achtung

Bei Schritt 2 entstehen für die Variablen x, y und P große Werte. Man muß daher die Ergebnisse                

modulo N nehmen .

Fals N selbst der Teiler ist, gibt es zwei Möglichkeiten:

N ist eine Primzahl

Man kann die Konstante c verändern (z.B.:von 1 auf 2,3,)

Wenn einem die Berechnung des Schrittes 3 zu Zeitaufwendig ist, kann man dies umgehen, indem

man erst nach z.B.: 10 -maliger Wiederholung des 2. Schrittes mit dem 3.Schritt beginnt .(Nur sehr

seltener Verlust des Teilers ) .

Hier ein Programmierbeispiel in Q-Basic :

CLS:CLEAR: X=1: Y=1: P=1: T=1 ' **** Programm zur Faktorenzerlegung von MEISEL Marcus ****

INPUT "zu teilende Zahl (N): ";N : MO=N: INPUT "Konstante C: "; C: CLS

X=X +C: GOSUB 9: Y=(Y +C) +C: GOSUB 10: P=P*(Y-X):GOSUB 11: V=V+1: IF V>=25 THEN V=25

LOCATE 10,0:PRINT "X Y P":LOCATE 10,V:PRINT X;:LOCATE 23,V:PRINT Y;: LOCATE 36,V:PRINT P;

IF P=0 THEN GOTO 19

IF P=1 THEN GOTO 3

INPUT "";I:GOTO 12

GOTO 3

X=X-(INT(X/MO))*MO:RETURN

Y=Y-(INT(Y/MO))*MO:RETURN

IF ABS ( P ) < > P THEN P= N - ABS ( P ) : P=P - ( INT ( P / MO ) ) * MO : RETURN

' *****ggT --ausrechnen

TE = P

RE = N - ( INT ( N / TE ) ) * TE: IF RE=1 THEN T=1: GOTO 3

IF RE=0 THEN T=TE: GOTO 17

N=TE:TE=RE:GOTO 14

' *****die Lösung

CLS:LOCATE 1,1:PRINT "Zahl ";MO;" u. ";P;" sind durch ";T;" Teilbar.";:BEEP:INPUT"" ;I:GOTO1

' *****eine Primzahl

CLS:LOCATE 1,1:PRINT "Zahl ";MO;"ist eine Primzahl! -oder ein anderes !!c!! prob.";:BEEP:INPUT"" ;I

GOTO1 ' ***** ENDE dieses ©Programs

Hier ein RechenBeispiel:

N=143 => X Y P

c=1 1 1 1

2 5 3

5 105 14

26 83 83

105 105 0=> Primzahl oder c anders wählen z.B.: c=2

N=143

c=2

=> X Y P

3 11 8

11 116 125

123 115 142

116 38 ggT(65,143)=13 =>13/143

Antwort : 143 ist durch 13 teilbar.

4.) Weitere Algorithmen :

Ein weiterer Algorithmus von J. M. Pollard:

Es wird ein Teiler p von einer Zahl N gesucht. Er ist gefunden, wenn

gilt : p-1 ist nur durch kleine Primfaktoren teilbar. Auf diesen

Rechengang wird in dem Artikel von Williams näher eingegangen.

Der SQUFOF-Algorithmus von D. Shanks:

Hierfür benötigt man die Theorie der Ideale in quadratischen

Zahlkörpern.Näheres : siehe Artikel von Williams.