MathLabs

算术与数论

二次剩余

在给定模数下是完全平方数的整数,通过二次互反律来研究。

直观直觉:哪些余数是完全平方数?

固定一个奇素数 pp,将余数 0,1,…,p−10,1,\dots,p-1 对 pp 逐一平方。由于 x2≡(p−x)2(modp)x^2 \equiv (p-x)^2 \pmod p,平方值成对重复出现,因此实际出现的非零余数只有大约一半。满足 gcd⁡(a,p)=1\gcd(a,p)=1 的数 aa,若同余式 x2≡a(modp)x^2 \equiv a \pmod{p} 有解,则称为模 pp 的“二次剩余”,否则称为“二次非剩余”。当 p=7p=7 时:将 1,2,3,4,5,61,2,3,4,5,6 平方得到 1,4,2,2,4,11,4,2,2,4,1,所以模 77 的二次剩余恰为 {1,2,4}\{1,2,4\},非剩余为 {3,5,6}\{3,5,6\}。

按一个小素数取模的平方映射有向图,二次剩余顶点被高亮显示。
模 m=11m = 11 的剩余弦:在 1010 个非零剩余中,恰好有一半(1,3,4,5,91, 3, 4, 5, 9)是二次剩余。

大学定义与勒让德符号

定义: 勒让德符号

对奇素数 pp 及满足 p∤ap \nmid a 的整数 aa,勒让德符号 (ap)\left(\dfrac{a}{p}\right) 在 aa 是模 pp 的二次剩余时取 +1+1,是二次非剩余时取 −1-1。按惯例当 p∣ap \mid a 时 (ap)=0\left(\dfrac{a}{p}\right)=0。勒让德符号对分子完全积性:(abp)=(ap)(bp)\left(\frac{ab}{p}\right) = \left(\frac{a}{p}\right)\left(\frac{b}{p}\right)。

(ap)≡ap−12(modp)(Euler’s criterion)\left(\dfrac{a}{p}\right) \equiv a^{\frac{p-1}{2}} \pmod{p}\qquad\text{(Euler’s criterion)}

欧拉判别法让我们无需做任何分解就能“计算” (ap)\left(\dfrac{a}{p}\right):把 aa 对 pp 取 p−12\frac{p-1}{2} 次方,再读出结果是 ±1\pm 1。它也立即给出 (−1p)=(−1)p−12\left(\dfrac{-1}{p}\right) = (-1)^{\frac{p-1}{2}},由此可知 −1-1 是二次剩余当且仅当 p≡1(mod4)p \equiv 1 \pmod 4。

(pq)(qp)=(−1)p−12⋅q−12(quadratic reciprocity)\left(\dfrac{p}{q}\right)\left(\dfrac{q}{p}\right) = (-1)^{\frac{p-1}{2}\cdot\frac{q-1}{2}}\qquad\text{(quadratic reciprocity)}
按 pp 的剩余类给出 (−1p)\left(\frac{-1}{p}\right) 与 (2p)\left(\frac{2}{p}\right) 的值
pp 的类(−1p)\left(\frac{-1}{p}\right)(2p)\left(\frac{2}{p}\right)示例
p≡1(mod4)p \equiv 1 \pmod 4+1+1取决于 p mod 8p \bmod 8p=13p=13
p≡3(mod4)p \equiv 3 \pmod 4−1-1取决于 p mod 8p \bmod 8p=7p=7
p≡1,7(mod8)p \equiv 1, 7 \pmod 8见上一行+1+1p=7p=7
p≡3,5(mod8)p \equiv 3, 5 \pmod 8见上一行−1-1p=13p=13

大学定理与证明

对奇素数 pp 及 gcd⁡(a,p)=1\gcd(a,p)=1,有 (ap)≡ap−12(modp)\left(\dfrac{a}{p}\right) \equiv a^{\frac{p-1}{2}} \pmod{p};等价地,aa 是模 pp 的二次剩余当且仅当 a(p−1)/2≡1(modp)a^{(p-1)/2} \equiv 1 \pmod p。

为什么成立?

它把一个存在性问题(是否存在 xx 使其平方为 aa)转化为一次模幂运算,正是这一点使勒让德符号既可计算又具有积性。

证明

由于 (Z/pZ)∗(\mathbb{Z}/p\mathbb{Z})^{*} 是阶为 p−1p-1 的循环群(素数模下的标准事实),固定一个原根 gg,写 a≡gk(modp)a \equiv g^k \pmod p,其中 0≤k≤p−20 \le k \le p-2。模 pp 的二次剩余恰好是 gg 的偶数次幂:若 a=g2ma=g^{2m},则 x=gmx=g^m 满足 x2≡ax^2\equiv a;反之每个平方 x2=(gj)2=g2jx^2=(g^j)^2=g^{2j} 都是偶数次幂。因此 aa 是二次剩余当且仅当 kk 为偶数。

