Short integer solution problem

From HandWiki
Revision as of 16:16, 6 February 2024 by Gametune (talk | contribs) (simplify)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Short integer solution (SIS) and ring-SIS problems are two average-case problems that are used in lattice-based cryptography constructions. Lattice-based cryptography began in 1996 from a seminal work by Miklós Ajtai[1] who presented a family of one-way functions based on SIS problem. He showed that it is secure in an average case if the shortest vector problem [math]\displaystyle{ \mathrm{SVP}_\gamma }[/math] (where [math]\displaystyle{ \gamma = n^c }[/math] for some constant [math]\displaystyle{ c\gt 0 }[/math]) is hard in a worst-case scenario. Average case problems are the problems that are hard to be solved for some randomly selected instances. For cryptography applications, worst case complexity is not sufficient, and we need to guarantee cryptographic construction are hard based on average case complexity.

Lattices

A full rank lattice [math]\displaystyle{ \mathfrak{L} \subset \R^n }[/math] is a set of integer linear combinations of [math]\displaystyle{ n }[/math] linearly independent vectors [math]\displaystyle{ \{b_1,\ldots,b_n\} }[/math], named basis:

[math]\displaystyle{ \mathfrak{L}(b_1,\ldots,b_n) = \left\{ \sum_{i=1}^n z_i b_i: z_i \in \Z \right\} = \{ B\boldsymbol{z}: \boldsymbol{z} \in \Z^n \} }[/math]

where [math]\displaystyle{ B \in \R ^{n\times n} }[/math] is a matrix having basis vectors in its columns.

Remark: Given [math]\displaystyle{ B_1,B_2 }[/math] two bases for lattice [math]\displaystyle{ \mathfrak{L} }[/math], there exist unimodular matrices [math]\displaystyle{ U_1 }[/math] such that [math]\displaystyle{ B_1 = B_2U_1^{-1}, B_2 = B_1U_1 }[/math].

Ideal lattice

Definition: Rotational shift operator on [math]\displaystyle{ \R^n (n \geq 2) }[/math] is denoted by [math]\displaystyle{ \operatorname{rot} }[/math], and is defined as:

[math]\displaystyle{ \forall \boldsymbol{x} = (x_1,\ldots,x_{n-1},x_n) \in \R^n: \operatorname{rot}(x_1,\ldots,x_{n-1},x_n) = (x_n,x_1,\ldots,x_{n-1}) }[/math]

Cyclic lattices

Micciancio introduced cyclic lattices in his work in generalizing the compact knapsack problem to arbitrary rings.[2] A cyclic lattice is a lattice that is closed under rotational shift operator. Formally, cyclic lattices are defined as follows:

Definition: A lattice [math]\displaystyle{ \mathfrak{L} \subseteq \Z^n }[/math] is cyclic if [math]\displaystyle{ \forall \boldsymbol{x} \in \mathfrak{L}: \operatorname{rot}(\boldsymbol{x}) \in \mathfrak{L} }[/math].

Examples:[3]

  1. [math]\displaystyle{ \Z^n }[/math] itself is a cyclic lattice.
  2. Lattices corresponding to any ideal in the quotient polynomial ring [math]\displaystyle{ R = \Z[x]/(x^n-1) }[/math] are cyclic:

consider the quotient polynomial ring [math]\displaystyle{ R = \Z[x]/(x^n-1) }[/math], and let [math]\displaystyle{ p(x) }[/math] be some polynomial in [math]\displaystyle{ R }[/math], i.e. [math]\displaystyle{ p(x) = \sum_{i=0}^{n-1}a_ix^i }[/math] where [math]\displaystyle{ a_i \in \Z }[/math] for [math]\displaystyle{ i = 0,\ldots, n-1 }[/math].

Define the embedding coefficient [math]\displaystyle{ \Z }[/math]-module isomorphism [math]\displaystyle{ \rho }[/math] as:

[math]\displaystyle{ \begin{align} \quad \rho: R & \rightarrow \Z^n \\[4pt] \rho(x) = \sum_{i=0}^{n-1}a_ix^i & \mapsto (a_0,\ldots,a_{n-1}) \end{align} }[/math]

Let [math]\displaystyle{ I \subset R }[/math] be an ideal. The lattice corresponding to ideal [math]\displaystyle{ I \subset R }[/math], denoted by [math]\displaystyle{ \mathfrak{L}_I }[/math], is a sublattice of [math]\displaystyle{ \Z^n }[/math], and is defined as

