Analysis of BKZ Guillaume Hanrot, Xavier Pujol, Damien Stehl e - - PowerPoint PPT Presentation

analysis of bkz
SMART_READER_LITE
LIVE PREVIEW

Analysis of BKZ Guillaume Hanrot, Xavier Pujol, Damien Stehl e - - PowerPoint PPT Presentation

Analysis of BKZ Guillaume Hanrot, Xavier Pujol, Damien Stehl e ENSL, LIP, CNRS, INRIA, Universit e de Lyon, UCBL May 5, 2011 Analysis of BKZ 1/32 b b b b b b b b b b b b b b b b b Lattices a 1 a 2 (SVP) Analysis of BKZ


slide-1
SLIDE 1

Analysis of BKZ

Guillaume Hanrot, Xavier Pujol, Damien Stehl´ e

ENSL, LIP, CNRS, INRIA, Universit´ e de Lyon, UCBL

May 5, 2011

Analysis of BKZ 1/32

slide-2
SLIDE 2

Lattices

a1 a2

b b b b b b b b b b b b b b b b b

(SVP)

Analysis of BKZ 2/32

slide-3
SLIDE 3

Lattices

b1

b b b b b b b b b b b b b b b b b

(SVP)Shortest vector problem (SVP)(SVP)

Analysis of BKZ 2/32

slide-4
SLIDE 4

Lattices

b1 b2

b b b b b b b b b b b b b b b b b

(SVP)Lattice reduction(SVP)

Analysis of BKZ 2/32

slide-5
SLIDE 5

Lattices

b1 b2

b b b b b b b b b b b b b b b b b

(SVP)Determinant(SVP)

Analysis of BKZ 2/32

slide-6
SLIDE 6

Lattices

b1 b2

b b b b b b b b b b b b b b b b b

Hermite factor: HF(b1, . . . , bn) = b1 (det L)1/n Goal of lattice reduction: find a basis with small HF. If b1 is a shortest vector, then HF(b1, . . . , bn) ≤ √γn, with γn = Hermite constant ≤ n.

Analysis of BKZ 2/32

slide-7
SLIDE 7

Lattices

b1 b2

b b b b b b b b b b b b b b b b b

Hermite factor: HF(b1, . . . , bn) = b1 (det L)1/n Goal of lattice reduction: find a basis with small HF. If b1 is a shortest vector, then HF(b1, . . . , bn) ≤ √γn, with γn = Hermite constant ≤ n.

Analysis of BKZ 2/32

slide-8
SLIDE 8

Lattices

b1 b2

b b b b b b b b b b b b b b b b b

Hermite factor: HF(b1, . . . , bn) = b1 (det L)1/n Goal of lattice reduction: find a basis with small HF. If b1 is a shortest vector, then HF(b1, . . . , bn) ≤ √γn, with γn = Hermite constant ≤ n.

Analysis of BKZ 2/32

slide-9
SLIDE 9

Lattice reduction

Lattice reduction and shortest vector problem: The security of lattice-based cryptosystems relies on the hardness of (variants of) SVP. SVP and lattice reduction are interdependent problems. Hierarchy of lattice reductions in dimension n: HKZ BKZβ LLL Hermite factor √γn ≃ (γβ(1 + ǫ))

n−1 2(β−1)

(γ2(1 + ǫ))

n−1 2

Time 2O(n) 2O(β)×? Poly(n)

HKZ = Hermite-Korkine-Zolotareff BKZ = Block Korkine-Zolotareff

Analysis of BKZ 3/32

slide-10
SLIDE 10

Lattice reduction

Lattice reduction and shortest vector problem: The security of lattice-based cryptosystems relies on the hardness of (variants of) SVP. SVP and lattice reduction are interdependent problems. Hierarchy of lattice reductions in dimension n: HKZ BKZβ LLL Hermite factor √γn ≃ (γβ(1 + ǫ))

n−1 2(β−1)

(γ2(1 + ǫ))

n−1 2

Time 2O(n) 2O(β)×? Poly(n)

HKZ = Hermite-Korkine-Zolotareff BKZ = Block Korkine-Zolotareff

Analysis of BKZ 3/32