现在计算 a(p−1)/2≡gk(p−1)/2(modp)a^{(p-1)/2} \equiv g^{k(p-1)/2} \pmod p。若 kk 为偶数,记 k=2mk=2m,由费马小定理得 gk(p−1)/2=gm(p−1)=(gp−1)m≡1m=1(modp)g^{k(p-1)/2}=g^{m(p-1)}=(g^{p-1})^m\equiv 1^m=1\pmod p,与剩余情形吻合。

若 kk 为奇数,则 gk(p−1)/2g^{k(p-1)/2} 是 gk(p−1)=(gp−1)k≡1(modp)g^{k(p-1)}=(g^{p-1})^k\equiv 1\pmod p 的平方根,故 a(p−1)/2≡±1a^{(p-1)/2}\equiv \pm 1。它不能等于 11:否则 gg 的阶必须整除 k(p−1)/2k(p-1)/2,但 ord⁡(g)=p−1\operatorname{ord}(g)=p-1 且 kk 为奇数意味着除非 (p−1)/2(p-1)/2 本身已是倍数,否则 k(p−1)/2k(p-1)/2 不是 p−1p-1 的倍数,这与 gg 是阶恰为 p−1p-1 的原根矛盾。所以恰好当 kk 为奇数即 aa 为非剩余时,a(p−1)/2≡−1(modp)a^{(p-1)/2}\equiv -1\pmod p。

综合两种情形:a(p−1)/2≡1(modp)a^{(p-1)/2}\equiv 1\pmod p 当且仅当 aa 是二次剩余,否则 ≡−1\equiv -1,这正是 (ap)≡ap−12(modp)\left(\dfrac{a}{p}\right) \equiv a^{\frac{p-1}{2}} \pmod{p} 的断言。

对不同的奇素数 pp 与 qq:(pq)(qp)=(−1)p−12⋅q−12\left(\dfrac{p}{q}\right)\left(\dfrac{q}{p}\right) = (-1)^{\frac{p-1}{2}\cdot\frac{q-1}{2}}。

为什么成立?

它说明 pp 是否为模 qq 的平方数,与 qq 是否为模 pp 的平方数,其实是“同一个”问题,差别仅在于一个只依赖 p,qp,q 模 44 的符号——把一个看似困难的双变量问题变成一条规则。

证明

我们借助高斯引理作为垫脚石:对奇素数 pp 及 gcd⁡(a,p)=1\gcd(a,p)=1,考察 a,2a,…,p−12aa, 2a, \dots, \frac{p-1}{2}a 模 pp 的最小正剩余,设其中超过 p/2p/2 的个数为 μ\mu。高斯引理断言 (ap)=(−1)μ\left(\dfrac{a}{p}\right)=(-1)^{\mu};这可由将每个这样的“大”剩余 rr 与 p−r≤p/2p-r\le p/2 配对、并以两种方式相乘这 p−12\frac{p-1}{2} 个数时追踪符号得到。

艾森斯坦的精细化把 μ\mu 表示为格点计数:当 a=qa=q 为奇数时,通过比较 ⌊kq/p⌋\lfloor kq/p\rfloor(小于 kqkq 的 pp 的倍数个数)与 kq mod pkq \bmod p 距 p/2p/2 的远近,可证明 μ≡∑k=1(p−1)/2⌊kqp⌋(mod2)\mu \equiv \sum_{k=1}^{(p-1)/2} \left\lfloor \frac{kq}{p} \right\rfloor \pmod 2。将同样的计数论证对称地应用,得到 (qp)=(−1)S(q,p)\left(\frac{q}{p}\right)=(-1)^{S(q,p)} 与 (pq)=(−1)S(p,q)\left(\frac{p}{q}\right)=(-1)^{S(p,q)},其中 S(q,p)=∑k=1(p−1)/2⌊kqp⌋S(q,p)=\sum_{k=1}^{(p-1)/2}\left\lfloor \frac{kq}{p}\right\rfloor,S(p,q)=∑k=1(q−1)/2⌊kpq⌋S(p,q)=\sum_{k=1}^{(q-1)/2}\left\lfloor \frac{kp}{q}\right\rfloor。

从几何上看,S(q,p)+S(p,q)S(q,p)+S(p,q) 统计满足 1≤x≤p−121\le x\le \frac{p-1}{2}、1≤y≤q−121\le y\le \frac{q-1}{2} 且严格位于直线 qx=pyqx=py 下方的格点 (x,y)(x,y)(共 S(q,p)S(q,p) 个),加上严格位于其上方的格点(由 p,qp,q 的对称角色,共 S(p,q)S(p,q) 个)。由于 gcd⁡(p,q)=1\gcd(p,q)=1 且 x<px<p,没有格点恰好落在该直线上,故这两个计数合起来恰好覆盖整个 p−12⋅q−12\frac{p-1}{2}\cdot\frac{q-1}{2} 个格点的矩形。