[math]\displaystyle{ \mathfrak{L}_I := \rho(I) = \left\{ (a_0,\ldots,a_{n-1}) \mid \sum_{i=0}^{n-1}a_ix^i \in I \right\} \subset \Z^n. }[/math]

Theorem: [math]\displaystyle{ \mathfrak{L} \subset \Z^n }[/math] is cyclic if and only if [math]\displaystyle{ \mathfrak{L} }[/math] corresponds to some ideal [math]\displaystyle{ I }[/math] in the quotient polynomial ring [math]\displaystyle{ R = \Z[x]/(x^n-1) }[/math].

proof: [math]\displaystyle{ \Leftarrow) }[/math] We have:

[math]\displaystyle{ \mathfrak{L} = \mathfrak{L}_I := \rho(I) = \left\{ (a_0,\ldots,a_{n-1}) \mid \sum_{i=0}^{n-1}a_ix^i \in I\right\} }[/math]

Let [math]\displaystyle{ (a_0,\ldots,a_{n-1}) }[/math] be an arbitrary element in [math]\displaystyle{ \mathfrak{L} }[/math]. Then, define [math]\displaystyle{ p(x) = \sum_{i=0}^{n-1}a_ix^i \in I }[/math]. But since [math]\displaystyle{ I }[/math] is an ideal, we have [math]\displaystyle{ xp(x) \in I }[/math]. Then, [math]\displaystyle{ \rho(xp(x)) \in \mathfrak{L}_I }[/math]. But, [math]\displaystyle{ \rho(xp(x)) = \operatorname{rot}(a_0,\ldots,a_{n-1}) \in \mathfrak{L}_I }[/math]. Hence, [math]\displaystyle{ \mathfrak{L} }[/math] is cyclic.

[math]\displaystyle{ \Rightarrow) }[/math]

Let [math]\displaystyle{ \mathfrak{L} \subset \Z^n }[/math] be a cyclic lattice. Hence [math]\displaystyle{ \forall (a_0,\ldots,a_{n-1}) \in \mathfrak{L}: \operatorname{rot}(a_0,\ldots,a_{n-1}) \in \mathfrak{L} }[/math].

Define the set of polynomials [math]\displaystyle{ I: = \left\{ \sum_{i=0}^{n-1}a_ix^i \mid (a_0,\ldots,a_{n-1}) \in \mathfrak{L} \right\} }[/math]:

  1. Since [math]\displaystyle{ \mathfrak{L} }[/math] a lattice and hence an additive subgroup of [math]\displaystyle{ \Z^n }[/math], [math]\displaystyle{ I \subset R }[/math] is an additive subgroup of [math]\displaystyle{ R }[/math].
  2. Since [math]\displaystyle{ \mathfrak{L} }[/math] is cyclic, [math]\displaystyle{ \forall p(x) \in I: xp(x) \in I }[/math].

Hence, [math]\displaystyle{ I \subset R }[/math] is an ideal, and consequently, [math]\displaystyle{ \mathfrak{L} = \mathfrak{L}_I }[/math].

Ideal lattices[4]

Let [math]\displaystyle{ f(x) \in \Z[x] }[/math] be a monic polynomial of degree [math]\displaystyle{ n }[/math]. For cryptographic applications, [math]\displaystyle{ f(x) }[/math] is usually selected to be irreducible. The ideal generated by [math]\displaystyle{ f(x) }[/math] is:

[math]\displaystyle{ (f(x)) := f(x) \cdot\Z[x] = \{ f(x)g(x): \forall g(x) \in \Z[x] \}. }[/math]

The quotient polynomial ring [math]\displaystyle{ R = \Z[x]/(f(x)) }[/math] partitions [math]\displaystyle{ \Z[x] }[/math] into equivalence classes of polynomials of degree at most [math]\displaystyle{ n-1 }[/math]:

[math]\displaystyle{ R = \Z[x]/(f(x)) = \left\{ \sum_{i=0}^{n-1}a_ix^i: a_i \in \Z \right\} }[/math]

where addition and multiplication are reduced modulo [math]\displaystyle{ f(x) }[/math].

Consider the embedding coefficient [math]\displaystyle{ \Z }[/math]-module isomorphism [math]\displaystyle{ \rho }[/math]. Then, each ideal in [math]\displaystyle{ R }[/math] defines a sublattice of [math]\displaystyle{ \Z^n }[/math] called ideal lattice.

