# Kettenbruchmethode

Die Kettenbruchmethode (Abk.: CFRAC) berechnet zwei Teiler einer natürlichen Zahl, die keine Primzahl ist. Durch wiederholte Anwendung lässt sich so die Primfaktorzerlegung dieser Zahl ermitteln.
Die Kettenbruchmethode wurde 1931 von Derrick Henry Lehmer und dem Mathematik-Liebhaber Ralph Ernest Powers veröffentlicht, der auch durch zahlentheoretische Rechenergebnisse bekannt ist. Über Jahrzehnte hinweg wurde sie kaum verwendet, da sie auf den zur Verfügung stehenden Rechenmaschinen häufig nach mehreren Stunden noch keine Faktorisierung fand. Im Jahr 1970 programmierten Michael Morrison und John Brillhart die Kettenbruchmethode auf einem Großrechner und berechneten damit die Faktorisierung der siebten Fermat-Zahl. In den achtziger Jahren war die Kettenbruchmethode das Standardverfahren für die Faktorisierung großer Zahlen. Das waren damals Zahlen mit bis zu fünfzig Dezimalstellen. 1990 wurde mit dem quadratischen Sieb ein neuer Algorithmus vorgestellt, der die Nachfolge der Kettenbruchmethode antrat. Die Laufzeit der Kettenbruchmethode in O-Notation ist $O\left(e^{\sqrt {2\log N\log \log N}}\right)$, wobei N die zu faktorisierende Zahl ist.

## Funktionsweise

Der Algorithmus 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). Solche Kongruenzen führen zu einer Faktorisierung von N (genauer beschrieben im Artikel Dixons Faktorisierungsmethode.)
Bei der Kettenbruchmethode nutzt man die Kettenbruchentwicklung von ${\sqrt {N}}$, um solche Kongruenzen zu finden. Die Kettenbruchentwicklung liefert beste Näherungen der Wurzel von N als Brüche der Form ${\frac {p}{q}}$. Jede Näherung liefert eine Kongruenz x$^2$ ≡ y (mod N) für x = p und y = x$^2$ − q$^2$ · N. Da der Bruch eine beste Näherung für ${\sqrt {N}}$ ist, ist y klein und mit hoher Wahrscheinlichkeit glatt und man braucht nur wenige solcher Kongruenzen.

## Algorithmus nach Morrison und Brillhart

Die hier beschriebene Variante der Kettenbruchmethode entspricht derjenigen, die von Morrison und Brillhart in ihrem Aufsatz „A Method of Factoring and the Factorization of $F_{7}$“ veröffentlicht wurde. Der Algorithmus besteht aus drei wesentlichen Schritten. Die zu faktorisierende Zahl wird im Weiteren mit $N$ bezeichnet.

### Schritt A

Wähle eine natürliche Zahl $k\geq 1$ und berechne die Kettenbruchentwicklung von ${\sqrt {kN}}$. Die Kettenbruchentwicklung kann an einer beliebigen Stelle $n-1$ abgebrochen werden, sodass man einen Kettenbruch der Form
    $[q_{0};q_{1},q_{2},\ldots ,q_{n-1},{\frac {{\sqrt {kN}}+P_{n}}{Q_{n}}}]$
erhält. Man berechnet nun die Paare $(A_{i-1},Q_{i})$, wobei $A_{i}$ der Zähler und $B_{i}$ der Nenner der $i$-ten Konvergente sind und sich $Q_{i}$ nach der Formel
    $Q_{i}=(-1)^{i}\cdot (A_{i-1}^{2}-kNB_{i-1}^{2})$
berechnet.

### Schritt B