因此 S(q,p)+S(p,q)=p−12⋅q−12S(q,p)+S(p,q) = \frac{p-1}{2}\cdot\frac{q-1}{2},将两个勒让德符号公式相乘得到 (pq)(qp)=(−1)S(p,q)+S(q,p)=(−1)p−12⋅q−12\left(\frac{p}{q}\right)\left(\frac{q}{p}\right)=(-1)^{S(p,q)+S(q,p)}=(-1)^{\frac{p-1}{2}\cdot\frac{q-1}{2}},这正是 (pq)(qp)=(−1)p−12⋅q−12\left(\dfrac{p}{q}\right)\left(\dfrac{q}{p}\right) = (-1)^{\frac{p-1}{2}\cdot\frac{q-1}{2}}。

大学实际应用与典型例题

二次剩余不仅是数论中的趣味话题:它们是一种公钥加密方案(Goldwasser–Micali)、一种广泛使用的伪随机比特生成器(Blum–Blum–Shub)以及用于设计具有尖锐自相关峰值的扩频雷达和类GPS测距码的勒让德序列的基础。

例题: 用互反律判断二次剩余

1010 是模 1313 的二次剩余吗?

解答

写 10=2⋅510=2\cdot 5,由积性得 (1013)=(213)(513)\left(\frac{10}{13}\right)=\left(\frac{2}{13}\right)\left(\frac{5}{13}\right)。

对于 (213)\left(\frac{2}{13}\right):由于 13≡5(mod8)13 \equiv 5 \pmod 8,关于 22 的补充公式给出 (213)=−1\left(\frac{2}{13}\right)=-1。

对于 (513)\left(\frac{5}{13}\right):由于 5≡1(mod4)5\equiv 1\pmod 4,互反律给出 (513)(135)=+1\left(\frac{5}{13}\right)\left(\frac{13}{5}\right)=+1,故 (513)=(135)=(35)\left(\frac{5}{13}\right)=\left(\frac{13}{5}\right)=\left(\frac{3}{5}\right)(因为 13≡3(mod5)13\equiv 3\pmod 5)。模 55 的平方数为 1,41,4,而 33 不在其中,故 (35)=−1\left(\frac{3}{5}\right)=-1,从而 (513)=−1\left(\frac{5}{13}\right)=-1。

相乘:(1013)=(−1)(−1)=+1\left(\frac{10}{13}\right)=(-1)(-1)=+1,所以 1010 确实是模 1313 的二次剩余——事实上 62=36≡10(mod13)6^2=36\equiv 10\pmod{13}。

例题: 扩频测距码中的勒让德序列

GPS式测距系统需要在零位移处有尖锐自相关峰、其余位移处相关几乎为零的二值码。用模素数 1111 的勒让德符号构造一个长度为 1111 的勒让德序列,并定性检验其自相关性。

解答

模 1111 的二次剩余为 {1,3,4,5,9}\{1,3,4,5,9\}(1,…,51,\dots,5 的平方),于是定义:若 ii 是剩余则 si=+1s_i=+1,若 ii 是非剩余 {2,6,7,8,10}\{2,6,7,8,10\} 则 si=−1s_i=-1,按惯例 s0=−1s_0=-1,得到 i=0,…,10i=0,\dots,10 的 ±1\pm 1 序列 s=(−1,+1,−1,+1,+1,+1,−1,−1,−1,+1,−1)s=(-1,+1,-1,+1,+1,+1,-1,-1,-1,+1,-1)。

关键的结构性事实是:对非零位移 k≢0k\not\equiv 0,自相关 ∑isisi+k\sum_i s_i s_{i+k} 借助勒让德符号的积性可化归为一个与 ∑x(x(x+k)11)\sum_x \left(\frac{x(x+k)}{11}\right) 密切相关的和,而经典的特征和估计表明,对模一个奇素数的任意非零位移 kk,该和都等于 −1-1——远小于 k=0k=0 时的峰值 1010。

这种二值自相关(每个非零位移都取固定的旁瓣值)正是使勒让德序列适合接收机同步的性质:将接收信号与已知序列的每个循环位移相关时,只有在真正对齐处才会产生明显巨大的尖峰,这正是GPS接收机锁定卫星码时序的方式。

当 x2≡a(modp)x^2\equiv a\pmod p 无解时,(ap)\left(\frac{a}{p}\right) 称为什么?

根据欧拉判别法,(ap)\left(\frac{a}{p}\right) 模 pp 同余于 aa 的第几次幂?

奇素数 pp 模 44 属于哪一类时,−1-1 是二次剩余?

哪种密码学构件直接基于判定模合数的二次剩余性的困难性?

参考文献

  1. Wikipedia contributors (2024). Quadratic reciprocity
  2. Kenneth Ireland, Michael Rosen (1990). A Classical Introduction to Modern Number Theory