Definition: [math]\displaystyle{ \mathfrak{L}_I }[/math], the lattice corresponding to an ideal [math]\displaystyle{ I }[/math], is called ideal lattice. More precisely, consider a quotient polynomial ring [math]\displaystyle{ R = \Z[x]/(p(x)) }[/math], where [math]\displaystyle{ (p(x)) }[/math] is the ideal generated by a degree [math]\displaystyle{ n }[/math] polynomial [math]\displaystyle{ p(x) \in \Z[x] }[/math]. [math]\displaystyle{ \mathfrak{L}_I }[/math], is a sublattice of [math]\displaystyle{ \Z^n }[/math], and is defined as:

[math]\displaystyle{ \mathfrak{L}_I := \rho(I) = \left\{ (a_0,\ldots,a_{n-1}) \mid \sum_{i=0}^{n-1}a_i x^i \in I \right\} \subset \Z^n. }[/math]

Remark:[5]

  1. It turns out that [math]\displaystyle{ \text{GapSVP} \gamma }[/math] for even small [math]\displaystyle{ \gamma = \operatorname{poly(n)} }[/math] is typically easy on ideal lattices. The intuition is that the algebraic symmetries causes the minimum distance of an ideal to lie within a narrow, easily computable range.
  2. By exploiting the provided algebraic symmetries in ideal lattices, one can convert a short nonzero vector into [math]\displaystyle{ n }[/math] linearly independent ones of (nearly) the same length. Therefore, on ideal lattices, [math]\displaystyle{ \mathrm{SVP} \gamma }[/math] and [math]\displaystyle{ \mathrm{SIVP} \gamma }[/math] are equivalent with a small loss.[6] Furthermore, even for quantum algorithms, [math]\displaystyle{ \mathrm{SVP} \gamma }[/math] and [math]\displaystyle{ \mathrm{SIVP} \gamma }[/math] are believed to be very hard in the worst-case scenario.

Short integer solution problem

The Short Integer Solution (SIS) problem is an average case problem that is used in lattice-based cryptography constructions. Lattice-based cryptography began in 1996 from a seminal work by Ajtai[1] who presented a family of one-way functions based on the SIS problem. He showed that it is secure in an average case if [math]\displaystyle{ \mathrm{SVP}_{\gamma} }[/math] (where [math]\displaystyle{ \gamma = n^c }[/math] for some constant [math]\displaystyle{ c\gt 0 }[/math]) is hard in a worst-case scenario. Along with applications in classical cryptography, the SIS problem and its variants are utilized in several post-quantum security schemes including CRYSTALS-Dilithium CRYSTALS-Dilithium and Falcon Falcon. [7][8]

SISq,n,m,β

Let [math]\displaystyle{ A \in \Z^{n\times m}_q }[/math] be an [math]\displaystyle{ n\times m }[/math] matrix with entries in [math]\displaystyle{ \Z_q }[/math] that consists of [math]\displaystyle{ m }[/math] uniformly random vectors [math]\displaystyle{ \boldsymbol{a_i} \in \Z^n_q }[/math]: [math]\displaystyle{ A = [\boldsymbol{a_1}|\cdots|\boldsymbol{a_m}] }[/math]. Find a nonzero vector [math]\displaystyle{ \boldsymbol{x} \in \Z^m }[/math] such that for some norm [math]\displaystyle{ \|\cdot\| }[/math]:

  • [math]\displaystyle{ 0 \lt \|\boldsymbol{x}\| \leq \beta }[/math],
  • [math]\displaystyle{ f_A(\boldsymbol{x}) := A\boldsymbol{x} = \boldsymbol{0} \in \Z^n_q }[/math].

A solution to SIS without the required constraint on the length of the solution ([math]\displaystyle{ \|\boldsymbol{x}\| \leq \beta }[/math]) is easy to compute by using Gaussian elimination technique. We also require [math]\displaystyle{ \beta \lt q }[/math], otherwise [math]\displaystyle{ \boldsymbol{x} = (q,0,\ldots,0) \in \Z^m }[/math] is a trivial solution.

In order to guarantee [math]\displaystyle{ f_A(\boldsymbol{x}) }[/math] has non-trivial, short solution, we require:

  • [math]\displaystyle{ \beta \geq \sqrt{n\log q} }[/math], and
  • [math]\displaystyle{ m \geq n\log q }[/math]

