A lone figure in a dark hall of towering bookshelves under a skylit oculus

The Geometry of SIS: Why Short Integer Solutions Are Hard to Find, and How n, m, q and β Set the Cost

Scope. A geometric account of the short integer solution problem (SIS), the assumption, in its module form, behind most commitments in the companion article on lattice-based polynomial commitment schemes, for readers who know what a lattice is and want the reasons it is believed hard. Theorems are labelled as theorems; the rest is the standard heuristic toolkit of lattice cryptanalysis, stated with its hypotheses, because that toolkit is what sets concrete parameters.

Standing on. Ajtai's definition and worst-case connection[1], the reduction and smoothing parameter of Micciancio and Regev[2], their analysis of attacks on q-ary lattices[3] and Peikert's survey[4]. The article's own contributions are small and labelled (§9): a two-sided bound that makes the existence threshold sharp, the attack length read as the height of a triangle of fixed area, and a check of the resulting one-line rule against a production commitment scheme's parameter table, 5,994 of whose 6,100 cells it reproduces.

Verification. Every number was computed by script; the existence bounds were checked on all 1,057,884 matrices of six small parameter sets. The figures recompute live.

Most lattice-based polynomial commitments start from the same line. The commitment to a vector s with small entries is t = As mod q for a public random matrix A, and it binds because two short openings of one t would hand over their difference z = s − s′: a short nonzero integer vector with Az ≡ 0. Finding such a vector is the short integer solution problem, and the binding of the Module-SIS schemes among them, Greyhound, Hachi and Akita included, rests on that problem at the parameters they choose.

It is an odd thing to rest on. The equation is linear, Gaussian elimination solves it at once, and short solutions exist in astronomical numbers; the only obstacle is the word short. What follows shows, mostly in pictures, why that metric condition costs an attacker time exponential in the dimension, how n, m, q and β enter through one combination, why random instances are as hard as the hardest lattice, and why quantum computers change only a constant.

The answer in one line

Lattice reduction with block size b leaves bases whose log-lengths fall with a fixed slope, 2 log2 δb per vector, and the determinant fixes the area under that profile at n log2 q. A right triangle with that slope and that area has height

log2(shortest vector reachable) ≈ 2 √(n · log2 q · log2 δb),

so breaking SIS at norm β needs log2 δb ≤ (log2 β)2/(4n log2 q). The block size that achieves this grows like n times a logarithm, and each block is a shortest-vector search costing 20.292b classically and 20.265b with a quantum computer. That is where the exponential comes from.

1 · The problem, and the two ways it becomes trivial

Given a uniformly random A ∈ ℤqn×m with m > n, find a nonzero z ∈ ℤm with

Az ≡ 0 (mod q),    ‖z‖ ≤ β.SISn,m,q,β

This is Ajtai's problem as Peikert's survey states it[1][4], in the Euclidean norm unless said otherwise. Without the norm the solutions form a lattice,

Λ⊥(A) = { z ∈ ℤm : Az ≡ 0 (mod q) } ⊇ qℤm.

For prime q and A of rank n, which fails with probability below qn−m/(q − 1), the map z ↦ Az mod q takes ℤm onto ℤqn with kernel Λ⊥, so the lattice has index and determinant qn. The solution set is the integer grid thinned to one point in qn: each row of A keeps one point in q.

Two observations bracket the problem[4]. Gaussian elimination writes A as [In | Ā] and hands over a basis at once: qe1, …, qen and the vectors (−Āej, ej), whose first n entries are uniform in (−q/2, q/2], giving length about q√(n/12). And with β ≥ q the problem is trivial, because qe1 is a solution. SIS lives in between. In the regimes that matter the lattice contains vectors far shorter than q (§2), and the task is to find one when the algebra provides only long ones.

drag to rotate
11
2
3.6
lattice points in the cube lattice points in the ball, and a shortest vector basis from Gaussian elimination; rings mark qei ball of radius β
integer points kept by A
1 in 121
lattice points in the cube
—
shortest length λ1
—
Gaussian heuristic (3qn/4π)1/3
—
integer points in the ball N(β)
—
solutions in the ball vs (N − 1)/qn
—
Figure 1. Λ⊥(A) for m = 3: the points of the cube [−q, q]3 with Az ≡ 0, for the A = [In | Ā] printed top left. Amber arrows are the elimination basis; the green arrow is a shortest vector, found by checking every point. The tiles compare the solutions in the ball with the (N − 1)/qn of §2. In three dimensions inspection finds the shortest vector; the difficulty appears in hundreds.

