同余方程

定义

同余方程

对正整数 mm 和一元整系数多项式 f(x)=∑i=0naixif(x)=\sum_{i=0}^n a_ix^i,其中未知数 x∈Zmx\in\mathbf{Z}_m,称形如

f(x)≡0(modm)(1)f(x)\equiv 0\pmod m\tag{1}

的方程为关于未知数 xx 的模 mm 的一元 同余方程(Congruence Equation).

若 an≢0(modm)a_n\not\equiv 0\pmod m,则称上式为 nn 次同余方程.

类似可定义同余方程组.

关于一次同余方程与方程组的相关内容请参见 线性同余方程 与 中国剩余定理.

本文首先研究同余方程的可解性和解集结构,之后会简要介绍高次同余方程的解法.

由 中国剩余定理 可知,求解关于模合数 mm 的同余方程可转化为求解模素数幂次的情况.故以下只介绍素数幂模同余方程和素数模同余方程的相关理论.

素数幂模同余方程

以下假设模数 m=pe (p∈P, e∈Z>1)m=p^e~(p\in\mathbf{P},~e\in\mathbf{Z}_{>1}).

注意到,若 x0x_0 是方程

f(x)≡0(modpe)f(x)\equiv 0\pmod{p^e}

的解,则 x0x_0 是方程

f(x)≡0(modpe−1)f(x)\equiv 0\pmod{p^{e-1}}

的解.这启发我们尝试通过较低的模幂次的解去构造较高的模幂次的解.我们有如下定理:

定理 1(Hensel 引理)

对素数 pp 和整数 e>1e>1,取整系数多项式 f(x)=∑i=0naixi (pe∤an)f(x)=\sum_{i=0}^na_ix^i~(p^e\nmid a_n),令 f′(x)=∑i=1niaixi−1f'(x)=\sum_{i=1}^nia_ix^{i-1} 为其导数.令 x0x_0 为方程

f(x)≡0(modpe−1)(2)f(x)\equiv 0\pmod{p^{e-1}}\tag{2}

的解,则:

  1. 若 f′(x0)≢0(modp)f'(x_0)\not\equiv 0\pmod p,则存在整数 tt 使得

    x=x0+pe−1t(3)x=x_0+p^{e-1}t \tag{3}

    是方程

    f(x)≡0(modpe)(4)f(x)\equiv 0\pmod{p^e} \tag{4}

    的解.

  2. 若 f′(x0)≡0(modp)f'(x_0)\equiv 0\pmod p 且 f(x0)≡0(modpe)f(x_0)\equiv 0\pmod{p^e},则对 t=0,1,…,p−1t=0,1,\dots,p-1,由式 (3)(3) 确定的 xx 均为方程 (4)(4) 的解.

  3. 若 f′(x0)≡0(modp)f'(x_0)\equiv 0\pmod p 且 f(x0)≢0(modpe)f(x_0)\not\equiv 0\pmod{p^e},则不能由式 (3)(3) 构造方程 (4)(4) 的解.

证明

我们假设式 (3)(3) 是方程 (4)(4) 的解,即

f(x0+pe−1t)≡0(modpe)f(x_0+p^{e-1}t)\equiv 0\pmod{p^e}

整理后可得

f(x0)+pe−1tf′(x0)≡0(modpe)f(x_0)+p^{e-1}tf'(x_0)\equiv 0\pmod{p^e}

于是

tf′(x0)≡−f(x0)pe−1(modp)(5)tf'(x_0)\equiv -\frac{f(x_0)}{p^{e-1}}\pmod p\tag{5}
  1. 若 f′(x0)≢0(modp)f'(x_0)\not\equiv 0\pmod p,则关于 tt 的方程 (5)(5) 有唯一解 t0t_0,代入式 (3)(3) 可验证其为方程 (4)(4) 的解.
  2. 若 f′(x0)≡0(modp)f'(x_0)\equiv 0\pmod p 且 f(x0)≡0(modpe)f(x_0)\equiv 0\pmod{p^e},则任意 tt 均能使方程 (5)(5) 成立,代入式 (3)(3) 可验证其均为方程 (4)(4) 的解.
  3. 若 f′(x0)≡0(modp)f'(x_0)\equiv 0\pmod p 且 f(x0)≢0(modpe)f(x_0)\not\equiv 0\pmod{p^e},则方程 (5)(5) 无解,从而不能由式 (3)(3) 构造方程 (4)(4) 的解.