Theorem: For any [math]\displaystyle{ m = \operatorname{poly}(n) }[/math], any [math]\displaystyle{ \beta \gt 0 }[/math], and any sufficiently large [math]\displaystyle{ q \geq \beta n^c }[/math] (for any constant [math]\displaystyle{ c \gt 0 }[/math]), solving [math]\displaystyle{ \operatorname{SIS}_{n,m,q,\beta} }[/math] with nonnegligible probability is at least as hard as solving the [math]\displaystyle{ \operatorname{GapSVP}_\gamma }[/math] and [math]\displaystyle{ \operatorname{SIVP}_\gamma }[/math] for some [math]\displaystyle{ \gamma = \beta \cdot O(\sqrt{n}) }[/math] with a high probability in the worst-case scenario.

R-SISq,n,m,β

The SIS problem solved over an ideal ring is also called the Ring-SIS or R-SIS problem.[2][9] This problem considers a quotient polynomial ring [math]\displaystyle{ R_q = \Z_q[x]/(f(x)) }[/math] with [math]\displaystyle{ f(x) = x^n-1 }[/math] for some integer [math]\displaystyle{ n }[/math] and with some norm [math]\displaystyle{ \|\cdot\| }[/math]. Of particular interest are cases where there exists integer [math]\displaystyle{ k }[/math] such that [math]\displaystyle{ n=2^k }[/math] as this restricts the quotient to cyclotomic polynomials. [10]

We then define the problem as follows:

Select [math]\displaystyle{ m }[/math] independent uniformly random elements [math]\displaystyle{ a_i \in R_q }[/math]. Define vector [math]\displaystyle{ \vec{\boldsymbol{a}}:=(a_1,\ldots,a_m) \in R_q^m }[/math]. Find a nonzero vector [math]\displaystyle{ \vec{\boldsymbol{z}}:=(z_1,\ldots,z_m) \in R^m }[/math] such that:

  • [math]\displaystyle{ 0\lt \|\vec{\boldsymbol{z}}\| \leq \beta }[/math],
  • [math]\displaystyle{ f_{\vec{\boldsymbol{a}}}(\vec{\boldsymbol{z}}) := \vec{\boldsymbol{a}}^{T}.\vec{\boldsymbol{z}} = \sum_{i=1}^m a_i.p_i = 0 \in R_q }[/math].

Recall that to guarantee existence of a solution to SIS problem, we require [math]\displaystyle{ m \approx n\log q }[/math]. However, Ring-SIS problem provide us with more compactness and efficacy: to guarantee existence of a solution to Ring-SIS problem, we require [math]\displaystyle{ m \approx \log q }[/math] .

Definition: The nega-circulant matrix of [math]\displaystyle{ b }[/math] is defined as:

[math]\displaystyle{ \text{for} \quad b = \sum_{i=0}^{n-1}b_ix^i \in R, \quad \operatorname{rot}(b) := \begin{bmatrix} b_0 & -b_{n-1} & \ldots & -b_1 \\[0.3em] b_1 & b_0 & \ldots & -b_2 \\[0.3em] \vdots & \vdots & \ddots & \vdots \\[0.3em] b_{n-1} & b_{n-2} & \ldots & b_0 \end{bmatrix} }[/math]