2 · Existence: a ball against a determinant

How small can β be before solutions stop existing? A nonzero z with every entry below q in absolute value has an entry zi that is a unit modulo the prime q, so Az = Σj zjaj is uniform on ℤqn when A is. With N(β) the number of integer points in the closed ball of radius β < q, linearity of expectation gives an identity, not an estimate:

𝔼A #{ z ≠ 0 : ‖z‖ ≤ β, Az ≡ 0 } = (N(β) − 1) / qn.

The ball holds about Vmβm integer points (Vm the volume of the unit ball), so the expected count passes 1 where the ball's volume meets the determinant,

β* = (qn/Vm)1/m ≈ √(m/2πe) · qn/m,Gaussian heuristic

the estimate of λ1(Λ⊥) that Micciancio and Regev use, reliable when m is neither close to n nor very large[3]. A second-moment argument turns the expectation into a sharp threshold.

Proposition · a two-sided existence window

Let q be an odd prime, 1 ≤ β < q/2 and N = N(β). Then

1 − 2⌊β⌋ qn/(N − 1) ≤ PrA[ a solution exists ] ≤ (N − 1)/(2qn).

Proof. Group the nonzero points of the ball into classes z ∼ cz (mod q), c ∈ ℤq×. Within a class the events Az ≡ 0 coincide; points of different classes are linearly independent modulo q, so their events are pairwise independent, each of probability q−n. A class contains z and −z, and at most 2⌊β⌋ points: since 2β < q, a point of the ball is fixed by its residue, and the members are indexed by the centred value of czi in [−β, β] ∖ {0} for a unit coordinate zi. So the number of classes K lies between (N − 1)/2⌊β⌋ and (N − 1)/2. Let X count the classes in the kernel. Markov's inequality gives Pr[X ≥ 1] ≤ 𝔼X = K/qn; pairwise independence gives Var X ≤ 𝔼X, and Chebyshev's inequality gives Pr[X = 0] ≤ 1/𝔼X. ∎

Where N(β) tracks the ball's volume, log N(β) grows like m log β, and the probability climbs from ε to 1 − ε within a factor 1 + O((log β + log 1/ε)/m) in β: a random q-ary lattice of high dimension with β* < q/2 has its shortest vector at β*, up to a vanishing relative error. For very wide matrices the ball holds far more integer points than its volume, and the threshold falls below β*.

move q, n and m
127
6
20
enumerated exactly:
Gaussian heuristic β*
—
smallest β with expected count ≥ 1
—
Pr ≤ 5% up to β =
—
Pr ≥ 95% from β =
—
exact enumeration
—
Figure 2. Top: log2 of the expected number of solutions, (N(β) − 1)/qn, with N counted exactly; it crosses zero near the Gaussian-heuristic radius β* (gold) unless m is large. Bottom: the proposition's band (blue) for the probability that a solution exists. The buttons load six fully enumerated parameter sets, whose exact probabilities appear as green dots inside the band. Raise m and the band closes into a step.

Each row multiplies the determinant by q; each column multiplies the ball's volume by about β√(2π/m), so columns win quickly. Once m > n log2 q, two of the 2m vectors in {0, 1}m collide under A, and their difference is a {0, ±1} solution of norm at most √(n log2 q + 1)[4]. Commitment matrices are far wider. Short solutions exist in enormous numbers; whatever makes SIS hard, it is not scarcity.

3 · Why search is exponential: the sphere in high dimension

In the plane, Lagrange and Gauss's reduction finds a shortest vector in a few steps. The difficulty is dimension, and its source is a fact about spheres. Two independent uniform directions on Sd−1 meet at an angle θ with density proportional to sind−2 θ, so in high dimension almost every pair is nearly orthogonal, and a cap of angle θ ≤ 90° holds (sin θ)d + o(d) of the sphere[5]. The angle that matters is 60°: two vectors of equal length have a difference shorter than either exactly when their angle is below 60°, and that cap holds (√3/2)d + o(d) = 2−0.2075d + o(d) of the sphere.

