# Faktorisierungsverfahren

Das Faktorisierungsproblem für ganze Zahlen ist eine Aufgabenstellung aus dem mathematischen Teilgebiet der Zahlentheorie. Dabei soll zu einer zusammengesetzten Zahl ein nichttrivialer Teiler ermittelt werden. Ist beispielsweise die Zahl 91 gegeben, so sucht man eine Zahl wie 7, die 91 teilt. Entsprechende Algorithmen, die dies bewerkstelligen, bezeichnet man als Faktorisierungsverfahren. Durch rekursive Anwendung von Faktorisierungsverfahren in Kombination mit Primzahltests kann die Primfaktorzerlegung einer ganzen Zahl berechnet werden.
Bis heute ist kein Faktorisierungsverfahren bekannt, das nichttriviale Teiler und damit die Primfaktorzerlegung einer Zahl effizient berechnet. Das bedeutet, dass ein enormer Rechenaufwand notwendig ist, um eine Zahl mit mehreren hundert Stellen zu faktorisieren. Diese Schwierigkeit wird in der Kryptografie ausgenutzt. Die Sicherheit von Verschlüsselungsverfahren wie dem RSA-Kryptosystem beruht darauf, dass die Faktorisierung des RSA-Moduls zum Entschlüsseln der Nachrichten schwierig ist; somit würde ein effizientes Faktorisierungsverfahren zum Brechen des RSA-Verfahrens führen. Es ist jedoch denkbar, dass man das RSA-Problem effizienter als das Faktorisierungsproblem lösen kann. Jedoch ist bisher kein solches Verfahren bekannt.
In der theoretischen Informatik werden Probleme in Komplexitätsklassen eingeteilt, die darüber Aufschluss geben, welchen Aufwand die Lösung eines Problems erfordert. Beim Faktorisierungsproblem für ganze Zahlen ist nicht bekannt, welcher Komplexitätsklasse es angehört: Zwar ist bekannt, dass das Problem (in seiner Entscheidungsvariante) in der Klasse NP liegt, aber unbekannt, ob es bereits in polynomieller Zeit lösbar ist. Das heißt, es ist nach aktuellem Wissensstand nicht auszuschließen, dass irgendwann ein Algorithmus entdeckt wird, der ganze Zahlen mit überschaubarem Aufwand faktorisieren kann.
Die besten bekannten Algorithmen sind das 1981 von Carl Pomerance erfundene Quadratische Sieb, das um 1990 von mehreren Mathematikern (u. a. John M. Pollard, Arjen Lenstra, Hendrik Lenstra Jr., Mark S. Manasse, Carl Pomerance) gemeinsam entwickelte Zahlkörpersieb und die Methode der elliptischen Kurven, die 1987 von Hendrik W. Lenstra, Jr. vorgestellt wurde.
Die RSA Factoring Challenge verfolgte bis zu ihrer Aussetzung im Jahre 2007 den aktuellen Forschungsstand auf dem Gebiet der Faktorisierungsverfahren. Daraus ergaben sich Anhaltspunkte für die notwendige Größe der im RSA-Kryptosystem verwandten Semiprimzahlen.
In der Praxis wird man, um eine Zahl zu faktorisieren, wie folgt vorgehen:
* Durch Probedivision kleine Faktoren finden/entfernen.
* Mit Hilfe eines Primzahltests herausfinden, ob die Zahl eine Primzahl oder eine Primpotenz ist.
* Mit der Methode der elliptischen Kurven nach vergleichsweise kleinen Primfaktoren (<10$^{30}$) suchen.
* Mit dem Quadratischen Sieb (für Zahlen mit weniger als 120 Dezimalstellen) oder dem Zahlkörpersieb faktorisieren.
Die ersten beiden Schritte werden dabei gelegentlich vertauscht.

## Überblick der Faktorisierungsverfahren

Im Folgenden bezeichnet $n$ immer eine zusammengesetzte Zahl, für die ein Teiler ermittelt werden soll.

### Probedivision

Das einfachste Verfahren zur Ermittlung eines Teilers von $n$ ist die Probedivision. Dabei wird $n$ durch alle Primzahlen beginnend mit der Zwei dividiert, bis sich eine Primzahl als deren Teiler erweist oder bis der Probedivisor größer als ${\sqrt {n}}$ geworden ist. Das Verfahren eignet sich sehr gut zur Bestimmung kleiner Primfaktoren, aber es ist sehr aufwändig, damit eine Zahl mit zwei oder mehr großen Primfaktoren vollständig zu zerlegen.

### Berechnung des größten gemeinsamen Teilers