slide-11
SLIDE 11

Lattice reduction

Lattice reduction and shortest vector problem: The security of lattice-based cryptosystems relies on the hardness of (variants of) SVP. SVP and lattice reduction are interdependent problems. Hierarchy of lattice reductions in dimension n: HKZ BKZβ LLL Hermite factor √γn ≃ (γβ(1 + ǫ))

n−1 2(β−1)

(γ2(1 + ǫ))

n−1 2

Time 2O(n) 2O(β)×? Poly(n)

HKZ = Hermite-Korkine-Zolotareff BKZ = Block Korkine-Zolotareff

Analysis of BKZ 3/32

slide-12
SLIDE 12

Lattice reduction

Lattice reduction and shortest vector problem: The security of lattice-based cryptosystems relies on the hardness of (variants of) SVP. SVP and lattice reduction are interdependent problems. Hierarchy of lattice reductions in dimension n: HKZ BKZβ LLL Hermite factor √γn ≃ (γβ(1 + ǫ))

n−1 2(β−1)

(γ2(1 + ǫ))

n−1 2

Time 2O(n) 2O(β)×? Poly(n)

HKZ = Hermite-Korkine-Zolotareff BKZ = Block Korkine-Zolotareff

Analysis of BKZ 3/32

slide-13
SLIDE 13

History of BKZ

Practice Theory Schnorr and Euchner (1994): algorithm for BKZ-reduction, without complexity analysis. Shoup: first public implementation of BKZ in NTL. Gama and Nguyen (2008): BKZ behaves badly when the block size is ≥ 25. Schnorr (1987): first hierarchies of algorithms between LLL and HKZ. Gama et al. (2006): Block-Rankin-reduction. Gama and Nguyen (2008): Slide-reduction.

Analysis of BKZ 4/32

slide-14
SLIDE 14

History of BKZ

Practice Theory Schnorr and Euchner (1994): algorithm for BKZ-reduction, without complexity analysis. Shoup: first public implementation of BKZ in NTL. Gama and Nguyen (2008): BKZ behaves badly when the block size is ≥ 25. Schnorr (1987): first hierarchies of algorithms between LLL and HKZ. Gama et al. (2006): Block-Rankin-reduction. Gama and Nguyen (2008): Slide-reduction.

Analysis of BKZ 4/32

slide-15
SLIDE 15

History of BKZ

Practice Theory Schnorr and Euchner (1994): algorithm for BKZ-reduction, without complexity analysis. Shoup: first public implementation of BKZ in NTL. Gama and Nguyen (2008): BKZ behaves badly when the block size is ≥ 25. Schnorr (1987): first hierarchies of algorithms between LLL and HKZ. Gama et al. (2006): Block-Rankin-reduction. Gama and Nguyen (2008): Slide-reduction.

Analysis of BKZ 4/32

slide-16
SLIDE 16

History of BKZ

Practice Theory Schnorr and Euchner (1994): algorithm for BKZ-reduction, without complexity analysis. Shoup: first public implementation of BKZ in NTL. Gama and Nguyen (2008): BKZ behaves badly when the block size is ≥ 25. Schnorr (1987): first hierarchies of algorithms between LLL and HKZ. Gama et al. (2006): Block-Rankin-reduction. Gama and Nguyen (2008): Slide-reduction.

Analysis of BKZ 4/32

slide-17
SLIDE 17

History of BKZ

Practice Theory Schnorr and Euchner (1994): algorithm for BKZ-reduction, without complexity analysis. Shoup: first public implementation of BKZ in NTL. Gama and Nguyen (2008): BKZ behaves badly when the block size is ≥ 25. Schnorr (1987): first hierarchies of algorithms between LLL and HKZ. Gama et al. (2006): Block-Rankin-reduction. Gama and Nguyen (2008): Slide-reduction.

Analysis of BKZ 4/32

slide-18
SLIDE 18

History of BKZ

Practice Theory Schnorr and Euchner (1994): algorithm for BKZ-reduction, without complexity analysis. Shoup: first public implementation of BKZ in NTL. Gama and Nguyen (2008): BKZ behaves badly when the block size is ≥ 25. Schnorr (1987): first hierarchies of algorithms between LLL and HKZ. Gama et al. (2006): Block-Rankin-reduction. Gama and Nguyen (2008): Slide-reduction.

