AKS-Primzahltest

Der AKS-Primzahltest (auch bekannt unter dem Namen Agrawal-Kayal-Saxena-Primzahltest) ist ein deterministischer Algorithmus, der für eine natürliche Zahl in polynomieller Laufzeit feststellt, ob sie prim ist oder nicht. Er wurde von den drei indischen Wissenschaftlern Manindra Agrawal, Neeraj Kayal und Nitin Saxena entdeckt und 2002 in einer Abhandlung mit dem Titel PRIMES is in P (deutsch sinngemäß: Das Primzahl-Problem gehört zur Komplexitätsklasse P) veröffentlicht.

Der später von anderen verbesserte Algorithmus unterscheidet sich wesentlich von allen vorher bekannten polynomiellen Primalitätsbeweis-Algorithmen: Er baut für den Nachweis der – bezogen auf die Länge der Eingangswerte – polynomiellen Laufzeit auf keinen unbewiesenen Hypothesen (wie beispielsweise der verallgemeinerten Riemannschen Vermutung) auf. Die asymptotische Laufzeit des ursprünglichen Algorithmus ist $\mathcal{O}\left((\log n)^{12+\varepsilon}\right)$ (Landau-Symbol $\mathcal{O}$), wobei n die zu testende Zahl ist.

Entstehungsgeschichte

1999 arbeitete Manindra Agrawal mit seinem Doktorvater Somenath Biswas an einer probabilistischen Methode, um die Gleichheit von Polynomen zu zeigen. Die beiden erarbeiteten daraus eine Methode für einen probabilistischen Primzahltest. Die Idee, die dahinter steckt und die sich später als sinnvoll herausstellte, ist folgender Hilfssatz:

Sei $a\in \mathbb Z, n\in \mathbb N, N>1$ und ggT(a,n)=1. Dann ist n eine Primzahl genau dann, wenn $(x+a)^n\equiv x^n+a\,\,(\text{mod }n)$.

Dabei ist x keine konkrete Zahl, sondern eine freie Variable. Die Koeffizienten der Potenzen von x sind zu vergleichen.

Für den so entstandenen Primzahltest galt, dass er nicht mit den aktuellen mithalten konnte. Im schlimmsten Falle musste man alle Koeffizienten berechnen, was ziemlich aufwendig sein konnte. Deshalb wurde die Idee zunächst nicht weiter verfolgt.

2001 nahmen die Studenten Rajat Bhattarcharjee und Prashant Pandey in ihrer Bachelorarbeit Primality Testing die Idee wieder auf. Sie erweiterten die Idee, die Polynome nicht nur modulo n sondern auch modulo x^r − 1 für ein r in der Größenordnung von log n zu berechnen. Dies hat den Vorteil, dass man $(x+a)^n \ (\text{mod } (x^r-1,n))$ in polynomieller Zeit berechnen kann. Nun gilt für eine Primzahl n, dass sie diese Kongruenz erfüllt, aber es erfüllen sie nun auch Zahlen, die keine Primzahlen sind.

Die beiden untersuchten diese Kongruenz für bestimmte a und r, um Bedingungen an a und r zu stellen, damit diese Kongruenz nur noch für Primzahlen gilt. Sie stellten nach einer Versuchsreihe die folgende Vermutung auf:

Ist $r\in \mathbb P$ kein Teiler von n und $(x+1)^n \equiv x^n+1\, (\text{mod }(x^r-1,n))$, dann ist n entweder prim oder $n^2\equiv 1\, (\text{mod }r)$.

2002 arbeiteten die beiden Studenten Neeraj Kayal und Nitin Saxena an ihrer Bachelorarbeit. Sie führten die Überlegungen ihrer Vorreiter weiter. Unter der Annahme, dass die Riemannsche Vermutung korrekt ist, konnten sie den obigen Satz beweisen. In einer leichten Vorahnung nannten sie dann ihre Bachelorarbeit: Towards a deterministic polynomial-time Primality Test.

Danach brachten sie den Algorithmus mit Manindra Agrawal in seine endgültige Form. Die dann veröffentlichte Schrift erfreute sich ziemlich schnell einer großen Beliebtheit. So wurde die Korrektheit innerhalb einer Woche bestätigt und die Webseite hatte über zwei Millionen Besucher in der ersten Woche.

Der Algorithmus

Im Folgenden sind n,a,b,r natürliche Zahlen. Die Eingabe ist die Zahl n>1.

1. if (n ist eine reine Potenz) return ZUSAMMENGESETZT
2. finde das kleinste r mit o_r(n)>log(n)2
3. if 1<ggT(a,n)<n für ein a ≤ r return ZUSAMMENGESETZT
4. if n ≤ r return PRIM
5. for a=1 to sqrt(phi(r))*log(n) do
      if (x+a)n ≠ xn+a (mod (xr-1,n)) return ZUSAMMENGESETZT
6. return PRIM

Erläuterungen:

• o_r(n) ist die Ordnung von n modulo r, das ist die kleinste Zahl k, für die $n^k \equiv 1\,(\text{mod }r)$ gilt.
• phi(x)=φ(x) ist die Eulersche Phi-Funktion oder Totient, das ist die Anzahl der zu x teilerfremden Zahlen im Bereich von 1 bis x.

Nach Agrawal, Kayal und Saxena

In den folgenden Monaten nach der Entdeckung erschienen neue Varianten (Lenstra 2002, Pomerance 2002, Berrizbeitia 2003, Cheng 2003, Bernstein 2003a/b, Lenstra und Pomerance 2003), die die AKS-Geschwindigkeit um Größenordnungen verbesserten. Wegen der großen Anzahl an Varianten sprechen Crandall und Papadopoulos in ihrem Aufsatz On the implementation of AKS-class primality tests (Über die Implementation von Primzahltests der AKS-Klasse) von März 2003 von der Klasse der AKS-Algorithmen, statt vom AKS-Algorithmus.

Der Algorithmus von Lenstra und Pomerance terminiert in $\mathcal{O}\left((\log N)^{6+\varepsilon}\right)$. Die Laufzeit des AKS-Algorithmus bewegt sich in der Praxis jedoch in ähnlichen Größenordnungen, da der Parameter r meist wenig oberhalb von log ²n gefunden werden kann.

Agrawal, Kayal und Saxena haben mit der obigen Vermutung einen ähnlichen Algorithmus aufgestellt:

Man suche zuerst ein $r\in\mathbb P$ mit $r\nmid(N^2-1)$ (so ein r liegt im Intervall 2,4log N). Mit diesem Algorithmus erhält man eine Laufzeit von $\mathcal{O}\left((\log n)^{3+\varepsilon}\right)$. Lenstra und Pomerance haben zu dieser Vermutung in "Remarks on Agrawal's Conjecture" eine Heuristik zum Finden von möglichen Gegenbeispielen angegeben. Ob es Zahlen wie die in deren Vermutung angenommenen gibt, ist jedoch bisher nicht bekannt.

Weblinks

• MathWorld: AKS Primality Test
• PRIMES is in P: Scientific paper describing AKS primality test (PDF; 195 kB)
• R. Crandall, Apple ACG, and J. Papadopoulos (March 18, 2003): On the implementation of AKS-class primality tests (PDF; 99 kB)
• Artikel von Folkmar Bornemann über den AKS-Primzahltest mit Fotos der drei indischen Forscher (PDF; 487 kB)
• Lenstra, Pomerance - Remarks on Agrawal's Conjecture