进而我们有推论:

推论 1

对 定理 1 的 pp,ee,f(x)f(x),x0x_0,

  1. 若 ss 是方程 f(x)≡0(modp)f(x)\equiv 0\pmod p 的解,且 f′(s)≢0(modp)f'(s)\not\equiv 0\pmod p,则存在 xs∈Zpex_s\in\mathbf{Z}_{p^e},xs≡s(modp)x_s\equiv s\pmod p 使得 xsx_s 是方程 (4)(4) 的解.
  2. 若方程 f(x)≡0(modp)f(x)\equiv 0\pmod p 与 f′(x)≡0(modp)f'(x)\equiv 0\pmod p 无公共解,则方程 (4)(4) 和方程 f(x)≡0(modp)f(x)\equiv 0\pmod p 的解数相同.

从而我们可以将素数幂模同余方程化归到素数模同余方程的情况.

素数模同余方程

以下令 p∈Pp\in\mathbf{P},整系数多项式 f(x)=∑i=0naixif(x)=\sum_{i=0}^na_ix^i,其中 p∤anp\nmid a_n,x∈Zpx\in\mathbf{Z}_p.

定理 2

若方程

f(x)≡0(modp)(6)f(x)\equiv 0\pmod p\tag{6}

有 kk 个不同的解 x1,x2,…,xk (k≤n)x_1,x_2,\dots,x_k~(k\leq n),则

f(x)≡g(x)∏i=1k(x−xi)(modp),f(x)\equiv g(x)\prod_{i=1}^k(x-x_i)\pmod p,

其中 deg⁡g=n−k\deg g=n-k 且 [xn−k]g(x)=an[x^{n-k}]g(x)=a_n.

证明

对 kk 应用数学归纳法.

  • 当 k=1k=1 时,做多项式带余除法,有 f(x)=(x−x1)g(x)+rf(x)=(x-x_1)g(x)+r,其中 r∈Zr\in\mathbf{Z}.

    由 f(x1)≡0(modp)f(x_1)\equiv 0\pmod p 可知 r≡0(modp)r\equiv 0\pmod p,从而 f(x)≡(x−x1)g(x)(modp)f(x)\equiv(x-x_1)g(x)\pmod p.

  • 假设命题对 k−1k-1(k>1k>1) 时的情况成立,现在设 f(x)f(x) 有 kk 个不同的解 x1,x2,…,xkx_1,x_2,\dots,x_k,则 f(x)≡(x−x1)h(x)(modp)f(x)\equiv(x-x_1)h(x)\pmod p,进而有

    (∀i=2,3,…,k),  0≡f(xi)≡(xi−x1)h(xi)(modp)(\forall i=2,3,\dots,k),~~0\equiv f(x_i)\equiv (x_i-x_1)h(x_i)\pmod p

    从而 h(x)h(x) 有 k−1k-1 个不同的解 x2,x3,…,xkx_2,x_3,\dots,x_k,由归纳假设有

    h(x)≡g(x)∏i=2k(x−xi)(modp)h(x)\equiv g(x)\prod_{i=2}^k(x-x_i)\pmod p

    其中 deg⁡g=n−k\deg g=n-k 且 [xn−k]g(x)=an[x^{n-k}]g(x)=a_n.

    因此命题得证.

推论 2

对素数 pp,

  • (∀x∈Z),  xp−1−1≡∏i=1p−1(x−i)(modp)(\forall x\in\mathbf{Z}),~~x^{p-1}-1 \equiv \prod_{i=1}^{p-1}(x-i)\pmod p.
  • (Wilson 定理)(p−1)!≡−1(modp)(p-1)! \equiv -1 \pmod p.