The fastest known algorithms for the shortest vector problem in high dimension are sieves, and this cap is their cost. A sieve keeps a list of lattice vectors of about equal length and replaces pairs by their differences, which are shorter only for pairs within 60°. If list vectors point in random directions, the standard heuristic, each finds a partner only when the list holds about the reciprocal of the cap, (4/3)d/2 = 20.2075d vectors, and comparing all pairs costs 20.415d[6]. Faster sieves search the pairs more cleverly: locality-sensitive filtering reaches 20.292d[7], quantum search inside the filters 20.2653d[5] and quantum random walks 20.2570d[8]. None goes below the list. The Dilithium designers call 20.2075d the paranoid floor, since the best solvers cover the sphere with that many cones of half-angle π/3[9].

Exponent per dimensionClosed formWhat it counts
0.2075log₂ √(4/3)Cones of half-angle 60° covering the sphere: the sieve's list, and the floor under every known sieve
0.2653log₂ √(13/9)The filtered sieve with Grover search
0.2925log₂ √(3/2)The sieve with locality-sensitive filtering
0.4150log₂ (4/3)Comparing every pair in the list

The constants are base-2 logarithms of √(4/3), √(3/2) and √(13/9), numbers of spherical geometry rather than of any implementation. Enumeration needs only polynomial memory but 2Θ(d log d) time, and no known algorithm, classical or quantum, approximates shortest vectors in arbitrary lattices within polynomial factors in time 2o(d)[4].

drag the sphere
60°
400
cap's share of the sphere at d
—
directions per useful partner
—
60° floor 20.2075d
—
all pairs 20.415d
—
filtered sieve 20.292d
—
quantum sieve 20.2653d
—
Figure 3. Top left: 600 uniform points on the sphere in three dimensions; the blue ones lie inside the amber cap of angle θ, a quarter of the sphere at 60°. Top right: two unit vectors at angle θ differ by 2 sin(θ/2), which is below 1 only when θ < 60°. Bottom: the exact cap fraction up to dimension 500 (blue) against its asymptote d log2 sin θ (dashed); in dimension 400 about one direction in 287 is a useful partner. The tiles give the sieve costs at the marked dimension, set by the 60° cap rather than by θ.

4 · Reduction is a slope, the determinant is an area

An attacker does not solve the shortest vector problem in dimension m. BKZ with block size b runs an exact solver in dimension b on projected blocks and leaves a basis whose Gram–Schmidt log-lengths ℓi = log2 ‖bi*‖ fall almost linearly, the geometric series assumption, with slope

2 log2 δb,    δb ≈ ( b/(2πe) · (πb)1/b )1/(2(b−1)).root-Hermite factor

The formula is the Dilithium specification's[9]; at b = 85 to 286 it reproduces Chen and Nguyen's simulations, δ85 = 1.0100 and δ286 = 1.0050[10], while LLL alone reaches about 1.0219[11]. Two more facts decide everything.

The area is conserved. Σi ℓi = log2 det Λ⊥ = n log2 q for every basis. Reduction moves length from the front of the profile to the back, never creating or destroying any.

The attacker chooses the base. Columns of A can be ignored by fixing their coordinates to zero[4], so an attacker may reduce the lattice of any w ≤ m columns, whose determinant is still qn for w > n.

The reduced profile on w columns is therefore a line of fixed slope and fixed area, and the attacker wants its first value ℓ1, the log-length of the vector it gets, as small as possible. At a given height, a line that stops exactly at zero encloses the most area, so the best profile is a right triangle of height h, base h/(2 log2 δb) and area h2/(4 log2 δb). Setting the area to n log2 q,

h = 2 √(n log2 q · log2 δb)  at  w* = √(n log2 q / log2 δb) columns,

which is Micciancio and Regev's estimate[3]. Fewer columns force the line to start higher; more buy nothing, and an attacker drops them, as Chen and Nguyen did on Darmstadt's SIS challenges of dimension 525 to 800, reducing sublattices of dimension 180 to 230[10]. Run on the q-ary basis, the profile also never rises above log2 q, since the qei stay available: when h would exceed log2 q they stay in front, the Z-shaped profile of the Dilithium analysis[9], and nothing shorter than q is found.