Die Probedivision kann durch den Euklidischen Algorithmus oder andere Verfahren zur Bestimmung des größten gemeinsamen Teilers so erweitert werden, dass man alle Primfaktoren von $n$ aus einem bestimmten Intervall findet. Dazu verwendet man das Produkt $m$ aller Primzahlen des Intervalls und berechnet den größten gemeinsamen Teiler der beiden Zahlen $m$ und $n$. Dieser ist das Produkt von Primfaktoren, die aus dem gewählten Intervall stammen, und man kann aus ihm die einzelnen Primfaktoren zurückgewinnen. Der Vorteil dieses Verfahrens liegt darin, dass man die Probedivision dann nur noch auf den Quotienten $n/\operatorname {ggT} (n,m)$ anwenden muss, der viel kleiner als $n$ ist.

### Faktorisierungsmethode von Fermat

Ein Verfahren, das sich besonders gut eignet, um Teiler in der Nähe von ${\sqrt {n}}$ zu finden, ist die Faktorisierungsmethode von Fermat. Dieser Algorithmus funktioniert nur für ungerade $n$ und nutzt aus, dass sich diese als Differenzen zweier Quadratzahlen darstellen lassen. Er berechnet zuerst die kleinste ganze Zahl $x$, die größer oder gleich ${\sqrt {n}}$ ist. Anschließend berechnet der Algorithmus die Differenzen $x^{2}-n$, $(x+1)^{2}-n$, $(x+2)^{2}-n$ …, bis eine dieser Differenzen eine Quadratzahl ist. Aus dieser werden Teiler von $n$ berechnet.

### Weitere Verfahren

* Faktorisierungsmethode von Lehman
* Kettenbruchmethode
* Pollard-p−1-Methode
* Pollard-p+1-Methode
* Pollard-Rho-Methode
* Dixons Faktorisierungsmethode (random squares method)
* Methode der Klassengruppen von Daniel Shanks, mit Varianten von Martin Seysen (class group relations method) und Lenstra/Schnorr (random class groups method)
* SQUFOF, square form factorization von Shanks
* Quadratisches Sieb
* Zahlkörpersieb
* Methode der elliptischen Kurven
* Schnorrs Gitter-basiertes Faktorisierungsverfahren

### Shor-Algorithmus

Eine besondere Stellung unter den Faktorisierungsverfahren nimmt der Shor-Algorithmus ein. Er kann nicht auf klassischen Rechnern ausgeführt werden, sondern benötigt einen Quantencomputer. Auf diesem kann er jedoch in Polynomialzeit einen Faktor von $n$ berechnen. Allerdings können noch keine Quantencomputer gebaut werden, die über eine für die Faktorisierung großer Zahlen ausreichende Registergröße verfügen. Die Funktion des Shor-Algorithmus beruht darauf, die Ordnung eines Elements der primen Restklassengruppe $(\mathbb {Z} /n\mathbb {Z} )^{\times }$ mit Hilfe der Quanten-Fouriertransformation zu bestimmen.
Nach Bekanntwerden des Shor-Algorithmus entwickelten Physiker technische Systeme und Versuchsanordnungen, die auf klassischem Weg, ohne Überlagerung von Quantenzuständen, die Faktorisierung natürlicher Zahlen ermöglichen. Dazu gehören z. B. Kernspinresonanz, kalte Atome, ultrakurze Lichtpulse und Mehrweg-Interferometrie.

## Geschichte
### Faktorisierungsverfahren der Antike

Seit Euklid von Alexandria ca. 300 Jahre vor Christus in seinem Hauptwerk, den Elementen, den Fundamentalsatz der Arithmetik formuliert und bewiesen hatte, war bekannt, dass jede natürliche Zahl eine eindeutige Primfaktorzerlegung besitzt. Mit der Methode der Probedivision, die im Wesentlichen ebenfalls schon Euklid bekannt war, hatte man schon sehr früh ein Verfahren gefunden, diese zu bestimmen; wenngleich es für größere Zahlen ungeeignet ist, da dann zu viel Zeit benötigt wird.

### 17. bis 19. Jahrhundert

Im Jahre 1643 beschrieb Pierre de Fermat in einem Brief (der Adressat ist nicht bekannt, vermutlich Marin Mersenne oder Bernard Frénicle de Bessy) die heutzutage nach ihm benannte Faktorisierungsmethode von Fermat, die darauf basiert, dass man die zu faktorisierende Zahl als Differenz zweier Quadrate darstellt. Diese Methode, die vom Zeitaufwand eher schlechter als die Probedivision ist, bildet die Grundlage für nahezu alle modernen Faktorisierungsverfahren.

### 20. Jahrhundert, vor der Einführung von Computern

1926 veröffentlichte Maurice Kraitchik eine Arbeit, in der er einige Verbesserungen der Faktorisierungsmethode von Fermat vorschlägt. Insbesondere betrachtet er neben der zu faktorisierenden Zahl n auch deren Vielfache, anders ausgedrückt, er sucht eine Kongruenz der Form x$^2$ ≡ y$^2$ (mod n). Um diese zu finden, multipliziert er geeignete Kongruenzen der Form x$^2$ ≡ y (mod n), die sich leicht und effektiv finden lassen (beschrieben im Artikel Quadratisches Sieb).
Derrick Henry Lehmer und Ralph Ernest Powers schlugen 1931 die sogenannte Kettenbruchmethode vor, um Kongruenzen der Form x$^2$ ≡ y (mod n) zu finden.