Analysis of BKZ 4/32

slide-19
SLIDE 19

Slide-reduction: Outputs a basis whose theoretical quality is equivalent to BKZ. Polynomial number of calls to a SVP oracle. Not as efficient as BKZ in practice.

Analysis of BKZ 5/32

slide-20
SLIDE 20

Progress made during the execution of BKZ

1.012 1.013 1.014 1.015 1.016 1.017 1.018 1.019 1.02 1.021 20 40 60 80 100 Hermite factor Number of tours Quality of BKZ output BKZ BKZ’

Experience on 64 LLL-reduced knapsack-like matrices (n = 108, β = 24).

Analysis of BKZ 6/32

slide-21
SLIDE 21

Progress made during the execution of BKZ

1.012 1.013 1.014 1.015 1.016 1.017 1.018 1.019 1.02 1.021 200 400 600 800 1000 1200 Hermite factor Number of tours Quality of BKZ output BKZ BKZ’

Experience on 64 LLL-reduced knapsack-like matrices (n = 108, β = 24).

Analysis of BKZ 6/32

slide-22
SLIDE 22

Our result

γβ = Hermite constant ≤ β. L a lattice with basis (b1, . . . , bn).

Theorem

After O n3 β2

  • log n

ǫ + log log max bi (det L)1/n

  • calls to HKZβ,

BKZβ returns a basis C of L such that: HF(C) ≤ (1 + ǫ)γβ

n−1 2(β−1)+ 3 2 Analysis of BKZ 7/32

slide-23
SLIDE 23

1

Reminders on lattice reduction

2

Analysis of BKZ in the sandpile model

3

Analysis of BKZ

4

Applications to LLL

5

Conclusion

Analysis of BKZ 8/32

slide-24
SLIDE 24

1

Reminders on lattice reduction

2

Analysis of BKZ in the sandpile model

3

Analysis of BKZ

4

Applications to LLL

5

Conclusion

Analysis of BKZ 9/32

slide-25
SLIDE 25

Gram-Schmidt orthogonalization

b1, . . . , bn linearly independent. The Gram-Schmidt orthogona- lization b∗

1, . . . , b∗ n is defined by:

For all i > j, µi,j =

(bi,b∗

j )

b∗

j 2 .

For all i, b∗

i = bi − j<i µi,jb∗ j .

A basis is size-reduced if all the |µi,j| are ≤ 1

2. b2 b3 b1

Analysis of BKZ 10/32

slide-26
SLIDE 26

Gram-Schmidt orthogonalization

b1, . . . , bn linearly independent. The Gram-Schmidt orthogona- lization b∗

1, . . . , b∗ n is defined by:

For all i > j, µi,j =

(bi,b∗

j )

b∗

j 2 .

For all i, b∗

i = bi − j<i µi,jb∗ j .

A basis is size-reduced if all the |µi,j| are ≤ 1

2. = b∗

2

b∗

3

b2 b3 b1 b∗

1

Analysis of BKZ 10/32

slide-27
SLIDE 27

Gram-Schmidt orthogonalization

b1, . . . , bn linearly independent. The Gram-Schmidt orthogona- lization b∗

1, . . . , b∗ n is defined by:

For all i > j, µi,j =

(bi,b∗

j )

b∗

j 2 .

For all i, b∗

i = bi − j<i µi,jb∗ j .

A basis is size-reduced if all the |µi,j| are ≤ 1

2. = b∗

2

b∗

3

b2 b3 b1 b∗

1

Analysis of BKZ 10/32

slide-28
SLIDE 28

LLL

B is δ-LLL-reduced if: It is size-reduced; δb∗

i 2 ≤ b∗ i+12 + µ2 i+1,ib∗ i 2 for all i < n.

→ xi ≤ 1

2 log γ2 + xi+1 − log δ

(xi = log b∗

i )

Analysis of BKZ 11/32

slide-29
SLIDE 29

LLL

B is δ-LLL-reduced if: It is size-reduced; δb∗