定理 3(Lagrange 定理)

方程 (6)(6) 至多有 nn 个不同解.

证明

假设 f(x)f(x) 有 n+1n+1 个不同解 x1,x2,…,xn+1x_1,x_2,\dots,x_{n+1},则由 定理 2,对 x1,x2,…,xnx_1,x_2,\dots,x_n 有

f(x)≡an∏i=1n(x−xi)(modp)f(x)\equiv a_n\prod_{i=1}^n(x-x_i)\pmod p

令 x=xn+1x=x_{n+1},则

0≡f(xn+1)≡an∏i=1n(xn+1−xi)(modp)0\equiv f(x_{n+1})\equiv a_n\prod_{i=1}^n(x_{n+1}-x_i)\pmod p

而右侧显然不是 pp 的倍数,因此假设矛盾.

推论 3

若同余方程 ∑i=0nbixi≡0(modp)\sum_{i=0}^nb_ix^i\equiv 0\pmod p 的解数大于 nn,则

(∀i=0,1,…,n),  p∣bi.(\forall i=0,1,\dots,n),~~p\mid b_i.
定理 4

方程 (6)(6) 若解的个数不为 pp,则必存在满足 $\deg r [!note]+ 证明

不妨设 n≥pn\geq p,对 f(x)f(x) 做多项式带余除法

f(x)=g(x)(xp−x)+r(x)f(x)=g(x)\left(x^p-x\right)+r(x)

其中 deg⁡r由[Fermat小定理](/knowledge−pages/math−number−theory−fermat)知对任意整数\deg r 由 [Fermat 小定理](/knowledge-pages/math-number-theory-fermat) 知对任意整数 x有有x^p\equiv x\pmod p$,从而

  • 若 r(x)≡0(modp)r(x)\equiv 0\pmod p,则由 推论 2 可知 f(x)f(x) 有 pp 个不同的解.
  • 若 r(x)≢0(modp)r(x)\not\equiv 0\pmod p,则由 f(x)≡r(x)(modp)f(x)\equiv r(x)\pmod p 可知 f(x)f(x) 和 r(x)r(x) 的解集相同.

我们可以通过这个定理对同余方程降次.

定理 5

设 n≤pn\leq p,则方程

xn+∑i=0n−1aixi≡0(modp)(7)x^n+\sum_{i=0}^{n-1}a_ix^i\equiv 0\pmod p\tag{7}

有 nn 个解当且仅当存在整系数多项式 q(x)q(x),r(x) (deg⁡r<n)r(x)~(\deg r < n) 使得

xp−x=f(x)q(x)+pr(x).(8)x^p-x=f(x)q(x)+pr(x). \tag{8}
证明
  • 必要性:由多项式除法知存在整系数多项式 q(x)q(x),r1(x) (deg⁡r1<n)r_1(x)~(\deg r_1 < n) 使得

    xp−x=f(x)q(x)+r1(x).x^p-x=f(x)q(x)+r_1(x).

    若方程 (7)(7) 有 nn 个解,则 r1≡0(modp)r_1\equiv 0\pmod p 也有 nn 个相同的解,进而由 推论 3 可知存在整系数多项式 r(x)r(x) 满足 r1(x)=pr(x)r_1(x)=pr(x),从而命题得证.

  • 充分性:若式 (8)(8) 成立,则由 Fermat 小定理 可知,对任意整数 xx,

    0≡xp−x≡f(x)q(x)(modp).0\equiv x^p-x\equiv f(x)q(x)\pmod p.

    即方程 f(x)q(x)≡0(modp)f(x)q(x)\equiv 0\pmod p 有 pp 个解.

    设方程 (7)(7) 的解数为 ss,则由 Lagrange 定理 可知 s≤ns\leq n.

    又由于 deg⁡q=p−n\deg q=p-n,则由 Lagrange 定理 可知 q(x)≡0(modp)q(x)\equiv 0\pmod p 的解数不超过 p−np-n,而方程 f(x)q(x)≡0(modp)f(x)q(x)\equiv 0\pmod p 的解集是 f(x)≡0(modp)f(x)\equiv 0\pmod p 解集和 q(x)≡0(modp)q(x)\equiv 0\pmod p 解集的并集,故 s+(p−n)≥ps+(p-n)\geq p,有 s≥ns\geq n.

    因此 s=ns=n.