The attack in one condition

Under these heuristics, BKZ-b finds a solution of norm at most β when

min( log2 q,  minn < w ≤ m [ w log2 δb + n log2 q / w ] ) ≤ log2 β,

and the core-SVP convention prices it as one call to the sieve in dimension b: 20.292b classically, 20.265b with quantum search[12][9]. BKZ makes many such calls; charging one is a convention, not a proven bound[9].

move the parameters and the block size
64
16
6.0
1,024
150
root-Hermite factor δb
—
slope per vector (bits)
—
best width w*
—
length reached vs β
—
BKZ-b on this instance
—
smallest b that breaks it
—
its cost, bits classical / quantum
—
Figure 4. Top: the Gram–Schmidt profile of the best sublattice before reduction (amber: n vectors of length q, then vectors of length 1) and after BKZ-b (blue): the same area at a constant slope, clipped at log2 q. Its first value is the vector the attacker obtains (gold); the green dashed line is the target. Bottom: the attacker's valley, the length reached against the number of columns kept, with its floor at w* (violet) and the available columns ending at m (grey). Where it dips below the target, shaded green, BKZ-b breaks the instance.

5 · How n, m, q and β set the price

Squaring the condition puts the parameter story in one inequality. With m ≥ w* and β < q, BKZ-b breaks SISn,m,q,β when

(log2 β)2 ≥ 4n · log2 q · log2 δb.

ParameterExistence (§2)Attack (§4)Net effect on the cost
n, rowsMultiplies the determinant by q per rowEnlarges the triangle's area; its height grows as √nThe block size needed grows at least linearly in n, and the cost exponentially in the block size
q, modulusAs n, through log qEnlarges the area, and lifts the trivial ceiling β < qHarder at fixed β, easier at fixed β = qc: log β against √(n log q) is what counts
β, norm boundLarger ball, more solutionsThe target line; the attack wins when the triangle's height falls below itAt fixed security, n must grow as (log β)2
m, columnsMore columns, far more solutionsIrrelevant once m ≥ w*; below it the attacker is pinned to the valley's left wallMatters only when small; wide matrices are free

The exponential shows in numbers. At q = 232 and β = 216 the attack needs block size 252 at n = 256, 685 at n = 512 and 1,703 at n = 1,024: about 74, 200 and 497 bits of classical core-SVP. Each doubling of n about doubles b/log2(b/2πe), which the condition makes roughly proportional to n log2 q/(log2 β)2.

The modulus cuts both ways. At fixed β a larger q enlarges the area: with n = 512 and β = 212 the block size needed is 583, 1,465 and 3,497 at q = 216, 232 and 264. At β = √q a larger q makes the problem easier, 1,703, 685 and 252, because the target grows faster than the square root of the area. Hardness is log β against √(n log q), never q alone. Figure 5 lays out the plane of parameters.

click the upper plot to probe a point
32
∞
128
breakable within the budget secure at the budget under one solution expected (finite m) Akita table, ring degree 32 ring degree 64
probe
—
block size needed
—
bits, classical / quantum
—
rows needed for the budget
—
Figure 5. Top: rows n against norm bits log2 β at the chosen modulus. Red instances fall to the budget, blue ones do not, and grey ones expect fewer than one solution at the chosen m (§2); thin lines are the 64- to 256-bit contours. Amber diamonds and gold squares are Akita's production table for q = 232 − 99 at 128-bit quantum security, the smallest rank secure at every width times the ring degree (§8); at the defaults they trace the white contour. Bottom: bits against n for three norm bounds, bending upward because b grows like n log n.

6 · Why a random instance is as hard as the worst case

The case that random matrices hide no shortcut rests on a theorem with a geometric engine.

Theorem · Ajtai; Micciancio–Regev; Gentry–Peikert–Vaikuntanathan; Micciancio–Peikert