i 2 ≤ b∗ i+12 + µ2 i+1,ib∗ i 2 for all i < n.

→ xi ≤ 1

2 log γ2 + xi+1 − log δ

(xi = log b∗

i )

x5 x4 x3

1 2 log γ2 − log δ

x2 x1

Analysis of BKZ 11/32

slide-30
SLIDE 30

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

b b

b3 b2 b1

b

Analysis of BKZ 12/32

slide-31
SLIDE 31

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

b3 b2 b1

b b b b b b b b b b b b b

Analysis of BKZ 12/32

slide-32
SLIDE 32

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

b b

b3 b2 b1

b

Analysis of BKZ 12/32

slide-33
SLIDE 33

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

b(2)

3

b(2)

2

= b∗

2

b

Analysis of BKZ 12/32

slide-34
SLIDE 34

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

b(2)

3

b(2)

2

= b∗

2

b b b b b b b b b b

Analysis of BKZ 12/32

slide-35
SLIDE 35

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

b(2)

3

b(2)

2

= b∗

2

b

Analysis of BKZ 12/32

slide-36
SLIDE 36

HKZ

B is HKZ-reduced if: It is size-reduced. b∗

i = shortest vector of L(b(i) i , . . . , b(i) n ).

For i < n, HF(b(i)

i , . . . , b(i) n ) ≤ √γn−i+1

Worst-case HKZ profile: xi = log b∗

i

= O(log2(n − i)) x10 x9 x8 x7 x6 x5 x4 x3 x2 x1

Analysis of BKZ 12/32

slide-37
SLIDE 37

BKZ

Algorithm (BKZβ, modified version)

Input: B of dimension n. Repeat ... times For i from 1 to n − β + 1 do Size-reduce B. HKZ-reduce the projected sublattice (b(i)

i , . . . , b(i) i+β−1).

Report the transformation on B. Termination?

Analysis of BKZ 13/32

slide-38
SLIDE 38

BKZ

Algorithm (BKZβ, modified version)

Input: B of dimension n. Repeat ... times For i from 1 to n − β + 1 do Size-reduce B. HKZ-reduce the projected sublattice (b(i)

i , . . . , b(i) i+β−1).

Report the transformation on B. Termination?

Analysis of BKZ 13/32

slide-39
SLIDE 39

1

Reminders on lattice reduction

2

Analysis of BKZ in the sandpile model

3

Analysis of BKZ

4

Applications to LLL

5

Conclusion

Analysis of BKZ 14/32

slide-40
SLIDE 40

Sandpile model

We consider only xi = log b∗

i for i ≤ n.

Each HKZ-reduction gives a worst-case profile. → The initial xi’s fully determine the xi’s after a call to HKZ. The sandpile execution of BKZ is deterministic.

Analysis of BKZ 15/32

slide-41
SLIDE 41

Sandpile model

We consider only xi = log b∗

i for i ≤ n.

Each HKZ-reduction gives a worst-case profile. → The initial xi’s fully determine the xi’s after a call to HKZ. The sandpile execution of BKZ is deterministic.

Analysis of BKZ 15/32

slide-42
SLIDE 42

Sandpile model

We consider only xi = log b∗

i for i ≤ n.

Each HKZ-reduction gives a worst-case profile. → The initial xi’s fully determine the xi’s after a call to HKZ. The sandpile execution of BKZ is deterministic.

Analysis of BKZ 15/32

slide-43
SLIDE 43

x1 x2 x3 x4 x5 x6 x7 x8 x9

Analysis of BKZ 16/32

slide-44
SLIDE 44

x5 x6 x7 x8 x9 x1 x2 x3 x4

Analysis of BKZ 16/32

slide-45
SLIDE 45

x1 x6 x7 x8 x9 x2 x3 x4 x5

Analysis of BKZ 16/32

slide-46
SLIDE 46

x1 x2 x7 x8 x9 x3 x4 x5 x6

Analysis of BKZ 16/32

slide-47
SLIDE 47

x1 x2 x3 x8 x9 x4 x5 x6 x7

Analysis of BKZ 16/32

slide-48
SLIDE 48

x1 x2 x3 x4 x9 x5 x6 x7 x8