Aus den in Schritt A erzeugten $(A,Q)$-Paaren werden diejenigen ausgewählt, bei denen alle Primfaktoren von $Q$ einer vorher festgelegten Faktorbasis entstammen. Die Faktorbasis enthält üblicherweise die Zahl −1 und alle Primzahlen, die Teiler von Q sein können, bis zu einer bestimmten Schranke. Man braucht nur die Primzahlen p mit $p=2$ oder $(kN)^{(p-1)/2}\,{\bmod {\;}}p\leq 1$ in die Faktorbasis aufzunehmen, denn Q kann nur durch diese p teilbar sein.
Durch mehrfache Probedivision lässt sich feststellen, ob das jeweilige $Q$ ein Produkt von Zahlen aus der Faktorbasis ist, und nebenbei erhält man die Primfaktorzerlegung von $Q$.
Mit den ausgewählten $(A,Q)$-Paaren füllt man eine Tabelle. Diese Tabelle enthält eine Spalte für jede Zahl der Faktorbasis. In die Faktorbasisspalten wird die Zahl 1 eingetragen, wenn der jeweilige Faktor in der Primzahlzerlegung von $Q$ ungerade oft vorkommt. Andernfalls trägt man 0 ein. Der Tabelleneintrag für das Paar $(375,-220=-1\cdot 2^{2}\cdot 5\cdot 11)$ und die Faktorbasis $\{-1,2,3,5,7,11,13\}$ sieht beispielsweise wie folgt aus:
| | −1 | 2 | 3 | 5 | 7 | 11 | 13 |
| (375, −220) | 1 | 0 | 0 | 1 | 0 | 1 | 0 |
Dann bestimmt man eine Auswahl von Tabellenzeilen, für die die Summe der Einträge in jeder Tabellenspalte gerade ist, sodass das Produkt der $Q$ der ausgewählten Zeilen jeden Faktor mit einer geraden Potenz enthält und somit eine Quadratzahl ist. Das ist sicher dann möglich, wenn die Zahl der Tabellenzeilen mindestens um eins größer ist als die Zahl der Faktoren in der Faktorbasis und damit der Tabellenspalten.

### Schritt C

Man multipliziert die A aller Paare und die Q aller Paare der ausgewählten Tabellenzeilen. So erhält man die Legendre-Kongruenz
    $x^{2}=\prod A^{2}\equiv \prod Q=y^{2}\;({\bmod {\,}}N)$.
Man berechnet dann $d_{1}=\operatorname {ggT} (x-y,N)$ und $d_{2}=\operatorname {ggT} (x+y,N)$.
Wenn $x\not \equiv y\;({\bmod {\,}}N)$ und $x\not \equiv -y\;({\bmod {\,}}N)$, dann sind $d_{1}$ und $d_{2}$ echte Teiler von N. Anderenfalls kann eine andere Auswahl von Tabellenzeilen erfolgreich sein. Ggfs. muss man noch weitere Zeilen berechnen.

## Geschichte

Der erste Schritt auf dem Weg zur Kettenbruchmethode war die 1643 von Pierre de Fermat beschriebene Faktorisierungsmethode von Fermat. Dieses Rechenverfahren findet zwei Quadratzahlen $x^{2}$ und $y^{2}$, sodass $x^{2}-y^{2}=N$ gilt. $N$ ist auch hier wieder die zu faktorisierende Zahl. Aufgrund der 3. Binomischen Formel gilt
    $N=x^{2}-y^{2}=(x+y)(x-y)$
und $x+y$ und $x-y$ sind deshalb Teiler von $N$.
1798 veröffentlichte Adrien-Marie Legendre sein Buch „Essai sur la théorie des nombres“. In diesem findet sich mit den Legendre-Kongruenzen eine Weiterentwicklung der Methode von Fermat. Die Differenz der Quadratzahlen $x^{2}$ und $y^{2}$ muss nicht mehr gleich $N$ sein, sondern kann ein beliebiges Vielfaches von $N$ sein. Die Wurzeln $x$ und $y$ der Quadratzahlen müssen zudem drei Bedingungen erfüllen: $0<x,y<N$, $x\neq y$ und $x+y\neq N$. Unter diesen Voraussetzungen ist $N$ ein Teiler von $x^{2}-y^{2}=(x+y)(x-y)$ und sowohl der größte gemeinsame Teiler $\operatorname {ggT} (N,x+y)$ als auch $\operatorname {ggT} (N,x-y)$ sind Faktoren von $N$.