For any m = poly(n), any β > 0 and any sufficiently large q ≥ β · poly(n), solving SISn,m,q,β on a random A with non-negligible probability is at least as hard as solving GapSVPγ and SIVPγ on arbitrary n-dimensional lattices, for some γ = β · poly(n)[4]. Micciancio and Regev reach γ = β · Õ(√n) with q as small as β · Õ(n√m)[2]; Gentry, Peikert and Vaikuntanathan lower the modulus to β · Õ(√n)[13], and Micciancio and Peikert to β · nε for any constant ε > 0, with γ then also depending on ‖z‖∞[14].

The reduction must turn an arbitrary lattice L into a uniformly random matrix, and it does so by blurring. Sample lattice vectors vi with a wide Gaussian spread, reduce them modulo qL and read their coordinates modulo q as the columns ai. Once the spread exceeds the smoothing parameter of qL, defined below, the residues are within statistical distance 2ε of uniform[13]. The reason is Poisson summation: a Gaussian of width s summed over a lattice has Fourier coefficient e−πs2‖w‖2 at each dual lattice point w, so

ηε(L) = min { s : Σw ∈ L* ∖ {0} e−πs2‖w‖2 ≤ ε }smoothing parameter

is the blur at which all the nonzero frequencies together carry weight ε. A continuous blur is then within a factor 1 ± ε of uniform at every point and within statistical distance ε/2 of it[2]: the lattice has disappeared. Micciancio and Regev bound ηε(L) by √(ln(2n(1 + 1/ε))/π) · λn(L)[2]; for ℤ2 the bound is nearly exact, 1.3821 against a true η0.01 of 1.3813.

blur the lattice
0.80
ℤ²
ε = nonzero dual weight
—
largest |density − 1| on screen
—
distance from uniform ≤ ε/2
—
η0.01 of this lattice
—
s / η0.01
—
Figure 6. Left: a two-dimensional lattice of determinant 1 (amber points, one cell outlined) blurred by a Gaussian of width s (white circle), coloured by density relative to uniform. Right: the dual lattice, each point drawn with area proportional to its Fourier weight; the violet origin always carries weight 1. As s grows every other weight collapses and the density goes flat; past η0.01 it is flat to within 1%, whichever lattice produced it.

The reduction keeps a set S of lattice vectors and uses it to sample the vi, with spread ‖S‖ · poly(n). The oracle sees a uniform A and returns a short z; then Σ zivi lies in qL, so v = Σ zivi/q is a vector of L of length at most ‖S‖ · β · poly(n)/q, below ‖S‖/2 once q ≥ β · poly(n). Repeating shrinks S to within γ of λn(L)[4]. A random SIS instance is a blurred photograph of an arbitrary lattice, and an algorithm that solves the photographs sharpens every lattice.

What this buys is qualitative: no noticeable share of weak instances short of an algorithm for worst-case lattice problems within polynomial factors, which the best known algorithms, classical or quantum, take time 2Θ(n) to solve[4]. It does not buy parameters; the reduction loses polynomial factors, so concrete security comes from §3 to §5.

7 · A quantum computer changes a constant, not the exponent

Shor's algorithm breaks factoring and discrete logarithms by finding a hidden period in an abelian group[15], and Regev observed that almost all known super-polynomial quantum speed-ups solve abelian cases of that hidden-subgroup problem[16]. SIS offers no period: its solution lattice is public, and what is hidden is a short element, a metric search. For general lattices the known bridge to hidden subgroups runs through the dihedral group, which is not abelian: Regev reduced poly(n)-unique-SVP to a dihedral coset problem[16], for which no polynomial-time quantum algorithm is known.

What quantum computers do provide is search. The quantum sieves of §3 lower the exponent from 0.292b to 0.2653b or 0.2570b (both assuming quantum RAM), a smaller constant in a cost still exponential in the block size. The reduction of §6 is classical and black-box, so it holds for quantum solvers too: a quantum algorithm for random SIS would be one for worst-case SIVP within polynomial factors, and none is known. A 2024 preprint claimed such an algorithm by way of the learning-with-errors problem; its author's update reports a bug in Step 9 of the algorithm and states that the claim does not hold[17].

Ideal lattices have one documented soft spot, and abelian structure is its source: quantum class-group and principal-ideal computations find vectors within a factor exp(Õ(√n)) of the shortest in ideals of cyclotomic fields in polynomial time, under a number-theoretic hypothesis[18]. That factor is subexponential, far above the polynomial factors in the worst-case reductions for Ring- and Module-SIS[4][19], and the Dilithium analysis notes that moving from rings to modules of rank above one takes a scheme further from those lattices[9].