Analysis of BKZ 16/32

slide-49
SLIDE 49

x1 x2 x3 x4 x5 x6 x7 x8 x9

Analysis of BKZ 16/32

slide-50
SLIDE 50

x5 x6 x7 x8 x9 x1 x2 x3 x4

Analysis of BKZ 16/32

slide-51
SLIDE 51

Matricial interpretation

x1 x2 x3 x4 x5 x6 x7 x8 x9 X = (x1, . . . , xn)T X0.5 ← A1X X1 ← A1X + Γ1 X2 ← A2X1 + Γ2 . . . Xk = AkXk + Γk with k = n − β + 1 A full tour: X ′ ← AX + Γ

Analysis of BKZ 17/32

slide-52
SLIDE 52

Matricial interpretation

x5 x6 x7 x8 x9 x1 x2 x3 x4 X = (x1, . . . , xn)T X0.5 ← A1X X1 ← A1X + Γ1 X2 ← A2X1 + Γ2 . . . Xk = AkXk + Γk with k = n − β + 1 A full tour: X ′ ← AX + Γ

Analysis of BKZ 17/32

slide-53
SLIDE 53

Matricial interpretation

x5 x6 x7 x8 x9 x1 x2 x3 x4 X = (x1, . . . , xn)T X0.5 ← A1X X1 ← A1X + Γ1 X2 ← A2X1 + Γ2 . . . Xk = AkXk + Γk with k = n − β + 1 A full tour: X ′ ← AX + Γ

Analysis of BKZ 17/32

slide-54
SLIDE 54

Matricial interpretation

x1 x6 x7 x8 x9 x2 x3 x4 x5 X = (x1, . . . , xn)T X0.5 ← A1X X1 ← A1X + Γ1 X2 ← A2X1 + Γ2 . . . Xk = AkXk + Γk with k = n − β + 1 A full tour: X ′ ← AX + Γ

Analysis of BKZ 17/32

slide-55
SLIDE 55

Matricial interpretation

x1 x2 x3 x4 x5 x6 x7 x8 x9 X = (x1, . . . , xn)T X0.5 ← A1X X1 ← A1X + Γ1 X2 ← A2X1 + Γ2 . . . Xk = AkXk + Γk with k = n − β + 1 A full tour: X ′ ← AX + Γ

Analysis of BKZ 17/32

slide-56
SLIDE 56

Matricial interpretation

x1 x2 x3 x4 x5 x6 x7 x8 x9 X = (x1, . . . , xn)T X0.5 ← A1X X1 ← A1X + Γ1 X2 ← A2X1 + Γ2 . . . Xk = AkXk + Γk with k = n − β + 1 A full tour: X ′ ← AX + Γ

Analysis of BKZ 17/32

slide-57
SLIDE 57

Expected properties of the model

X ← AX + Γ

Well-reduced output: → study of fixed points (X ∞ = AX ∞ + Γ). Convergence in a polynomial number of steps: → study of eigenvalues of ATA (so that AkX2 is bounded).

Analysis of BKZ 18/32

slide-58
SLIDE 58

Expected properties of the model

X ← AX + Γ

Well-reduced output: → study of fixed points (X ∞ = AX ∞ + Γ). Convergence in a polynomial number of steps: → study of eigenvalues of ATA (so that AkX2 is bounded).

Analysis of BKZ 18/32

slide-59
SLIDE 59

Expected properties of the model

X ← AX + Γ

Well-reduced output: → study of fixed points (X ∞ = AX ∞ + Γ). Convergence in a polynomial number of steps: → study of eigenvalues of ATA (so that AkX2 is bounded).

Analysis of BKZ 18/32

slide-60
SLIDE 60

Fixed point X ∞ - Uniqueness X ∞ = AX ∞ + Γ

What matters is the rank of A. The solutions of AX ∞ = X ∞ are vectors in Span(1, . . . , 1). Unique solution if we consider only {X| xi = 0}.

Analysis of BKZ 19/32

slide-61
SLIDE 61

Fixed point X ∞ - Uniqueness X ∞ = AX ∞ + Γ

