Solovay-Strassen-Test

Der Solovay-Strassen-Test (nach Robert Solovay und Volker Strassen) ist ein probabilistischer Primzahltest. Der Test prüft für eine ungerade Zahl n, ob sie eine Primzahl oder das Produkt mehrerer Faktoren ist. Für den Fall, dass der Test als Antwort liefert, dass es sich bei n nicht um eine Primzahl handelt, liefert der Test jedoch im Allgemeinen keinen Faktor der Zahl n.

Der Test ist ein probabilistischer Algorithmus. Das heißt, das Testergebnis kann mit einer gewissen Wahrscheinichkeit falsch sein. Diese Wahrscheinlichkeit kann aber durch hinreichend häufige Wiederholung beliebig gering gering gehalten werden.

Der Test fällt in die Klasse der Monte-Carlo-Algorithmen d. h. für den Test wird ein Zufallsexperiment benötigt.

Er basiert ähnlich wie der Miller-Rabin-Test auf einem Satz, mit dem festgestellt werden kann, ob eine Zahl eine Primzahl ist, oder nicht.

Vorgehensweise

Zu einer ungeraden Zahl n, von der wir feststellen wollen, ob sie eine Primzahl ist, wird zufällig eine natürliche Zahl a mit 1 < a < n gewählt.
Es wird

g := ggT(a,n)

berechnet. Falls g > 1 ist, so ist g ein echter Teiler von n, und somit ist n keine Primzahl (und der Test ist dann hier beendet).

Wenn diese Hürde überwunden ist, wird

b:=a^{\frac {n-1}{2}}\mod n

berechnet. Falls 1 < b < n-1, kann n nach einem Satz von Euler (eine etwas strengere Variante des kleinen Fermatschen Satzes) keine Primzahl sein (und der Test ist dann hier beendet). Ist andererseits b = 1 oder b = n-1, kann es sich bei n um eine Eulersche Pseudoprimzahl handeln, und der folgende Schrit muss noch ausgeführt werden.

Es wird das Jacobi-Symbol

j := J(a,n)

berechnet. Falls n eine Primzahl ist, so muss hier

j\equiv b{\pmod {n}}

gelten. Gilt dies nicht, kann also n keine Primzahl sein (und der Test ist dann hier beendet).

Nur für den Fall, dass die Berechnungen nacheinander

g = 1

b = 1 oder b = -1

j\equiv b{\pmod {n}}

ergeben, liefert der Test, dass n eine Primzahl ist.

Bewertung

Liefert der Solovay-Strassen-Test, dass n keine Primzahl ist, so ist diese Aussage korrekt.
Liefert der Test als Antwort, dass n eine Primzahl ist, so ist dies jedoch keine Garantie dafür.

Um die Güte des Solovay-Strassen-Tests zu bewerten, muss unterschieden werden, ob n eine Primzahl ist oder nicht.

Wenn n eine Primzahl ist, liefert der Test immer das korrekte Ergebnis.
Wenn n keine Primzahl ist, ist die Wahrscheinlichket, im ersten Schritt des Tests ein a zu wählen, sodass der Test ein falsches Ergebnis liefert, kleiner als 1/2 (siehe unten: Falsche Zeugen).

Um die Güte des Tests für Nichtprimzahlen zu erhöhen, wird der Test mit unabhängig gewählten zufälligen a-Werten hinreichend oft wiederholt. Wenn der Test k mal wiederholt wird, dann ist die Wahrscheinlichkeit dass in allen k Wiederholungen das Ergebnis Primzahl lautet (obwohl n keine Primzahl ist), kleiner als $1/2^{k}$ . Dies ist eine pessimistische Schätzung - in den meisten Fällen wird die Güte wesentlich besser sein.

Effizienz

Der Solovay-Strassen-Test ist effizient, da der ggT, die Potenzen und das Jacobi-Symbol effizient berechnet werden können.

Beispiel

Am Beispiel der zusammengesetzten Zahl n = 91 (einer Pseudoprimzahl) werden die möglichen Abläufe des Tests gezeigt:

Ist 91 eine Primzahl?

Fall 1: a = 7

Fall 2: a=23

Fall 3: a= 17

Fall 4: a=29

g = ggT(7,91) = 7
=>
keine Primzahl

g = ggT(23,91) = 1
b =  $23^{(91-1)/2}$  mod 91 = 64
=>
keine Primzahl

g = ggT(17,91) = 1
b =  $17^{(91-1)/2}$  mod 91 = -1
j = J(17,91) = 1
j = b
=>
keine Primzahl

g = ggT(29,91) = 1
b =  $29^{(91-1)/2}$  mod 91 = 1
j = J(29,91) = 1
j = b
=>
91 ist eine Primzahl

Die Wahrscheinlichkeit an:

Fall 1 zu scheitern beträgt bei 91 ca. 20%
Fall 2 zu scheitern beträgt bei 91 ca. 61%
Fall 3 zu scheitern beträgt bei 91 ca. 9% (geschätzt)

Falsche Zeugen

Sei n > 2 eine ungerade Nichtprimzahl. Eine Zahl a mit ggT(a,n) = 1 heißt falscher Zeuge für die Primalität von $n$ bezüglich des Solovay-Strassen-Tests, falls

$a^{\frac {n-1}{2}}\equiv J(a,n){\pmod {n}}$

Die Menge der falschen Zeugen bildet eine Untergruppe zur multiplikativen Gruppe $(\mathbb {Z} /n)^{*}$ mit Ordnung $\leq {\frac {1}{2}}\phi (n)$ ( $\phi$ bezeichnet die Eulersche φ-Funktion).

Da $\phi (n)<n$ gilt, sind höchstens die Hälfte aller zur Auswahl stehenden Zahlen a falsche Zeugen.