二次剩余学习笔记

Poemsunny

已于 2025-02-07 14:25:53 修改

阅读量829

点赞数 17

文章标签：学习

于 2024-07-28 20:43:39 首次发布

本文链接：https://blog.csdn.net/weixin_73071245/article/details/140754884

版权

二次剩余

二次剩余，是对下面这个方程求解
$x^2 \equiv n \pmod p\\$
（其中 $p$ 为奇素数）

ps.以下内容都是建立在 $n\not=0$ 的基础上的，因为 $n = 0$ 时显然，且答案是在 $\bmod\ p$ 意义下的。

我们使用欧拉准则来判断这个方程是否有解。
由费马小定理，得
$n^{p-1} \equiv 1 \pmod p\\ \because p为奇素数\\ \therefore n^{2(\frac{p-1}{2})}\equiv 1\pmod p\\ \therefore n^{\frac{p-1}{2}}\equiv 1\pmod p或n^{\frac{p-1}{2}}\equiv -1\pmod p$
若这个方程有解，则有
$x^{p-1}\equiv 1\pmod p或x^{p-1}\equiv -1\pmod p$
后者显然不成立，所以 $n^{\frac{p-1}{2}}\equiv -1\pmod p$ 时方程无解。

下面，我们证 $n^{\frac{p-1}{2}}\equiv 1\pmod p$ 时方程有解。

我们将 $n$ 用 $g^k$ 表示，其中 $g$ 为 $p$ 的原根，则有
$g^{k(\frac{p-1}{2})}\equiv 1\pmod p\\ \because p-1|k(\frac{p-1}{2})\\ \therefore 2|k，即 k 为偶数\\ \therefore 原方程的一个解即为g^{\frac{k}{2}}$
这样我们就证明了 $n^{\frac{p-1}{2}}\equiv 1\pmod p$ 是原方程有解的充要条件。

那么，这个方程究竟有几个解呢？我们先假设这个方程有多个解，那么对于任意两个解 $x_1,x_2$ ，有
$x_1^{2} \equiv x_{2}^{2} \pmod p\\ \therefore x_1^2-x_2^2\equiv 0\pmod p\\ \therefore (x_1+x_2)(x_1-x_2)\equiv 0\pmod p\\ \because x_1,x_2\leq p\\ \therefore x_1-x_2\not\equiv 0\\ \therefore x_1+x_2\equiv 0\\$
由此可得，方程有且仅有两个解，且它们互为相反数。

因为是在 $\bmod \ p$ 意义下，所以 $x$ 的相反数为 $p - x$ 。

那么对于 $1$ 至 $p - 1$ 的 $\frac{p-1}{2}$ 对相反数，每一对都对应一个二次剩余，且互不相等，也就是说， $1$ 至 $p - 1$ 中，有 $\frac{p-1}{2}$ 个数是二次剩余，而有 $\frac{p-1}{2}$ 个数不是二次剩余。

我们用 Cipolla 算法对这个方程进行求解。

先在 $[1, p - 1]$ 中随机一个 $a$ ，使 $a^2-n$ 是非二次剩余，因为 $[1, p - 1]$ 中非二次剩余的数有 $\frac{p-1}{2}$ 个，所以我们的期望随机次数是 $2$ 次。

接下来，我们定义
$i^2 \equiv a^2-n \pmod p$
（我知道这个方程无解，但是这里的 $i$ 其实是类似复数的一个定义，我们定义这样的一个数，使得它的平方为 $a^2-n$ ）

方程的解即为 $\pm(a+i)^{\frac{p+1}{2}}$ 。

我们考虑证明。
将 $x$ 代入，得
$a+i)^{p+1}$
我们可以将其拆成
$a+i)^p(a+i)$
由二项式定理，得
$(a+i)^p=\sum_{j=0}^{p}{p \choose j}a^ji^{p-j}$
在 $\bmod \ p$ 意义下，因为 $j\in[1,p-1]$ 时组合数中都含有 $p$ 这一项，所以 $\bmod \ p$ 为 $0$ ，因此
$(a+i)^p\equiv a^p+i^p\equiv a+i^p \pmod p$
我们将 $i^p$ 进行化简，得
$i^p\equiv i(i^2)^{\frac{p-1}{2}}\equiv i(a^2-n)^{\frac{p-1}{2}}\equiv -i\pmod p$
原式即为
$(a-i)(a+i)\equiv a^2-i^2\equiv n\pmod p$