← 返回 资料库 › 算术与数论 › 同余 算术与数论
二次剩余 在给定模数下是完全平方数的整数,通过二次互反律来研究。
直观 直觉:哪些余数是完全平方数? 固定一个奇素数 p p p ,将余数 0 , 1 , … , p − 1 0,1,\dots,p-1 0 , 1 , … , p − 1 对 p p p 逐一平方。由于 x 2 ≡ ( p − x ) 2 ( m o d p ) x^2 \equiv (p-x)^2 \pmod p x 2 ≡ ( p − x ) 2 ( mod p ) ,平方值成对重复出现,因此实际出现的非零余数只有大约一半。满足 gcd ( a , p ) = 1 \gcd(a,p)=1 g cd( a , p ) = 1 的数 a a a ,若同余式 x 2 ≡ a ( m o d p ) x^2 \equiv a \pmod{p} x 2 ≡ a ( mod p ) 有解,则称为模 p p p 的“二次剩余”,否则称为“二次非剩余”。当 p = 7 p=7 p = 7 时:将 1 , 2 , 3 , 4 , 5 , 6 1,2,3,4,5,6 1 , 2 , 3 , 4 , 5 , 6 平方得到 1 , 4 , 2 , 2 , 4 , 1 1,4,2,2,4,1 1 , 4 , 2 , 2 , 4 , 1 ,所以模 7 7 7 的二次剩余恰为 { 1 , 2 , 4 } \{1,2,4\} { 1 , 2 , 4 } ,非剩余为 { 3 , 5 , 6 } \{3,5,6\} { 3 , 5 , 6 } 。
模 m = 11 m = 11 m = 11 的剩余弦:在 10 10 10 个非零剩余中,恰好有一半(1 , 3 , 4 , 5 , 9 1, 3, 4, 5, 9 1 , 3 , 4 , 5 , 9 )是二次剩余。 大学 定义与勒让德符号 定义: 勒让德符号
对奇素数 p p p 及满足 p ∤ a p \nmid a p ∤ a 的整数 a a a ,勒让德符号 ( a p ) \left(\dfrac{a}{p}\right) ( p a ) 在 a a a 是模 p p p 的二次剩余时取 + 1 +1 + 1 ,是二次非剩余时取 − 1 -1 − 1 。按惯例当 p ∣ a p \mid a p ∣ a 时 ( a p ) = 0 \left(\dfrac{a}{p}\right)=0 ( p a ) = 0 。勒让德符号对分子完全积性:( a b p ) = ( a p ) ( b p ) \left(\frac{ab}{p}\right) = \left(\frac{a}{p}\right)\left(\frac{b}{p}\right) ( p ab ) = ( p a ) ( p b ) 。
( a p ) ≡ a p − 1 2 ( m o d p ) (Euler’s criterion) \left(\dfrac{a}{p}\right) \equiv a^{\frac{p-1}{2}} \pmod{p}\qquad\text{(Euler’s criterion)} ( p a ) ≡ a 2 p − 1 ( mod p ) (Euler’s criterion) 欧拉判别法让我们无需做任何分解就能“计算” ( a p ) \left(\dfrac{a}{p}\right) ( p a ) :把 a a a 对 p p p 取 p − 1 2 \frac{p-1}{2} 2 p − 1 次方,再读出结果是 ± 1 \pm 1 ± 1 。它也立即给出 ( − 1 p ) = ( − 1 ) p − 1 2 \left(\dfrac{-1}{p}\right) = (-1)^{\frac{p-1}{2}} ( p − 1 ) = ( − 1 ) 2 p − 1 ,由此可知 − 1 -1 − 1 是二次剩余当且仅当 p ≡ 1 ( m o d 4 ) p \equiv 1 \pmod 4 p ≡ 1 ( mod 4 ) 。
( p q ) ( q p ) = ( − 1 ) p − 1 2 ⋅ q − 1 2 (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)} ( q p ) ( p q ) = ( − 1 ) 2 p − 1 ⋅ 2 q − 1 (quadratic reciprocity) 按 p p p 的剩余类给出 ( − 1 p ) \left(\frac{-1}{p}\right) ( p − 1 ) 与 ( 2 p ) \left(\frac{2}{p}\right) ( p 2 ) 的值 p p p 的类( − 1 p ) \left(\frac{-1}{p}\right) ( p − 1 ) ( 2 p ) \left(\frac{2}{p}\right) ( p 2 ) 示例 p ≡ 1 ( m o d 4 ) p \equiv 1 \pmod 4 p ≡ 1 ( mod 4 ) + 1 +1 + 1 取决于 p m o d 8 p \bmod 8 p mod 8 p = 13 p=13 p = 13 p ≡ 3 ( m o d 4 ) p \equiv 3 \pmod 4 p ≡ 3 ( mod 4 ) − 1 -1 − 1 取决于 p m o d 8 p \bmod 8 p mod 8 p = 7 p=7 p = 7 p ≡ 1 , 7 ( m o d 8 ) p \equiv 1, 7 \pmod 8 p ≡ 1 , 7 ( mod 8 ) 见上一行 + 1 +1 + 1 p = 7 p=7 p = 7 p ≡ 3 , 5 ( m o d 8 ) p \equiv 3, 5 \pmod 8 p ≡ 3 , 5 ( mod 8 ) 见上一行 − 1 -1 − 1 p = 13 p=13 p = 13
大学 定理与证明 对奇素数 p p p 及 gcd ( a , p ) = 1 \gcd(a,p)=1 g cd( a , p ) = 1 ,有 ( a p ) ≡ a p − 1 2 ( m o d p ) \left(\dfrac{a}{p}\right) \equiv a^{\frac{p-1}{2}} \pmod{p} ( p a ) ≡ a 2 p − 1 ( mod p ) ;等价地,a a a 是模 p p p 的二次剩余当且仅当 a ( p − 1 ) / 2 ≡ 1 ( m o d p ) a^{(p-1)/2} \equiv 1 \pmod p a ( p − 1 ) /2 ≡ 1 ( mod p ) 。
为什么成立? 它把一个存在性问题(是否存在 x x x 使其平方为 a a a )转化为一次模幂运算,正是这一点使勒让德符号既可计算又具有积性。
证明 由于 ( Z / p Z ) ∗ (\mathbb{Z}/p\mathbb{Z})^{*} ( Z / p Z ) ∗ 是阶为 p − 1 p-1 p − 1 的循环群(素数模下的标准事实),固定一个原根 g g g ,写 a ≡ g k ( m o d p ) a \equiv g^k \pmod p a ≡ g k ( mod p ) ,其中 0 ≤ k ≤ p − 2 0 \le k \le p-2 0 ≤ k ≤ p − 2 。模 p p p 的二次剩余恰好是 g g g 的偶数次幂:若 a = g 2 m a=g^{2m} a = g 2 m ,则 x = g m x=g^m x = g m 满足 x 2 ≡ a x^2\equiv a x 2 ≡ a ;反之每个平方 x 2 = ( g j ) 2 = g 2 j x^2=(g^j)^2=g^{2j} x 2 = ( g j ) 2 = g 2 j 都是偶数次幂。因此 a a a 是二次剩余当且仅当 k k k 为偶数。
现在计算 a ( p − 1 ) / 2 ≡ g k ( p − 1 ) / 2 ( m o d p ) a^{(p-1)/2} \equiv g^{k(p-1)/2} \pmod p a ( p − 1 ) /2 ≡ g k ( p − 1 ) /2 ( mod p ) 。若 k k k 为偶数,记 k = 2 m k=2m k = 2 m ,由费马小定理得 g k ( p − 1 ) / 2 = g m ( p − 1 ) = ( g p − 1 ) m ≡ 1 m = 1 ( m o d p ) g^{k(p-1)/2}=g^{m(p-1)}=(g^{p-1})^m\equiv 1^m=1\pmod p g k ( p − 1 ) /2 = g m ( p − 1 ) = ( g p − 1 ) m ≡ 1 m = 1 ( mod p ) ,与剩余情形吻合。
若 k k k 为奇数,则 g k ( p − 1 ) / 2 g^{k(p-1)/2} g k ( p − 1 ) /2 是 g k ( p − 1 ) = ( g p − 1 ) k ≡ 1 ( m o d p ) g^{k(p-1)}=(g^{p-1})^k\equiv 1\pmod p g k ( p − 1 ) = ( g p − 1 ) k ≡ 1 ( mod p ) 的平方根,故 a ( p − 1 ) / 2 ≡ ± 1 a^{(p-1)/2}\equiv \pm 1 a ( p − 1 ) /2 ≡ ± 1 。它不能等于 1 1 1 :否则 g g g 的阶必须整除 k ( p − 1 ) / 2 k(p-1)/2 k ( p − 1 ) /2 ,但 ord ( g ) = p − 1 \operatorname{ord}(g)=p-1 ord ( g ) = p − 1 且 k k k 为奇数意味着除非 ( p − 1 ) / 2 (p-1)/2 ( p − 1 ) /2 本身已是倍数,否则 k ( p − 1 ) / 2 k(p-1)/2 k ( p − 1 ) /2 不是 p − 1 p-1 p − 1 的倍数,这与 g g g 是阶恰为 p − 1 p-1 p − 1 的原根矛盾。所以恰好当 k k k 为奇数即 a a a 为非剩余时,a ( p − 1 ) / 2 ≡ − 1 ( m o d p ) a^{(p-1)/2}\equiv -1\pmod p a ( p − 1 ) /2 ≡ − 1 ( mod p ) 。
综合两种情形:a ( p − 1 ) / 2 ≡ 1 ( m o d p ) a^{(p-1)/2}\equiv 1\pmod p a ( p − 1 ) /2 ≡ 1 ( mod p ) 当且仅当 a a a 是二次剩余,否则 ≡ − 1 \equiv -1 ≡ − 1 ,这正是 ( a p ) ≡ a p − 1 2 ( m o d p ) \left(\dfrac{a}{p}\right) \equiv a^{\frac{p-1}{2}} \pmod{p} ( p a ) ≡ a 2 p − 1 ( mod p ) 的断言。
对不同的奇素数 p p p 与 q q q :( p q ) ( q p ) = ( − 1 ) p − 1 2 ⋅ q − 1 2 \left(\dfrac{p}{q}\right)\left(\dfrac{q}{p}\right) = (-1)^{\frac{p-1}{2}\cdot\frac{q-1}{2}} ( q p ) ( p q ) = ( − 1 ) 2 p − 1 ⋅ 2 q − 1 。
为什么成立? 它说明 p p p 是否为模 q q q 的平方数,与 q q q 是否为模 p p p 的平方数,其实是“同一个”问题,差别仅在于一个只依赖 p , q p,q p , q 模 4 4 4 的符号——把一个看似困难的双变量问题变成一条规则。
证明 我们借助高斯引理作为垫脚石:对奇素数 p p p 及 gcd ( a , p ) = 1 \gcd(a,p)=1 g cd( a , p ) = 1 ,考察 a , 2 a , … , p − 1 2 a a, 2a, \dots, \frac{p-1}{2}a a , 2 a , … , 2 p − 1 a 模 p p p 的最小正剩余,设其中超过 p / 2 p/2 p /2 的个数为 μ \mu μ 。高斯引理断言 ( a p ) = ( − 1 ) μ \left(\dfrac{a}{p}\right)=(-1)^{\mu} ( p a ) = ( − 1 ) μ ;这可由将每个这样的“大”剩余 r r r 与 p − r ≤ p / 2 p-r\le p/2 p − r ≤ p /2 配对、并以两种方式相乘这 p − 1 2 \frac{p-1}{2} 2 p − 1 个数时追踪符号得到。
艾森斯坦的精细化把 μ \mu μ 表示为格点计数:当 a = q a=q a = q 为奇数时,通过比较 ⌊ k q / p ⌋ \lfloor kq/p\rfloor ⌊ k q / p ⌋ (小于 k q kq k q 的 p p p 的倍数个数)与 k q m o d p kq \bmod p k q mod p 距 p / 2 p/2 p /2 的远近,可证明 μ ≡ ∑ k = 1 ( p − 1 ) / 2 ⌊ k q p ⌋ ( m o d 2 ) \mu \equiv \sum_{k=1}^{(p-1)/2} \left\lfloor \frac{kq}{p} \right\rfloor \pmod 2 μ ≡ ∑ k = 1 ( p − 1 ) /2 ⌊ p k q ⌋ ( mod 2 ) 。将同样的计数论证对称地应用,得到 ( q p ) = ( − 1 ) S ( q , p ) \left(\frac{q}{p}\right)=(-1)^{S(q,p)} ( p q ) = ( − 1 ) S ( q , p ) 与 ( p q ) = ( − 1 ) S ( p , q ) \left(\frac{p}{q}\right)=(-1)^{S(p,q)} ( q p ) = ( − 1 ) S ( p , q ) ,其中 S ( q , p ) = ∑ k = 1 ( p − 1 ) / 2 ⌊ k q p ⌋ S(q,p)=\sum_{k=1}^{(p-1)/2}\left\lfloor \frac{kq}{p}\right\rfloor S ( q , p ) = ∑ k = 1 ( p − 1 ) /2 ⌊ p k q ⌋ ,S ( p , q ) = ∑ k = 1 ( q − 1 ) / 2 ⌊ k p q ⌋ S(p,q)=\sum_{k=1}^{(q-1)/2}\left\lfloor \frac{kp}{q}\right\rfloor S ( p , q ) = ∑ k = 1 ( q − 1 ) /2 ⌊ q k p ⌋ 。
从几何上看,S ( q , p ) + S ( p , q ) S(q,p)+S(p,q) S ( q , p ) + S ( p , q ) 统计满足 1 ≤ x ≤ p − 1 2 1\le x\le \frac{p-1}{2} 1 ≤ x ≤ 2 p − 1 、1 ≤ y ≤ q − 1 2 1\le y\le \frac{q-1}{2} 1 ≤ y ≤ 2 q − 1 且严格位于直线 q x = p y qx=py q x = p y 下方的格点 ( x , y ) (x,y) ( x , y ) (共 S ( q , p ) S(q,p) S ( q , p ) 个),加上严格位于其上方的格点(由 p , q p,q p , q 的对称角色,共 S ( p , q ) S(p,q) S ( p , q ) 个)。由于 gcd ( p , q ) = 1 \gcd(p,q)=1 g cd( p , q ) = 1 且 x < p x<p x < p ,没有格点恰好落在该直线上,故这两个计数合起来恰好覆盖整个 p − 1 2 ⋅ q − 1 2 \frac{p-1}{2}\cdot\frac{q-1}{2} 2 p − 1 ⋅ 2 q − 1 个格点的矩形。
因此 S ( q , p ) + S ( p , q ) = p − 1 2 ⋅ q − 1 2 S(q,p)+S(p,q) = \frac{p-1}{2}\cdot\frac{q-1}{2} S ( q , p ) + S ( p , q ) = 2 p − 1 ⋅ 2 q − 1 ,将两个勒让德符号公式相乘得到 ( p q ) ( q p ) = ( − 1 ) S ( p , q ) + S ( q , p ) = ( − 1 ) p − 1 2 ⋅ q − 1 2 \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}} ( q p ) ( p q ) = ( − 1 ) S ( p , q ) + S ( q , p ) = ( − 1 ) 2 p − 1 ⋅ 2 q − 1 ,这正是 ( p q ) ( q p ) = ( − 1 ) p − 1 2 ⋅ q − 1 2 \left(\dfrac{p}{q}\right)\left(\dfrac{q}{p}\right) = (-1)^{\frac{p-1}{2}\cdot\frac{q-1}{2}} ( q p ) ( p q ) = ( − 1 ) 2 p − 1 ⋅ 2 q − 1 。
大学 实际应用与典型例题 二次剩余不仅是数论中的趣味话题:它们是一种公钥加密方案(Goldwasser–Micali)、一种广泛使用的伪随机比特生成器(Blum–Blum–Shub)以及用于设计具有尖锐自相关峰值的扩频雷达和类GPS测距码的勒让德序列的基础。
例题: 用互反律判断二次剩余
10 10 10 是模 13 13 13 的二次剩余吗?
解答 写 10 = 2 ⋅ 5 10=2\cdot 5 10 = 2 ⋅ 5 ,由积性得 ( 10 13 ) = ( 2 13 ) ( 5 13 ) \left(\frac{10}{13}\right)=\left(\frac{2}{13}\right)\left(\frac{5}{13}\right) ( 13 10 ) = ( 13 2 ) ( 13 5 ) 。
对于 ( 2 13 ) \left(\frac{2}{13}\right) ( 13 2 ) :由于 13 ≡ 5 ( m o d 8 ) 13 \equiv 5 \pmod 8 13 ≡ 5 ( mod 8 ) ,关于 2 2 2 的补充公式给出 ( 2 13 ) = − 1 \left(\frac{2}{13}\right)=-1 ( 13 2 ) = − 1 。
对于 ( 5 13 ) \left(\frac{5}{13}\right) ( 13 5 ) :由于 5 ≡ 1 ( m o d 4 ) 5\equiv 1\pmod 4 5 ≡ 1 ( mod 4 ) ,互反律给出 ( 5 13 ) ( 13 5 ) = + 1 \left(\frac{5}{13}\right)\left(\frac{13}{5}\right)=+1 ( 13 5 ) ( 5 13 ) = + 1 ,故 ( 5 13 ) = ( 13 5 ) = ( 3 5 ) \left(\frac{5}{13}\right)=\left(\frac{13}{5}\right)=\left(\frac{3}{5}\right) ( 13 5 ) = ( 5 13 ) = ( 5 3 ) (因为 13 ≡ 3 ( m o d 5 ) 13\equiv 3\pmod 5 13 ≡ 3 ( mod 5 ) )。模 5 5 5 的平方数为 1 , 4 1,4 1 , 4 ,而 3 3 3 不在其中,故 ( 3 5 ) = − 1 \left(\frac{3}{5}\right)=-1 ( 5 3 ) = − 1 ,从而 ( 5 13 ) = − 1 \left(\frac{5}{13}\right)=-1 ( 13 5 ) = − 1 。
相乘:( 10 13 ) = ( − 1 ) ( − 1 ) = + 1 \left(\frac{10}{13}\right)=(-1)(-1)=+1 ( 13 10 ) = ( − 1 ) ( − 1 ) = + 1 ,所以 10 10 10 确实是模 13 13 13 的二次剩余——事实上 6 2 = 36 ≡ 10 ( m o d 13 ) 6^2=36\equiv 10\pmod{13} 6 2 = 36 ≡ 10 ( mod 13 ) 。
例题: 扩频测距码中的勒让德序列
GPS式测距系统需要在零位移处有尖锐自相关峰、其余位移处相关几乎为零的二值码。用模素数 11 11 11 的勒让德符号构造一个长度为 11 11 11 的勒让德序列,并定性检验其自相关性。
解答 模 11 11 11 的二次剩余为 { 1 , 3 , 4 , 5 , 9 } \{1,3,4,5,9\} { 1 , 3 , 4 , 5 , 9 } (1 , … , 5 1,\dots,5 1 , … , 5 的平方),于是定义:若 i i i 是剩余则 s i = + 1 s_i=+1 s i = + 1 ,若 i i i 是非剩余 { 2 , 6 , 7 , 8 , 10 } \{2,6,7,8,10\} { 2 , 6 , 7 , 8 , 10 } 则 s i = − 1 s_i=-1 s i = − 1 ,按惯例 s 0 = − 1 s_0=-1 s 0 = − 1 ,得到 i = 0 , … , 10 i=0,\dots,10 i = 0 , … , 10 的 ± 1 \pm 1 ± 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) s = ( − 1 , + 1 , − 1 , + 1 , + 1 , + 1 , − 1 , − 1 , − 1 , + 1 , − 1 ) 。
关键的结构性事实是:对非零位移 k ≢ 0 k\not\equiv 0 k ≡ 0 ,自相关 ∑ i s i s i + k \sum_i s_i s_{i+k} ∑ i s i s i + k 借助勒让德符号的积性可化归为一个与 ∑ x ( x ( x + k ) 11 ) \sum_x \left(\frac{x(x+k)}{11}\right) ∑ x ( 11 x ( x + k ) ) 密切相关的和,而经典的特征和估计表明,对模一个奇素数的任意非零位移 k k k ,该和都等于 − 1 -1 − 1 ——远小于 k = 0 k=0 k = 0 时的峰值 10 10 10 。
这种二值自相关(每个非零位移都取固定的旁瓣值)正是使勒让德序列适合接收机同步的性质:将接收信号与已知序列的每个循环位移相关时,只有在真正对齐处才会产生明显巨大的尖峰,这正是GPS接收机锁定卫星码时序的方式。
常见错误. 不要把仅对素数 p p p 定义、且 = 1 =1 = 1 就意味着 a a a 确实是模 p p p 平方数的勒让德 符号,与对奇合数 n n n 通过其素因子上勒让德符号相乘定义的雅可比 符号混淆:( a n ) Jacobi = 1 \left(\frac{a}{n}\right)_{\text{Jacobi}}=1 ( n a ) Jacobi = 1 并不 意味着 a a a 是模 n n n 的二次剩余——可能 a a a 对 n n n 的每个素因子都是非剩余,且对 n n n 的素因子为非剩余的情形出现偶数次,使雅可比符号为 + 1 +1 + 1 ,而 a a a 根本不是模 n n n 的平方数。 历史注记
高斯将二次互反律称为他的“theorema aureum”(黄金定理);他在19岁那年即1796年找到了第一个证明,并在一生中陆续发表了八种不同的证明,其中包括后来被艾森斯坦完善的格点计数论证。勒让德更早就猜到了这条定律,但未能给出完全一般的证明。
卡尔·弗里德里希·高斯
当 x 2 ≡ a ( m o d p ) x^2\equiv a\pmod p x 2 ≡ a ( mod p ) 无解时,( a p ) \left(\frac{a}{p}\right) ( p a ) 称为什么?
二次剩余 二次非剩余 原根 完全平方数
根据欧拉判别法,( a p ) \left(\frac{a}{p}\right) ( p a ) 模 p p p 同余于 a a a 的第几次幂?
a p − 1 a^{p-1} a p − 1 a p a^{p} a p a p − 1 2 a^{\frac{p-1}{2}} a 2 p − 1 a 2 a^{2} a 2 奇素数 p p p 模 4 4 4 属于哪一类时,− 1 -1 − 1 是二次剩余?
p ≡ 1 ( m o d 4 ) p \equiv 1 \pmod 4 p ≡ 1 ( mod 4 ) p ≡ 3 ( m o d 4 ) p \equiv 3 \pmod 4 p ≡ 3 ( mod 4 ) p ≡ 2 ( m o d 4 ) p \equiv 2 \pmod 4 p ≡ 2 ( mod 4 ) 从不 哪种密码学构件直接基于判定模合数的二次剩余性的困难性?
Goldwasser–Micali 密码系统 RSA 签名填充方案 AES 分组密码 SHA-256 哈希函数