14 LWE 问题与带噪线性方程 整理日期:2026-10-16。
前置笔记:09 高斯函数与离散高斯分布 、12 q 元格与模线性方程 、13 SIS 问题与短整数关系 。
本文采用整数误差形式,向量按列排列,长度使用欧氏范数。
1. 从带噪方程到 LWE 样本
带误差学习问题 (Learning with Errors,LWE)研究从带有随机误差的模线性方程中恢复秘密,或识别这些方程的分布。
设 q ≥ 2 q\ge2 q ≥ 2 ,秘密为
s ∈ Z q n , s\in\mathbb Z_q^n, s ∈ Z q n ,
误差分布 χ \chi χ 定义在 Z \mathbb Z Z 上,通常集中在 0 0 0 附近。
固定秘密 s s s 后,一个 LWE 样本按以下方式生成:
a ← U ( Z q n ) , e ← χ , b = ⟨ a , s ⟩ + e m o d q . a\leftarrow U(\mathbb Z_q^n),
\qquad
e\leftarrow\chi,
\qquad
\boxed{b=\langle a,s\rangle+e\bmod q.} a ← U ( Z q n ) , e ← χ , b = ⟨ a , s ⟩ + e mod q .
公开的是 ( a , b ) (a,b) ( a , b ) ,秘密 s s s 和误差 e e e 不公开。部分文献直接将误差分布定义在 Z q \mathbb Z_q Z q 上,对应这里的整数误差取模后的分布。样本定义见 Regev《The Learning with Errors Problem》,CCC 2010,§1,PDF 第 2 页 。
将 m m m 个样本合并,记
A = ( a 1 , … , a m ) ∈ Z q n × m , e = ( e 1 , … , e m ) T , A=(a_1,\ldots,a_m)\in\mathbb Z_q^{n\times m},
\qquad
e=(e_1,\ldots,e_m)^{\mathsf T}, A = ( a 1 , … , a m ) ∈ Z q n × m , e = ( e 1 , … , e m ) T ,
便得到
b = A T s + e m o d q . \boxed{
b=A^{\mathsf T}s+e\bmod q.
} b = A T s + e mod q .
这里 n n n 是秘密的维数,m m m 是样本数。所有样本共享同一个秘密;各个 a i a_i a i 和 e i e_i e i 独立抽取。
本文采用均匀秘密版本:每个实例开始时抽取一次 s ← U ( Z q n ) s\leftarrow U(\mathbb Z_q^n) s ← U ( Z q n ) ,随后固定它生成全部样本。矩阵形式及随机实例的构造见 Micciancio《Random Lattices and Lattice-Based Cryptography》,Fall 2019,PDF 第 1 页,Definition 2 及其后的分布说明 。
2. 噪声的作用与尺度
如果 e = 0 e=0 e = 0 ,且 q q q 为素数、A A A 满行秩,那么
A T s = b A^{\mathsf T}s=b A T s = b
可以通过有限域上的线性消元唯一求解。
加入未知误差后,对方程作线性组合时,误差也随之组合。模 q q q 上的可逆变换一般不保持整数误差的短性,因此直接消元不再自动恢复秘密。
“小误差”描述的是取模前的整数大小。例如,模 q q q 的剩余类 q − 1 q-1 q − 1 可以来自误差 − 1 -1 − 1 。
一个常用的高斯误差构造是
Z ∼ N ( 0 , r 2 2 π ) , e = ⌊ Z ⌉ , Z\sim\mathcal N\left(0,\frac{r^2}{2\pi}\right),
\qquad
e=\lfloor Z\rceil, Z ∼ N ( 0 , 2 π r 2 ) , e = ⌊ Z ⌉ ,
其中 ⌊ Z ⌉ \lfloor Z\rceil ⌊ Z ⌉ 表示取最近整数。这里 r r r 沿用第 09 篇的高斯尺度约定,取整前的标准差为
σ = r 2 π . \sigma=\frac{r}{\sqrt{2\pi}}. σ = 2 π r .
常记 r = α q r=\alpha q r = α q ,用 α \alpha α 表示高斯尺度相对于模数的大小。这与 Regev 原始论文的周期化、离散化高斯约定相衔接。见 Regev《On Lattices, Learning with Errors, Random Linear Codes, and Cryptography》,2009-05-02 作者稿,§2,PDF 第 14–15 页,式 (7)–(8) 。
需要保留两个区分:
连续高斯取整与第 09 篇的离散高斯 D Z , r D_{\mathbb Z,r} D Z , r 通常不是同一个分布。
高斯误差没有固定的绝对上界;对误差长度的控制通常是高概率保证。
因此,LWE 的参数需要说明 n , m , q n,m,q n , m , q 和完整的误差分布 χ \chi χ ,不能只写“加入少量噪声”。
3. 搜索 LWE 与判定 LWE
搜索 LWE:恢复秘密。
给定按第一节生成的 ( A , b ) (A,b) ( A , b ) ,目标是恢复生成这些样本时使用的 s s s 。相应参数可记为
L W E n , m , q , χ . \mathrm{LWE}_{n,m,q,\chi}. LWE n , m , q , χ .
对选定的渐近参数族,搜索 LWE 困难假设要求:任意概率多项式时间算法恢复 s s s 的成功概率都是可忽略的。概率计入 A , s , e A,s,e A , s , e 的抽样和算法自身的随机性。
判定 LWE:区分两类分布。
考虑
D L W E = ( A , A T s + e m o d q ) \mathcal D_{\mathrm{LWE}}
=
\left(A,A^{\mathsf T}s+e\bmod q\right) D LWE = ( A , A T s + e mod q )
与
D U = ( A , u ) , u ← U ( Z q m ) . \mathcal D_{\mathrm U}
=
(A,u),
\qquad
u\leftarrow U(\mathbb Z_q^m). D U = ( A , u ) , u ← U ( Z q m ) .
两类分布中的 A A A 都均匀随机;第二类中的 u u u 与 A A A 独立。
判定任务是识别输入来自哪一类分布。区分器 D D D 的优势定义为
Adv D = ∣ Pr [ D ( D L W E ) = 1 ] − Pr [ D ( D U ) = 1 ] ∣ . \boxed{
\operatorname{Adv}_D
=
\left|
\Pr[D(\mathcal D_{\mathrm{LWE}})=1]
-
\Pr[D(\mathcal D_{\mathrm U})=1]
\right|.
} Adv D = ∣ Pr [ D ( D LWE ) = 1 ] − Pr [ D ( D U ) = 1 ] ∣ .
判定 LWE 困难假设要求所有概率多项式时间区分器的优势都可忽略。这里的要求是计算上难以区分 ,并不要求两类分布在统计上接近。
在素数模数 q ≤ poly ( n ) q\le\operatorname{poly}(n) q ≤ poly ( n ) 、允许使用多项式数量样本的标准设置下,可以利用判定算法构造搜索算法,即
搜索 LWE ≤ p o l y 判定 LWE . \text{搜索 LWE}\ \le_{\mathrm{poly}}\ \text{判定 LWE}. 搜索 LWE ≤ poly 判定 LWE .
这一归约有明确的模数、样本和成功概率条件,不能理解为任意固定样本数下的无条件等价。见 Regev《The Learning with Errors Problem》,§3,PDF 第 12–13 页,Lemma 3.1–3.2 。
4. 格解释与有界距离解码
沿用第十二篇的模像格
Λ q ( A ) = { v ∈ Z m : v ≡ A T t ( m o d q ) ,某个 t ∈ Z q n } . \Lambda_q(A)
=
\left\{
v\in\mathbb Z^m:
v\equiv A^{\mathsf T}t\pmod q
\text{,某个 }t\in\mathbb Z_q^n
\right\}. Λ q ( A ) = { v ∈ Z m : v ≡ A T t ( mod q ) ,某个 t ∈ Z q n } .
取观测向量 b b b 的整数代表 b ~ \widetilde b b 。由 LWE 方程直接得到
v = b ~ − e ∈ Λ q ( A ) , v=\widetilde b-e\in\Lambda_q(A), v = b − e ∈ Λ q ( A ) ,
因此
dist ( b ~ , Λ q ( A ) ) ≤ ∥ e ∥ 2 . \boxed{
\operatorname{dist}\bigl(\widetilde b,\Lambda_q(A)\bigr)
\le\|e\|_2.
} dist ( b , Λ q ( A ) ) ≤ ∥ e ∥ 2 .
也就是说,观测点位于模像格中某个格点的附近。这一格解释见 Micciancio,§2,PDF 第 4 页 。
若进一步满足
∥ e ∥ 2 < 1 2 λ 1 ( Λ q ( A ) ) , \boxed{
\|e\|_2
<
\frac12\lambda_1\bigl(\Lambda_q(A)\bigr),
} ∥ e ∥ 2 < 2 1 λ 1 ( Λ q ( A ) ) ,
则 v v v 是 b ~ \widetilde b b 的唯一最近格点。因为对任何其他格点 w ≠ v w\ne v w = v ,有
∥ b ~ − w ∥ 2 ≥ ∥ w − v ∥ 2 − ∥ e ∥ 2 ≥ λ 1 ( Λ q ( A ) ) − ∥ e ∥ 2 > ∥ e ∥ 2 . \|\widetilde b-w\|_2
\ge
\|w-v\|_2-\|e\|_2
\ge
\lambda_1\bigl(\Lambda_q(A)\bigr)-\|e\|_2
>
\|e\|_2. ∥ b − w ∥ 2 ≥ ∥ w − v ∥ 2 − ∥ e ∥ 2 ≥ λ 1 ( Λ q ( A ) ) − ∥ e ∥ 2 > ∥ e ∥ 2 .
这给出有界距离解码 (Bounded Distance Decoding,BDD)的视角:在目标点足够靠近格的承诺下,恢复其唯一最近格点。
恢复 v v v 后,若 q q q 为素数且 A A A 满行秩,再解
A T s ≡ v ( m o d q ) A^{\mathsf T}s\equiv v\pmod q A T s ≡ v ( mod q )
即可唯一恢复秘密。一般 LWE 定义本身并不自动保证上述距离条件;使用 BDD 解释时需要核对参数。
5. 与 SIS 的联系
对同一个矩阵 A A A :
SIS 在模核格 Λ q ⊥ ( A ) \Lambda_q^\perp(A) Λ q ⊥ ( A ) 中寻找短的非零向量。
搜索 LWE 从模像格 Λ q ( A ) \Lambda_q(A) Λ q ( A ) 附近的观测中恢复秘密。
第十二篇已经给出两类格的缩放对偶关系。还可以直接从方程看出它们之间的联系。
若找到一个 SIS 关系
A x ≡ 0 ( m o d q ) , Ax\equiv0\pmod q, A x ≡ 0 ( mod q ) ,
则对 LWE 样本有
x T b ≡ x T A T s + x T e = ( A x ) T s + x T e ≡ x T e ( m o d q ) . \begin{aligned}
x^{\mathsf T}b
&\equiv x^{\mathsf T}A^{\mathsf T}s+x^{\mathsf T}e\\
&=(Ax)^{\mathsf T}s+x^{\mathsf T}e\\
&\equiv x^{\mathsf T}e
\pmod q.
\end{aligned} x T b ≡ x T A T s + x T e = ( A x ) T s + x T e ≡ x T e ( mod q ) .
因此,模核中的关系可以消去秘密项,留下误差的整数线性组合。
这个组合的取模前大小满足
∣ x T e ∣ ≤ ∥ x ∥ 2 ∥ e ∥ 2 . |x^{\mathsf T}e|
\le
\|x\|_2\|e\|_2. ∣ x T e ∣ ≤ ∥ x ∥ 2 ∥ e ∥ 2 .
短关系有助于控制组合后的噪声。但能否据此区分 LWE 样本,还取决于关系长度、误差分布与模数,不能仅凭存在一个 SIS 解就得出结论。
6. 困难性来源与后续阅读
Regev 的原始结果在特定高斯误差和参数条件下,建立了从最坏情况近似格问题到随机 LWE 的量子归约 。
其含义是:若能高效求解相应的 LWE,就能通过量子算法高效解决对应近似因子的最坏情况 GapSVP 和 SIVP。这里的量子性属于归约过程,LWE 样本本身仍由经典运算生成。见 Regev,2009-05-02 作者稿,PDF 第 2 页,Theorem 1.1 。
本篇先建立样本分布与计算任务。后续结合离散高斯和平滑参数学习归约,再以 Regev 公钥加密理解:合法解密如何控制误差,安全性证明如何使用判定 LWE。
引用此笔记