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.

Eine Gitterbasis lässt sich als Matrix B = [ b 1 , … , b n ] ∈ R m {\displaystyle B=[b_{1},\dots ,b_{n}]\in \mathbb {R} ^{m}} angeben. Der kürzeste Gittervektor minimiert die Norm des Vektors ∑ i = 1 n t i b i {\displaystyle \textstyle \sum _{i=1}^{n}t_{i}b_{i}} , wobei t 1 , … , t n ∈ Z {\displaystyle t_{1},\dots ,t_{n}\in \mathbb {Z} } nicht alle gleich 0 sind. Der LLL-Algorithmus führt diese Minimierung bezüglich der euklidischen Norm durch.

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

1. 2 | r i , j | ≤ r j , j , 1 ≤ j < i ≤ n {\displaystyle 2|r_{i,j}|\leq r_{j,j},1\leq j<i\leq n} (Längenreduktion),
2. δ r i , i 2 ≤ r i + 1 , i 2 + r i + 1 , i + 1 2 , 1 ≤ i < n {\displaystyle \delta r_{i,i}^{2}\leq r_{i+1,i}^{2}+r_{i+1,i+1}^{2},1\leq i<n} (LLL-Eigenschaft)

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

Der Reduktionsfaktor δ {\displaystyle \delta } 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 δ = 1 {\displaystyle \delta =1} ist nicht bekannt, ob es einen effizienten Algorithmus gibt.

LLL-reduzierte Basen approximieren die kürzesten Vektoren mit einem Approximationsfaktor, der exponentiell in der Anzahl der Vektoren ist. Im Folgenden sei α := 1 δ − 1 4 {\displaystyle \alpha :={\tfrac {1}{\delta -{\frac {1}{4}}}}} . Offensichtlich ist α ≥ 4 3 {\displaystyle \alpha \geq {\tfrac {4}{3}}} . Der i {\displaystyle i} -te LLL-reduzierte Vektor b i {\displaystyle b_{i}} approximiert das i {\displaystyle i} -te sukzessive Minimum λ i {\displaystyle \lambda _{i}} auf folgende Weise: | | b i | | 2 λ i 2 ≤ α n − 1 {\displaystyle {\tfrac {||b_{i}||^{2}}{\lambda _{i}^{2}}}\leq \alpha ^{n-1}} , wobei das i {\displaystyle i} -te sukzessive Minimum die Länge des i {\displaystyle i} -ten kürzesten Vektors ist, der von keiner Menge kürzerer Vektoren linear abhängig ist.

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

Der Algorithmus folgt unmittelbar aus der Definition: Erzwinge Spalte für Spalte, dass B {\displaystyle B} 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.

Die Matrix R {\displaystyle R} 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 R {\displaystyle R} 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 k {\displaystyle k} -te Spalte von R {\displaystyle R} :

1. Für i = k − 1 , … , 1 {\displaystyle i=k-1,\dots ,1} wiederhole:
2. R i := R i − [ r i , k r k , k ] R k {\displaystyle R_{i}:=R_{i}-\left[{\frac {r_{i,k}}{r_{k,k}}}\right]R_{k}}

Dabei ist R i {\displaystyle R_{i}} der i {\displaystyle i} -te Spaltenvektor von R {\displaystyle R} und [ ⋅ ] {\displaystyle [\cdot ]} die Rundungsfunktion, die auf die nächste ganze Zahl rundet. Es ist wichtig, dass Schritt 2. mit fallendem und nicht mit steigendem i {\displaystyle i} ausgeführt wird.

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.^[1] 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.^[2]

Der Algorithmus findet Anwendung in vielen Bereichen, teilweise in modifizierter Version. Er lässt sich auch auf zum Z n {\displaystyle \mathbb {Z} ^{n}} 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^[3] 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^[4] oder sogar nahezu lineare Komplexität^[5] besitzen.

• 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.

• 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

1. ↑ 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
2. ↑ T. Kotnik, H. J. J. te Riele: Mertens conjecture revisited. ANTS 7 (Lecture Notes Computer Science 4076) (2006), S. 156–167
3. ↑ William R. Unger: Computing the character table of a finite group, Journal of Symbolic Computation 41 (2006) No. 8, 847–862
4. ↑ 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 (T)) SIAM J. Comput. 39 (2009) Nr. 3, S. 874–903
5. ↑ 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