What matters is the rank of A. The solutions of AX ∞ = X ∞ are vectors in Span(1, . . . , 1). Unique solution if we consider only {X| xi = 0}.

Analysis of BKZ 19/32

slide-62
SLIDE 62

Fixed point X ∞ - Uniqueness X ∞ = AX ∞ + Γ

What matters is the rank of A. The solutions of AX ∞ = X ∞ are vectors in Span(1, . . . , 1). Unique solution if we consider only {X| xi = 0}.

Analysis of BKZ 19/32

slide-63
SLIDE 63

Fixed point X ∞ - Existence

The last β vectors have the shape of an HKZ-reduced basis. Recursive formula for the previous vectors: x∞

i

= β 2(β − 1) log γβ +

i+β

  • j=i+1

x∞

j

β − 1. Asymptotically, line of slope − log γβ

β−1 .

i xi O((log β)2)

Analysis of BKZ 20/32

slide-64
SLIDE 64

Fixed point X ∞ - Existence

The last β vectors have the shape of an HKZ-reduced basis. Recursive formula for the previous vectors: x∞

i

= β 2(β − 1) log γβ +

i+β

  • j=i+1

x∞

j

β − 1. Asymptotically, line of slope − log γβ

β−1 .

i xi O((log β)2)

Analysis of BKZ 20/32

slide-65
SLIDE 65

Fixed point X ∞ - Existence

The last β vectors have the shape of an HKZ-reduced basis. Recursive formula for the previous vectors: x∞

i

= β 2(β − 1) log γβ +

i+β

  • j=i+1

x∞

j

β − 1. Asymptotically, line of slope − log γβ

β−1 .

i xi O((log β)2) ≃ (n − β)log γβ

β−1

Analysis of BKZ 20/32

slide-66
SLIDE 66

Eigenvalues of ATA

Method: study of the roots of the characteristic polynomial of ATA. Let χn(λ) = det(λIn − AT

n An).

Recurrence formula: χn+2(λ) = [2β(β − 1) + 1] λ − 1 β2 χn+1 − β − 1 β 2 λ2χn By a change of variable, it becomes a classical recurrence (Chebyshev polynomials): ψn+2(µ) = 2µψn+1(µ) − ψn(µ)

(change of variable: τ(µ) = 2β(β − 1)(µ − 1) et ψn(µ) = “

β β−1

”n−β ·

¯ χn(1−τ(µ)) τ(µ)

) Analysis of BKZ 21/32

slide-67
SLIDE 67

Eigenvalues of ATA

Method: study of the roots of the characteristic polynomial of ATA. Let χn(λ) = det(λIn − AT

n An).

Recurrence formula: χn+2(λ) = [2β(β − 1) + 1] λ − 1 β2 χn+1 − β − 1 β 2 λ2χn By a change of variable, it becomes a classical recurrence (Chebyshev polynomials): ψn+2(µ) = 2µψn+1(µ) − ψn(µ)

(change of variable: τ(µ) = 2β(β − 1)(µ − 1) et ψn(µ) = “

β β−1

”n−β ·

¯ χn(1−τ(µ)) τ(µ)

) Analysis of BKZ 21/32

slide-68
SLIDE 68

Eigenvalues of ATA

Method: study of the roots of the characteristic polynomial of ATA. Let χn(λ) = det(λIn − AT

n An).

Recurrence formula: χn+2(λ) = [2β(β − 1) + 1] λ − 1 β2 χn+1 − β − 1 β 2 λ2χn By a change of variable, it becomes a classical recurrence (Chebyshev polynomials): ψn+2(µ) = 2µψn+1(µ) − ψn(µ)

(change of variable: τ(µ) = 2β(β − 1)(µ − 1) et ψn(µ) = “

β β−1

”n−β ·

¯ χn(1−τ(µ)) τ(µ)

) Analysis of BKZ 21/32

slide-69
SLIDE 69

Explicit expression for ψn: ψn = Un−β+1 − β − 1 β Un−β with Un(cos x) = sin(nx)

sin x .

Studying this function leads to the following results:

1 is a simple root of the characteristic polynomial. The second largest eigenvalue of AT A is ≤ 1 − 1 2 β2 n2 .

