LLL-Algorithmus

Der LLL-Algorithmus ist ein nach Arjen Lenstra, Hendrik Lenstra und László Lovász benannter, 1982 veröffentlichter Algorithmus, der für ein Gitter eine Basis aus möglichst kurzen Vektoren berechnet. Diese Vektoren sind Approximationen für die kürzesten voneinander linear unabhängigen Vektoren des Gitters. Bei seiner Entdeckung war der LLL-Algorithmus der erste effiziente Gitterreduktionsalgorithmus.

Kurze Gittervektoren

Eine Gitterbasis lässt sich als Matrix <math>B = [b_1, \dots, b_n] \in \mathbb R^m</math> angeben. Der kürzeste Gittervektor minimiert die Norm des Vektors <math>\textstyle \sum_{i=1}^n t_i b_i</math>, wobei <math>t_1, \dots, t_n \in \mathbb Z</math> nicht alle gleich 0 sind. Der LLL-Algorithmus führt diese Minimierung bezüglich der euklidischen Norm durch.

LLL-Basen

Definition

Üblicherweise werden LLL-Basen über ihre QR-Zerlegung <math>B = QR</math> definiert. Eine LLL-Basis <math>B = QR</math> zum Reduktionsfaktor <math>\delta \in \left[\frac{1}{4}, 1\right]</math> wird über folgende zwei Eigenschaften definiert:

1. <math>2|r_{i, j}| \le r_{j, j}, 1 \le j < i \le n</math> (Längenreduktion),
2. <math>\delta r_{i, i}^2 \le r_{i+1, i}^2 + r_{i+1, i+1}^2, 1 \le i < n</math> (LLL-Eigenschaft)

Dabei sind <math>r_{i, j}</math> die Einträge der Matrix <math>R</math>. Es wird angenommen, dass die QR-Zerlegung von <math>B</math> so durchgeführt wird, dass die Hauptdiagonale von <math>R</math> keine negativen Elemente enthält.

Der Reduktionsfaktor <math>\delta</math> gibt die Stärke der Reduktion an: Je näher er an 1 liegt, desto kürzer sind die Vektoren der reduzierten Basis. Für <math>\delta=1</math> ist nicht bekannt, ob es einen effizienten Algorithmus gibt.

Eigenschaften

LLL-reduzierte Basen approximieren die kürzesten Vektoren mit einem Approximationsfaktor, der exponentiell in der Anzahl der Vektoren ist. Im Folgenden sei <math>\alpha := \tfrac1{\delta-\frac14}</math>. Offensichtlich ist <math>\alpha \ge \tfrac43</math>. Der <math>i</math>-te LLL-reduzierte Vektor <math>b_i</math> approximiert das <math>i</math>-te sukzessive Minimum <math>\lambda_i</math> auf folgende Weise: <math>\tfrac{||b_i||^2}{\lambda_i^2} \le \alpha^{n-1}</math>, wobei das <math>i</math>-te sukzessive Minimum die Länge des <math>i</math>-ten kürzesten Vektors ist, der von keiner Menge kürzerer Vektoren linear abhängig ist.

Algorithmus

1. Berechne die erste Spalte von <math>R</math>.
2. Setze <math>k := 2</math> (<math>k</math> ist die Spalte, die gerade bearbeitet werden soll; die ersten <math>k-1</math> Spalten sind immer LLL-reduziert)
3. Solange <math>k \le n</math> wiederhole: Berechne die <math>k</math>-te Spalte von <math>R</math>, längenreduziere sie (Algorithmus siehe unten). Prüfe, ob die LLL-Eigenschaft für die Spalten <math>k-1</math> und <math>k</math> erfüllt ist.
4. Falls ja: Setze <math>k := k+1</math>
5. Falls nein: Vertausche Spalten <math>k-1</math> und <math>k</math>, setze <math>k := \max(k-1, 2)</math>.

Der Algorithmus folgt unmittelbar aus der Definition: Erzwinge Spalte für Spalte, dass <math>B</math> eine LLL-Basis ist. Längenreduktion lässt sich leicht herbeiführen. Lediglich die LLL-Eigenschaft ist kompliziert, weil sie indirekt weit voneinander entfernte Spalten in Beziehung bringt. Darum ist es nötig, Spalten zu vertauschen und wieder einen Schritt zurückzugehen.

Längenreduktion

Die Matrix <math>R</math> ist eine obere Dreiecksmatrix mit positiver Diagonale. Ziel der Längenreduktion ist es, alle Nicht-Diagonalelemente möglichst (betragsmäßig) klein zu machen, ohne das Gitter zu verändern. Jede Zeile von <math>R</math> hat folgenden Aufbau: null oder mehr Nullen, das Diagonalelement, null oder mehr weitere Elemente. Dabei besitzen die Nullen schon den kleinstmöglichen Betrag. Die LLL-Eigenschaft sorgt dafür, dass das Diagonalelement nicht allzu groß sein kann. Die Längenreduktion bewirkt nun für alle folgenden Elemente, dass sie betragsmäßig höchstens halb so groß wie das Diagonalelement sind. Das lässt sich dadurch herbeiführen, dass für jedes zu große Nicht-Diagonalelement die Spalte, die das entsprechende Diagonalelement enthält, so oft addiert oder subtrahiert wird, bis das zu große Element betragsmäßig minimal ist. Der folgende Algorithmus längenreduziert die <math>k</math>-te Spalte von <math>R</math>:

1. Für <math>i = k-1, \dots, 1</math> wiederhole:
2. <math>R_i := R_i - \left[\frac{r_{i,k}}{r_{k,k}}\right] R_k</math>

Dabei ist <math>R_i</math> der <math>i</math>-te Spaltenvektor von <math>R</math> und <math>[\cdot]</math> die Rundungsfunktion, die auf die nächste ganze Zahl rundet. Es ist wichtig, dass Schritt 2. mit fallendem und nicht mit steigendem <math>i</math> ausgeführt wird.

Anwendungen

Eine erste hervorragende Anwendung erhielt der LLL-Algorithmus 1984 bei der Widerlegung der mit der Riemannschen Vermutung zusammenhängenden Mertensschen Vermutung durch Andrew Odlyzko und Herman te Riele.<ref>A. M. Odlyzko, H. J. J. te Riele: Disproof of the Mertens conjecture. (Memento vom 16. Februar 2013 im Internet Archive) J. reine angew. Math. 357 (1985), S. 138–160</ref> Dabei wurde ein Gitter der Dimension 70 konstruiert und seine Basis reduziert, das es erlaubte, ein Problem der mehrdimensionalen inhomogenen diophantischen Approximation mit ausgewählten nichttrivialen Nullstellen der Riemannschen Zetafunktion näherungsweise zu lösen. 2006 wurde die Berechnung an Gittern der Dimension um 100 wiederholt und die Ergebnisse verbessert.<ref>T. Kotnik, H. J. J. te Riele: Mertens conjecture revisited. ANTS 7 (Lecture Notes Computer Science 4076) (2006), S. 156–167</ref>

Der Algorithmus findet Anwendung in vielen Bereichen, teilweise in modifizierter Version. Er lässt sich auch auf zum <math>\mathbb Z ^n</math> isometrisch isomorphen Gittern anwenden, wobei die Reduktion dann bezüglich des dortigen Skalarprodukts durchgeführt wird. Ein eher unerwartetes Beispiel hierfür ist der Algorithmus von Unger<ref>William R. Unger: Computing the character table of a finite group, Journal of Symbolic Computation 41 (2006) No. 8, 847–862</ref> aus der Gruppentheorie, der zum Finden von irreduziblen Charakteren einer Gruppe verwendet wird.

Da in bestimmten Anwendungen des Algorithmus die anfangs vorliegenden Basisvektoren eine außerordentliche Länge aufweisen können, wurde angestrebt, seine Laufzeit in Abhängigkeit von dieser Länge zu verringern. Es wurden Varianten entwickelt, die eine quadratische Zeitkomplexität bezüglich der Bitlänge der Ausgangsdaten<ref>Phong Q. Nguyen, Damien Stehlé: @1@2Vorlage:Toter Link/ftp.di.ens.frAn LLL algorithm with quadratic complexity. (Seite nicht mehr abrufbar, festgestellt im November 2022. Suche im Internet Archive ) SIAM J. Comput. 39 (2009) Nr. 3, S. 874–903</ref> oder sogar nahezu lineare Komplexität<ref>Andrew Novocin, Damien Stehlé, Gilles Villard: An LLL-reduction algorithm with quasi-linear time complexity. STOC 2011 (Symposium on Theory of Computing), S. 403–412</ref> besitzen.

Literatur

• A. K. Lenstra, H. W. Lenstra, Jr., L. Lovász: Factoring polynomials with rational coefficients. In: Mathematische Annalen. Band 261, Nr. 4, 1982, S. 515–534, doi:10.1007/BF01457454. 
• Phong Q. Nguyen, Brigitte Valée (Hrsg.): The LLL algorithm. Survey and applications. Springer 2010, ISBN 978-3-642-02294-4.
• Murray R. Bremner: Lattice basis reduction: An introduction to the LLL algorithm and its applications. CRC Press 2011, ISBN 978-1-439-80702-6.

Weblinks

• Gitter und Kryptographie (PDF; 1015 KB). Johann Wolfgang Goethe-Universität Frankfurt am Main, 17. Oktober 2011 (Skript für die Vorlesungen von Prof. C. P. Schnorr, Wintersemester 2010)
• Phong Q. Nguyen: Lattice reduction algorithms: Theory and practice. eingeladener Powerpoint-Vortrag zur EUROCRYPT'11-Konferenz

Einzelnachweise

<references />