初等数论入门 Lesson 6 同余方程与中国剩余定理

剩余系

  • 完全剩余系:一组数 a1,a2,…,ama_1,a_2,\dots,a_m 称为模 mm 的一个完全剩余系,如果它们分别取自模 mm 的 mm 个不同的剩余类。

    有以下两种典型的剩余系:

    • 最小非负剩余系:{0,1,2,…,m−1}\{0, 1, 2, \dots, m - 1\}
    • 绝对最小剩余系:{−⌊m2⌋,…,−1,0,1,…,⌈m2⌉}\left\{ -\left\lfloor \frac{m}{2} \right\rfloor, \dots, -1, 0, 1, \dots, \left\lceil \frac{m}{2} \right\rceil \right\}

    有一个性质:

    • 如果 {a1,…,am}\{a_1,\dots,a_m\} 是完全剩余系,且gcd⁡(a,m)=1\gcd(a,m)=1,则 {aa1+b,aa2+b,…,aam+b}\{aa_1 + b, aa_2 + b, \dots, aa_m + b\} 也是一个模 mm 的完全剩余系。
  • 缩系(简化剩余系)

    在完全剩余系中,如果我们只挑与 mm 互素的代表元,就得到了缩系。

    性质:乘法封闭性

    • 如果 {r1,r2,…,rφ(m)}\{r_1,r_2,\dots,r_{\varphi(m)}\} 是模 mm 的缩系,且 a∈{ri}a\in \{r_i\},则 {ar1,ar2,…,arφ(m)}\{ar_1,ar_2,\dots,ar_{\varphi(m)}\} 仍然为模 mm 的缩系。

欧拉定理

设 m≥2m\ge 2,且 (a,m)=1(a,m)=1,则:

aφ(m)≡1(modm)\boxed{a^{\varphi(m)}\equiv 1\pmod m}

推导过程:

  1. 取 mm 的缩系 {r1,r2,…,rφ(m)}\{r_1,r_2,\dots,r_{\varphi(m)}\}

  2. 又 (a,m)=1(a,m)=1,所以 {ar1,ar2,…,arφ(m)}\{ar_1,ar_2,\dots,ar_{\varphi(m)}\} 仍为模 mm 的缩系

  3. 因此:

    r1r2⋯rφ(m)≡(ar1)(ar2)⋯(arφ(m))(modm)r_1r_2\cdots r_{\varphi(m)}\equiv (ar_1)(ar_2)\cdots(ar_{\varphi(m)})\pmod m

  4. 变化得:

    r1r2⋯rφ(m)≡aφ(m)(r1r2⋯rφ(m))(modm)r_1r_2\cdots r_{\varphi(m)}\equiv a^{\varphi(m)}(r_1r_2\cdots r_{\varphi(m)})\pmod m

  5. 又 ∀ri:(ri,m)=1\forall r_i:(r_i,m)=1, 故:

    aφ(m)≡1(modm)a^{\varphi(m)}\equiv 1\pmod m

    成立

费马小定理

当 pp 为素数,且 p∤ap\nmid a,则:

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

费马小定理本质上是欧拉定理 mm 为素数的一个特例。

应用:降幂:化简大指数

一元线性同余方程

求解:

ax≡b(modm)ax\equiv b\pmod m

等价于:

ax−b=myax+m′y=bax-b=my \\ ax + m'y =b

于是 解同余方程变成了解不定方程

我们规范 ax≡b(modm)ax\equiv b\pmod m 求解步骤如下,令 d=gcd⁡(a,m)d=\gcd(a,m)

  1. 判断:若 d∤bd\nmid b,则无解

  2. 约简(我们在此处约简,而不再解不定方程处约简):

    adx≡bd(modmd)\dfrac a d x\equiv \dfrac b d\pmod {\dfrac m d}

  3. 求逆元,我们要求一个 tt 满足:

    ad⋅t≡1(modmd)\dfrac a d \cdot t \equiv 1\pmod {\dfrac m d}

    这一步我们就使用拓展欧几里得算法求出:

    ad⋅t+md⋅s=1\dfrac a d \cdot t+\dfrac m d \cdot s =1

  4. 求特解与通解:

    • 特解:x0≡t⋅bd(modmd)x_0\equiv t \cdot \dfrac b d\pmod {\dfrac m d}
    • 通解:x≡x0+k⋅md(modm)x\equiv x_0+k\cdot \dfrac m d\pmod m,其中 k=0,1,…,d−1k=0,1,\dots,d-1

    因此模 mm 意义下共有 dd 个不同的解。

中国剩余定理

Chinese Remainder Theorem CRT

模数互素情况

已知 mim_i 两两互质,求解:

{x≡a1(modm1)x≡a2(modm2)⋯x≡ak(modmk)\begin{cases} x \equiv a_1 \pmod{m_1} \\ x \equiv a_2 \pmod{m_2} \\ \cdots \\ x \equiv a_k \pmod{m_k} \end{cases}

  1. 计算 M=m1m2⋯mkM=m_1m_2\cdots m_k,这是解的周期

  2. 对于每个 ii,令 Mi=MmiM_i=\dfrac M{m_i},必有 gcd⁡(Mi,mi)=1\gcd(M_i,m_i)=1

  3. 求解 tit_i:

    Miti≡1(modmi)M_it_i\equiv 1\pmod {m_i}

  4. 构造解:

    x=a1M1t1+a2M2t2+⋯+akMktkx=a_1M_1t_1+a_2M_2t_2+\cdots+a_kM_kt_k

CRT断言:该方程组在模 MM 意义下有唯一解

拓展中国剩余定理

模数不互质情况

我们先来解决两个方程的情况:

{x≡a1(modm1)x≡a2(modm2)\begin{cases} x\equiv a_1 \pmod{m_1}\\ x\equiv a_2\pmod {m_2} \end{cases}

令 d=gcd⁡(m1,m2)d=\gcd(m_1,m_2)

可解性定理:该方程组有解,当且仅当:

a1≡a2(modd)a_1\equiv a_2\pmod d

证明:

因为 d∣m1,  d∣m2d\mid m_1,\;d\mid m_2,故:

{x−a1≡0(modd)x−a2≡0(modd)\begin{cases} x-a_1\equiv 0 \pmod{d}\\ x-a_2\equiv 0\pmod {d} \end{cases}

因此:

a1−a2≡0(modd)a_1-a_2\equiv 0\pmod d

步骤:

  1. 由方程1:x=a1+m1kx=a_1+m_1k

  2. 代入方差2:a1+m1k≡a@(modm2)a_1+m_1k\equiv a_@\pmod {m_2}

    即:

    m1k≡a2−a1(modm2)m_1k\equiv a_2-a_1\pmod {m_2}

  3. 解关于 kk 的同余方差。由于 gcd⁡(m1,m2)=d∣(a2−a1)\gcd(m_1,m_2)=d|(a_2-a_1),故方程有解

  4. 设解得 k≡k0(modm2d)k\equiv k_0\pmod {\dfrac{m_2}d},代回得:

    x=a1+m1k0+m1m2d⋅tx=a_1+m_1k_0+\dfrac{m_1m_2}d\cdot t

  5. 合并后的方程:

    x≡x0(modlcm(m1,m2))x\equiv x_0\pmod {\text{lcm}(m_1,m_2)}

    其中 x0=a1+m1k0x_0=a_1+m_1k_0 (解)

多方程推广

策略:两两合并

若某一步 ai≢aj(modgcd⁡(mi,mj))a_i \not\equiv a_j \pmod{\gcd(m_i, m_j)},则无解