Analysis of BKZ 22/32

slide-70
SLIDE 70

Explicit expression for ψn: ψn = Un−β+1 − β − 1 β Un−β with Un(cos x) = sin(nx)

sin x .

Studying this function leads to the following results:

1 is a simple root of the characteristic polynomial. The second largest eigenvalue of AT A is ≤ 1 − 1 2 β2 n2 .

Analysis of BKZ 22/32

slide-71
SLIDE 71

Results on the sandpile model

The slope − log γβ

β−1 of the fixed point corresponds to a Hermite

factor

b1 (det L)1/n close to γ

n−1 2(β−1)

β

. Geometric convergence: X − X ∞ decreases by a constant factor every n2

β2 tours, i.e. n3 β2 calls to HKZβ. n3 β2 (log n ǫ + log log max bi (det L)1/n ) calls to HKZβ are enough to

  • btain X − X ∞ < ǫ.

Analysis of BKZ 23/32

slide-72
SLIDE 72

Results on the sandpile model

The slope − log γβ

β−1 of the fixed point corresponds to a Hermite

factor

b1 (det L)1/n close to γ

n−1 2(β−1)

β

. Geometric convergence: X − X ∞ decreases by a constant factor every n2

β2 tours, i.e. n3 β2 calls to HKZβ. n3 β2 (log n ǫ + log log max bi (det L)1/n ) calls to HKZβ are enough to

  • btain X − X ∞ < ǫ.

Analysis of BKZ 23/32

slide-73
SLIDE 73

Results on the sandpile model

The slope − log γβ

β−1 of the fixed point corresponds to a Hermite

factor

b1 (det L)1/n close to γ

n−1 2(β−1)

β

. Geometric convergence: X − X ∞ decreases by a constant factor every n2

β2 tours, i.e. n3 β2 calls to HKZβ. n3 β2 (log n ǫ + log log max bi (det L)1/n ) calls to HKZβ are enough to

  • btain X − X ∞ < ǫ.

Analysis of BKZ 23/32

slide-74
SLIDE 74

1

Reminders on lattice reduction

2

Analysis of BKZ in the sandpile model

3

Analysis of BKZ

4

Applications to LLL

5

Conclusion

Analysis of BKZ 24/32

slide-75
SLIDE 75

Comparison between the model and BKZ

When the determinant is fixed, there is no vector inequality on the xi’s between: a worst-case HKZ-reduced basis (equalities in Minkowski inequalities) an arbitrary HKZ-reduced basis (strict inequalities). x5 x4 x3 x2 x1 → The previous results cannot be transposed directly.

Analysis of BKZ 25/32

slide-76
SLIDE 76

Comparison between the model and BKZ

When the determinant is fixed, there is no vector inequality on the xi’s between: a worst-case HKZ-reduced basis (equalities in Minkowski inequalities) an arbitrary HKZ-reduced basis (strict inequalities). x5 x4 x3 x2 x1 → The previous results cannot be transposed directly.

Analysis of BKZ 25/32

slide-77
SLIDE 77

Change of basis

Obtaining information on the individual xi’s is difficult. The model can give some information on πi = 1

i

i

j=1 xj, the

mean of the first xj’s. New dynamical system: Π ← AΠ + Γ ( A = PAP−1) In the real world, we still have Π ← Π′ ≤ AΠ + Γ (coefficient-wise).

Analysis of BKZ 26/32

slide-78
SLIDE 78

Change of basis

Obtaining information on the individual xi’s is difficult. The model can give some information on πi = 1

i

i

j=1 xj, the

mean of the first xj’s. New dynamical system: Π ← AΠ + Γ ( A = PAP−1) In the real world, we still have Π ← Π′ ≤ AΠ + Γ (coefficient-wise).

Analysis of BKZ 26/32

slide-79
SLIDE 79

Change of basis

Obtaining information on the individual xi’s is difficult. The model can give some information on πi = 1

i

i

j=1 xj, the

mean of the first xj’s. New dynamical system: Π ← AΠ + Γ ( A = PAP−1) In the real world, we still have Π ← Π′ ≤ AΠ + Γ (coefficient-wise).

Analysis of BKZ 26/32