8 · What a lattice commitment pays: the norm budget

Binding is SIS at the collision norm. Two openings s ≠ s′ of one commitment, each of norm at most βo, give the solution s − s′ of norm at most 2βo. The schemes of the companion article accept relaxed openings and folded responses, so their collision norm also carries the challenges; Akita's security notes give its Euclidean form as 8γ√S, with γ bounding multiplication by a challenge and S the proved squared norm of the response[20]. Whatever the slack, it lands in log β, and the condition of §5 prices it:

n ≥ (log2 β)2 / (4 log2 q · log2 δb).the quadratic law

With the modulus fixed by the proof system's field, the scalar rank grows with the square of the number of bits in the collision norm, so each extra bit of slack costs more rank than the last; the table gives the rows at q = 232 − 99 and 128-bit quantum security, b = 128/0.265 ≈ 483. This is why every scheme re-decomposes its witness into small digits after each fold, and why the slack-free norm proofs of Serval, Maltese and Akita are worth their complexity.

collision normrows by the rulerank × 64 by the ruleAkita's table, ring degree 64
281002 × 642 × 64
2163997 × 647 × 64
22489614 × 6414 × 64
2281,22020 × 6420 × 64

The rule is two lines; production estimators are not. Akita's estimator prices every table cell offline, searching attack parameters under the ADPS16 core-SVP model at 0.265b and certifying, for each collision norm, ring degree and module rank, the widest matrix that keeps 128-bit quantum security[20]. Its checked-in Euclidean table for q = 232 − 99 has 6,100 cells: ring degrees 32 to 512, collision norms 20.5 to 230.5 in half-bit steps, ranks 1 to 20. Call a rank secure at every width when the valley floor of §4 stays above the target, and otherwise allow widths up to the valley's left crossing. The two-line rule then reproduces 5,994 cells exactly. Of the other 106, 95 allow one or two ring elements less width than the table, ten allow one more, and one is a rank the rule calls secure at every width where the table caps it at 36 ring elements; the rule's bound there is 223.9 rows against a boundary of 224.

Two qualifications travel with it. Like the rule, the estimator treats a Module-SIS matrix of rank r over a ring of degree D as unstructured SIS with n = rD; Akita's notes state that it “does not price attacks that exploit ring or module structure”[20], the position of §7. And it searches the same core-SVP and geometric-series model, so the agreement confirms the closed form, not the model: both stand or fall with the heuristics of §4.

9 · Provenance

ProvenanceWhat
Standard, from the cited papersSIS and its trivialities; the Gaussian heuristic; the sieve exponents; BKZ and the geometric series assumption; w* and the length 22√(n log q log δ); core-SVP; the worst-case theorem and smoothing; the quantum results.
This article's framingThe attack length as the height of a triangle of fixed area, a geometric proof of an optimum MR09 states without derivation; the two-sided existence window, from folklore ingredients; the quadratic law for commitment designers; all figures.
Checked against dataThe block-size formula against Chen and Nguyen's Table 2 (six rows, within 2·10−5); core-SVP against Dilithium's Table 1; MR09's Figure 2 optimum (w* ≈ 392); the existence window against every matrix of six small sets; the quadratic law against Akita's table (5,994 of 6,100 cells).

10 · What to keep

  • Existence is cheap and sharp. (N(β) − 1)/qn is exact, the threshold sits where N(β) reaches qn, and wide matrices have solutions to spare; only finding them is hard.
  • Search is expensive because of spheres. In dimension d a useful partner lies in a cap of measure 2−0.2075d + o(d), and every known sieve pays a power of that.
  • The attack is a triangle. Reduction fixes the slope and the determinant fixes the area, so log2(reachable length) ≈ 2√(n log2 q · log2 δb): n buys exponential cost, β is paid quadratically in bits, q cuts both ways, m stops mattering past w*.
  • Random instances are blurred worst cases, and quantum computers improve only the constant in the exponent.
  • For a commitment, the binding norm is the budget: rank grows with the square of the collision norm's bit length, which is why the schemes re-decompose after each fold.