对于非首一多项式,由于 Zp\mathbf{Z}_p 是域,故可以将其化为首一多项式,从而适用该定理.

定理 6

设 n∣p−1n\mid p-1,p∤ap\nmid a,则方程

xn≡a(modp)(9)x^n\equiv a\pmod p\tag{9}

有解当且仅当

ap−1n≡1(modp).a^{\frac{p-1}{n}}\equiv 1\pmod p.

而且,若 (9)(9) 有解,则解数为 nn.

Note

方程 (9)(9) 解集的具体结构可参见 k 次剩余.

证明
  • 必要性:若方程 (9)(9) 有解 x0x_0,则

    ap−1n≡(x0n)p−1n≡1(modp)a^{\frac{p-1}{n}}\equiv {\left(x_0^n\right)}^{\frac{p-1}{n}}\equiv 1\pmod p
  • 充分性:若 ap−1n≡1(modp)a^{\frac{p-1}{n}}\equiv 1\pmod p,则

    xp−x=x(xp−1−1)=x((xn)p−1n−ap−1n+ap−1n−1)=(xn−a)P(x)+x(ap−1n−1)\begin{aligned} x^p-x&=x\left(x^{p-1}-1\right)\\ &=x\left(\left(x^n\right)^{\frac{p-1}{n}}-a^{\frac{p-1}{n}}+a^{\frac{p-1}{n}}-1\right)\\ &=\left(x^n-a\right)P(x)+x\left(a^{\frac{p-1}{n}}-1\right)\\ \end{aligned}

    其中 P(x)P(x) 是某个整系数多项式,因此由 定理 5 可知方程 (9)(9) 有 nn 个解.

高次同余方程(组)的求解方法

首先我们可以借助 中国剩余定理 将求解 同余方程组 转为求解 同余方程,以及将求解模 合数 mm 的同余方程转化为求解模 素数幂次 的同余方程.之后我们借助 定理 1 将求解模 素数幂次 的同余方程转化为求解模 素数 的同余方程.

结合模素数同余方程的若干定理,我们只需考虑方程

xn+∑i=0n−1aixi≡0(modp)x^n+\sum_{i=0}^{n-1}a_ix^i\equiv 0\pmod p

的求法,其中 pp 是素数,n<pn<p.

我们可以通过将 xx 代换为 x−an−1nx-\dfrac{a_{n-1}}{n} 来消去 xn−1x^{n-1} 项,从而我们只需考虑方程

xn+∑i=0n−2aixi≡0(modp)(10)x^n+\sum_{i=0}^{n-2}a_ix^i\equiv 0\pmod p\tag{10}

的求法,其中 pp 是素数,n<pn<p.

  • 若 n=1n=1,则求法参见 线性同余方程.

  • 若 n=2n=2,则求法参见 二次剩余.

  • 若方程 (10)(10) 可化为

    xn≡a(modp),x^n\equiv a\pmod p,

    则求法参见 k 次剩余.

参考资料

  1. Congruence Equation -- from Wolfram MathWorld
  2. Lagrange's theorem (number theory) - Wikipedia
  3. 潘承洞,潘承彪.初等数论.
  4. 冯克勤.初等数论及其应用.
  5. 闵嗣鹤,严士健.初等数论.

本文内容来自 OI Wiki《同余方程》;作者与修订历史见原文。AIOJ 对排版、站内链接和资源路径做了格式适配。

内容依 CC BY-SA 4.0 与 SATA 许可共享;附加条款以原文声明为准。感谢 OI Wiki 作者与贡献者;请按 SATA 约定前往 OI Wiki 仓库 Star 支持。