### 20. Jahrhundert, nach Einführung von Computern

Mit der Einführung von Computern begann die intensive Erforschung von Faktorisierungsverfahren. Aufbauend auf der Idee von Lehmer und Powers entwickelte John Brillhart in den 1960er Jahren ein auf linearer Algebra über dem endlichen Körper F$_2$ basierendes Verfahren, um aus einer Liste von Kongruenzen der Form x$^2$ ≡ y (mod n) geeignete auswählen zu können. Zusammen mit Michael Morrison gelang ihm damit im Jahre 1975 die Faktorisierung der für damalige Zeit extrem großen 39-stelligen Fermat-Zahl F$_7$. Insbesondere war es damit zum ersten Mal gelungen, ein Faktorisierungsverfahren mit subexponentieller Laufzeit (der Stellenanzahl) zu finden.
Inspiriert durch das lineare Sieb von Richard Schroeppel konnte Carl Pomerance 1981 eine Beschleunigung des Verfahrens erreichen, indem er ein Siebverfahren an Stelle der bis dato benutzten Probedivision einführte. Da das Siebverfahren sich nicht für die Kettenbruchmethode eignete, ging Pomerance wieder zu dem ursprünglich von Kraitchik vorgeschlagenen Verfahren über. Hierdurch war es möglich geworden, Zahlen mit bis zu 100 Stellen zu faktorisieren; insbesondere gelang es damit 1994, die 129-stellige Zahl RSA-129 zu zerlegen. Dieses als Quadratisches Sieb bezeichnete Verfahren ist heute noch das schnellste bekannte Verfahren zur Faktorisierung von Zahlen mit weniger als 100 Stellen.
In den 1980er Jahren vermutete man, dass Methoden, die auf der Idee von Kraitchik basieren, nicht substanziell schneller als das quadratische Sieb sein können. Diese Vermutung wurde dadurch gestützt, dass es mittlerweile etliche Verfahren mit ähnlichen Laufzeiten gab, und durch ein Ergebnis aus der analytischen Zahlentheorie über glatte Zahlen.
Anfang der 1990er Jahre wurde diese Vermutung durch das Zahlkörpersieb eindrucksvoll widerlegt. Das Zahlkörpersieb wurde 1988 von John Pollard für spezielle Zahlen vorgeschlagen und danach von einer ganzen Reihe von Mathematikern (u. a. John Pollard, Arjen Lenstra, Hendrik Lenstra, Jr., Mark Manasse, Carl Pomerance, Joe Buhler, Len Adlemann) so verändert, dass es für beliebige Zahlen anwendbar wurde. Durch den Übergang zu algebraischen Zahlkörpern war es möglich geworden, die während der Rechnung benutzten Zahlen deutlich kleiner zu halten und damit die erwähnte Beschleunigung zu erreichen. Insbesondere gelang damit 1990 die vollständige Faktorisierung der 155-stelligen Fermat-Zahl F$_9$.

### 21. Jahrhundert

Mit dem Gittersieb (einer von Pollard vorgeschlagenen Variante des Zahlkörpersiebs) und anderen Methoden wurde 2005 die Faktorisierung der bislang größten aus zwei großen Primfaktoren zusammengesetzten Zahl (einer sogenannten Fastprimzahl) ohne spezielle Struktur nach zweijähriger Arbeit auf einem Rechnerpool fertiggestellt. Dabei handelt es sich um die Zahl RSA-200, eine 200-stellige Dezimalzahl, die gemeinsam mit vielen anderen Semiprimzahlen im Rahmen der RSA Factoring Challenge generiert wurde.
Im Jahr 2012 faktorisierte eine Gruppe von 500 Teilnehmern am BOINC-Projekt NFS@Home eine Zahl mit 211 Dezimalziffern und entschlüsselte damit eine Geheimbotschaft, die im Jahr 1997 von Donald Knuth in dem Buch The Art of Computer Programming als seinerzeit unlösbare Aufgabe gestellt wurde. Knuth ersetzte daraufhin das Problem unter Verwendung einer Semiprimzahl mit 318 Dezimalziffern.
Die Liste von Champions im nach Allan Joseph Champneys Cunningham benannten Cunningham-Projekt listet aktuelle Faktorisierungsrekorde für verschiedene Zerlegungsverfahren auf.

## Implementierungen

Das Programm ARIBAS von Otto Forster implementiert verschiedene der hier besprochenen Verfahren – sei es als Bestandteil der Laufzeitbibliothek oder in Ergänzung zum Buch des Autors über Algorithmische Zahlentheorie.