References

Published versions are cited by venue with a DOI or an author copy; preprints by ePrint or arXiv number. Each entry names the definition, lemma, table or figure the text relies on.

  1. M. Ajtai, "Generating Hard Instances of Lattice Problems," STOC 1996, 99–108 — doi.org/10.1145/237814.237838. The short integer solution problem, the function z ↦ Az, and the first worst-case to average-case connection.
  2. D. Micciancio, O. Regev, "Worst-Case to Average-Case Reductions Based on Gaussian Measures," SIAM Journal on Computing 37(1):267–302, 2007 (FOCS 2004) — cims.nyu.edu/~regev/papers/average.pdf. Definition 3.1 (the smoothing parameter), Lemma 3.3 (ηε(Λ) ≤ √(ln(2n(1 + 1/ε))/π) · λn(Λ)), Lemma 4.1 (statistical distance at most ½ρ1/s(L* ∖ {0})) and footnote 12 (the pointwise 1 ± ε form); the reduction with γ = β · Õ(√n), as summarised in Peikert's survey §4.1.2.
  3. D. Micciancio, O. Regev, "Lattice-based Cryptography," in D. J. Bernstein, J. Buchmann, E. Dahmen (eds.), Post-Quantum Cryptography, Springer, 2009, 147–191 — cims.nyu.edu/~regev/papers/pqc.pdf. Section 3: the estimate λ1(Λ⊥) ≈ qn/m√(m/2πe), the length min{q, qn/mδm} reached by lattice reduction (eq. 1), its minimum 22√(n log q log δ) at m = √(n log q/log δ) (eq. 2) and Figure 2.
  4. C. Peikert, "A Decade of Lattice Cryptography," Foundations and Trends in Theoretical Computer Science 10(4):283–424, 2016; ePrint 2015/939 — eprint.iacr.org/2015/939. Definition 4.1.1 and the observations after it (Gaussian elimination, β < q, monotonicity in m and n, the pigeonhole bound); Theorem 4.1.2, the improvements of MR04, GPV08 and MP13, and the overview of the reduction (§4.1.2); the state of classical and quantum algorithms (§2.2.2).
  5. T. Laarhoven, Search Problems in Cryptography: From Fingerprinting to Lattice Sieving, PhD thesis, Eindhoven University of Technology, 2016 — thijs.com/docs/phd-final.pdf. The angle density (sin θ)d−2 between random directions (eq. 10.13) and the cap measure poly(d) · (sin θ)d (Lemma 13.2); the Nguyen–Vidick list of (4/3)d/2 vectors and time 20.415d (§8.2.1); hypercone filtering at 20.292d classically and 20.2653d quantumly (§14.2.10).
  6. P. Q. Nguyen, T. Vidick, "Sieve Algorithms for the Shortest Vector Problem Are Practical," Journal of Mathematical Cryptology 2(2):181–207, 2008 — doi.org/10.1515/JMC.2008.009.
  7. A. Becker, L. Ducas, N. Gama, T. Laarhoven, "New Directions in Nearest Neighbor Searching with Applications to Lattice Sieving," SODA 2016, 10–24; ePrint 2015/1128 — eprint.iacr.org/2015/1128. Locality-sensitive filtering; heuristic time 20.292d + o(d).
  8. A. Chailloux, J. Loyer, "Lattice Sieving via Quantum Random Walks," ASIACRYPT 2021, 63–91; ePrint 2021/570 — eprint.iacr.org/2021/570. Heuristic time 20.2570d + o(d) in the QRAM model, against Laarhoven's 20.2653d + o(d) in the same model.
  9. S. Bai, L. Ducas, E. Kiltz, T. Lepoint, V. Lyubashevsky, P. Schwabe, G. Seiler, D. Stehlé, CRYSTALS-Dilithium: Algorithm Specifications and Supporting Documentation, version 3.1, 2021 — pq-crystals.org/dilithium/data/dilithium-specification-round3-20210208.pdf. Appendix C.1: the constants log2√(3/2), log2√(13/9) and the “paranoid” log2√(4/3) with its cone-covering argument, the slope formula, and Figure 11 (the three-zone profile); Table 1 (b = 423 against SIS: 123 classical and 112 quantum core-SVP bits); Appendix C.4 (algebraic attacks, and module rank above one); Appendix C.5 (why the sign of the o(b) term in core-SVP is unclear).
  10. Y. Chen, P. Q. Nguyen, "BKZ 2.0: Better Lattice Security Estimates," ASIACRYPT 2011, 1–20 — doi.org/10.1007/978-3-642-25385-0_1. Table 2 (block sizes 85, 106, 133, 168, 216 and 286 for root-Hermite factors 1.010 down to 1.005); §5.1 and Table 1 (Darmstadt lattice challenges of dimension 525–800, solved or improved through sublattices of dimension 180–230).
  11. N. Gama, P. Q. Nguyen, "Predicting Lattice Reduction," EUROCRYPT 2008, 31–51 — doi.org/10.1007/978-3-540-78967-3_3. Table 1: average root-Hermite factors 1.0219 for LLL and 1.0128 for BKZ-20 on random lattices.
  12. E. Alkim, L. Ducas, T. Pöppelmann, P. Schwabe, "Post-quantum Key Exchange — A New Hope," USENIX Security 2016, 327–343; ePrint 2015/1092 — eprint.iacr.org/2015/1092. The core-SVP hardness measure.
  13. C. Gentry, C. Peikert, V. Vaikuntanathan, "Trapdoors for Hard Lattices and New Cryptographic Constructions," STOC 2008, 197–206; ePrint 2007/432 — eprint.iacr.org/2007/432. Corollary 2.8: a discrete Gaussian on Λ reduced modulo a sublattice Λ′ is within statistical distance 2ε of uniform once s ≥ ηε(Λ′).
  14. D. Micciancio, C. Peikert, "Hardness of SIS and LWE with Small Parameters," CRYPTO 2013, 21–39; ePrint 2013/069 — eprint.iacr.org/2013/069.
  15. P. W. Shor, "Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer," SIAM Journal on Computing 26(5):1484–1509, 1997 — doi.org/10.1137/S0097539795293172.
  16. O. Regev, "Quantum Computation and Lattice Problems," SIAM Journal on Computing 33(3):738–760, 2004 (FOCS 2002); arXiv cs/0304005 — arxiv.org/abs/cs/0304005. The reduction from unique-SVP to the dihedral hidden subgroup problem, and the remark (§1) that almost all known super-polynomial quantum speed-ups solve abelian hidden subgroup problems.
  17. Y. Chen, "Quantum Algorithms for Lattice Problems," ePrint 2024/555 — eprint.iacr.org/2024/555. The author's note on the ePrint page: Step 9 of the algorithm contains a bug, and the claimed polynomial-time quantum algorithm for LWE does not hold.
  18. R. Cramer, L. Ducas, B. Wesolowski, "Short Stickelberger Class Relations and Application to Ideal-SVP," EUROCRYPT 2017, 324–348; ePrint 2016/885 — eprint.iacr.org/2016/885. Ideal-SVP in cyclotomic fields within exp(Õ(√n)) in quantum polynomial time under a number-theoretic hypothesis; the abstract notes that Ring-LWE-based systems do not seem directly affected.
  19. A. Langlois, D. Stehlé, "Worst-Case to Average-Case Reductions for Module Lattices," Designs, Codes and Cryptography 75(3):565–599, 2015 — doi.org/10.1007/s10623-014-9938-4. Module-SIS and Module-LWE with worst-case reductions over module lattices.
  20. Q. Dao, O. Bodaghi, A. Khajehpour, G. Vitto, M. Badakhshan, M. Georghiades, F. Liu, J. Zhang, J. Thaler, "Akita: A High-Performance Lattice-Based Polynomial Commitment Scheme," ePrint 2026/1983 — eprint.iacr.org/2026/1983; and the Akita repository at commit 3a45521 — github.com/LayerZero-Labs/akita: book chapters “Lattices and Module-SIS” and “Security model” (the ADPS16 quantum policy at 0.2650b, the collision bounds 8κ1Z∞ and 64γ2S, and the statement that the estimator does not price ring or module structure), and the generated table crates/akita-params/src/sis/generated_l2_sis_table/q32.rs.