When the quotient polynomial ring is [math]\displaystyle{ R = \Z[x]/(x^n+1) }[/math] for [math]\displaystyle{ n = 2^k }[/math], the ring multiplication [math]\displaystyle{ a_i.p_i }[/math] can be efficiently computed by first forming [math]\displaystyle{ \operatorname{rot}(a_i) }[/math], the nega-circulant matrix of [math]\displaystyle{ a_i }[/math], and then multiplying [math]\displaystyle{ \operatorname{rot}(a_i) }[/math] with [math]\displaystyle{ \rho(p_i(x)) \in Z^n }[/math], the embedding coefficient vector of [math]\displaystyle{ p_i }[/math] (or, alternatively with [math]\displaystyle{ \sigma(p_i(x)) \in Z^n }[/math], the canonical coefficient vector.

Moreover, R-SIS problem is a special case of SIS problem where the matrix [math]\displaystyle{ A }[/math] in the SIS problem is restricted to negacirculant blocks: [math]\displaystyle{ A = [\operatorname{rot}(a_1)|\cdots|\operatorname{rot}(a_m)] }[/math]. [10]

M-SISq,n,d,m,β

The SIS problem solved over a module lattice is also called the Module-SIS or M-SIS problem. Like R-SIS, this problem considers the quotient polynomial ring [math]\displaystyle{ R = \Z[x]/(f(x)) }[/math] and [math]\displaystyle{ R_q = \Z_q[x]/(f(x)) }[/math] for [math]\displaystyle{ f(x) = x^n-1 }[/math] with a special interest in cases where [math]\displaystyle{ n }[/math] is a power of 2. Then, let [math]\displaystyle{ M }[/math] be a module of rank [math]\displaystyle{ d }[/math] such that [math]\displaystyle{ M\subseteq R^d }[/math] and let [math]\displaystyle{ \|\cdot\| }[/math] be an arbitrary norm over [math]\displaystyle{ R_q^{m} }[/math].

We then define the problem as follows:

Select [math]\displaystyle{ m }[/math] independent uniformly random elements [math]\displaystyle{ a_i \in R_q^d }[/math]. Define vector [math]\displaystyle{ \vec{\boldsymbol{a}}:=(a_1,\ldots,a_m) \in R_q^{d\times m} }[/math]. Find a nonzero vector [math]\displaystyle{ \vec{\boldsymbol{z}}:=(z_1,\ldots,z_m) \in R^m }[/math] such that:

  • [math]\displaystyle{ 0 \lt \|\vec{\boldsymbol{z}}\| \leq \beta }[/math],
  • [math]\displaystyle{ f_{\vec{\boldsymbol{a}}}(\vec{\boldsymbol{z}}) := \vec{\boldsymbol{a}}^{T}.\vec{\boldsymbol{z}} = \sum_{i=1}^m a_i.p_i = 0 \in R_q^d }[/math].

While M-SIS is a less compact variant of SIS than R-SIS, the M-SIS problem is asymptotically at least as hard as R-SIS and therefore gives a tighter bound on the hardness assumption of SIS. This makes assuming the hardness of M-SIS a safer, but less efficient underlying assumption when compared to R-SIS. [10]

See also

References

  1. 1.0 1.1 Ajtai, Miklós. [Generating hard instances of lattice problems.] Proceedings of the twenty-eighth annual ACM symposium on Theory of computing. ACM, 1996.
  2. 2.0 2.1 Micciancio, Daniele. [Generalized compact knapsacks, cyclic lattices, and efficient one-way functions from worst-case complexity assumptions.] Foundations of Computer Science, 2002. Proceedings. The 43rd Annual IEEE Symposium on. IEEE, 2002.
  3. Fukshansky, Lenny, and Xun Sun. [On the geometry of cyclic lattices.] Discrete & Computational Geometry 52.2 (2014): 240–259.
  4. Craig Gentry. Fully Homomorphic Encryption Using Ideal Lattices. In the 41st ACM Symposium on Theory of Computing (STOC), 2009.
  5. Peikert, Chris. [A decade of lattice cryptography.] Cryptology ePrint Archive, Report 2015/939, 2015.
  6. Peikert, Chris, and Alon Rosen. [Efficient collision-resistant hashing from worst-case assumptions on cyclic lattices.] Theory of Cryptography. Springer Berlin Heidelberg, 2006. 145–166.
  7. Bai, Shi; Ducas, Léo; Kiltz, Eike; Lepoint, Tancrède; Lyubashevsky, Vadim; Schwabe, Peter; Seiler, Grego4; Stehlé, Damien (October 1, 2020). "CRYSTALS-Dilithium:Algorithm Specifications and Supporting Documentation". https://pq-crystals.org/dilithium/data/dilithium-specification-round3.pdf. 
  8. Fouque, Pierre-Alain; Hoffstein, Jeffrey; Kirchner, Paul; Lyubashevsky, Vadim; Pornin, Thomas; Prest, Thomas; Ricosset, Thomas; Seiler, Gregor et al. (October 1, 2020). "Falcon: Fast-Fourier Lattice-based Compact Signatures over NTRU". https://falcon-sign.info. 
  9. Lyubashevsky, Vadim, et al. [SWIFFT: A modest proposal for FFT hashing.] Fast Software Encryption. Springer Berlin Heidelberg, 2008.
  10. 10.0 10.1 10.2 Langlois, Adeline, and Damien Stehlé. [Worst-case to average-case reductions for module lattices.] Designs, Codes and Cryptography 75.3 (2015): 565–599.