slide-80
SLIDE 80

Change of basis

Obtaining information on the individual xi’s is difficult. The model can give some information on πi = 1

i

i

j=1 xj, the

mean of the first xj’s. New dynamical system: Π ← AΠ + Γ ( A = PAP−1) π1 π2 π3 π4 π5 In the real world, we still have Π ← Π′ ≤ AΠ + Γ (coefficient-wise).

Analysis of BKZ 26/32

slide-81
SLIDE 81

Results on BKZβ

Using the inequality Π′ ≤ AΠ + Γ recursively gives: Π[k] − Π∞ ≤ Ak(Π[0] − Π∞). The upper bound on the eigenvalues of ATA is used to bound the 2-norm of the right term. Π[k] − Π∞ ≤ (1 + log n)

1 2

  • 1 − β2

2n2 k

2

Π[0] − Π∞2

Analysis of BKZ 27/32

slide-82
SLIDE 82

Results on BKZβ

Using the inequality Π′ ≤ AΠ + Γ recursively gives: Π[k] − Π∞ ≤ Ak(Π[0] − Π∞). The upper bound on the eigenvalues of ATA is used to bound the 2-norm of the right term. Π[k] − Π∞ ≤ (1 + log n)

1 2

  • 1 − β2

2n2 k

2

Π[0] − Π∞2

Analysis of BKZ 27/32

slide-83
SLIDE 83

Meaning of the Πi’s: π1 = x1 = log b1 πn = n

i=1 xi = log det L

Π[k] − Π∞ ≤ (1 + log n)

1 2

  • 1 − β2

2n2 k

2

Π[0] − Π∞2 →

b1 (det L)1/n ≤ (1 + ǫ)γ

n−1 2(β−1)+ 3 2

β

in O( n2

β2 · n) calls to HKZβ.

Analysis of BKZ 28/32

slide-84
SLIDE 84

Meaning of the Πi’s: π1 = x1 = log b1 πn = n

i=1 xi = log det L

Π[k] − Π∞ ≤ (1 + log n)

1 2

  • 1 − β2

2n2 k

2

Π[0] − Π∞2 →

b1 (det L)1/n ≤ (1 + ǫ)γ

n−1 2(β−1)+ 3 2

β

in O( n2

β2 · n) calls to HKZβ.

Analysis of BKZ 28/32

slide-85
SLIDE 85

Differences between LLL and BKZ2

Swaps in LLL / HKZ2 = Gauss-reductions in BKZ2. → the complexity of both operations is O(size(B)). Different, non-adaptative order in BKZ2.

Analysis of BKZ 29/32

slide-86
SLIDE 86

Differences between LLL and BKZ2

Swaps in LLL / HKZ2 = Gauss-reductions in BKZ2. → the complexity of both operations is O(size(B)). Different, non-adaptative order in BKZ2.

Analysis of BKZ 29/32

slide-87
SLIDE 87

Differences between LLL and BKZ2

Swaps in LLL / HKZ2 = Gauss-reductions in BKZ2. → the complexity of both operations is O(size(B)). Different, non-adaptative order in BKZ2.

Analysis of BKZ 29/32

slide-88
SLIDE 88

Quasi-linear LLL

In BKZ2: Each Gauss-reduction costs O(log max bi). Poly(n) × log log maxi

b∗

i

(det L)1/n Gauss-reductions.

A basis such that

b1 (det L)1/n ≤

  • 4

3 n−1

(1 + ǫ) is returned. With more work, it is possible to obtain an LLL-reduced basis.

Analysis of BKZ 30/32

slide-89
SLIDE 89

1

Reminders on lattice reduction

2

Analysis of BKZ in the sandpile model

3

Analysis of BKZ

4

Applications to LLL

5

Conclusion

Analysis of BKZ 31/32

slide-90
SLIDE 90

Conclusion

The optimal quality that can be proven for BKZβ is reached in a polynomial number of calls to HKZβ. Binary complexity of BKZ2? Adaptive strategies. In practice, the algorithm reaches better approximation factors than expected. → For how long is it interesting to continue the execution

  • nce we go beyond the theorical factor?

Analysis of BKZ 32/32