基本概念和一次同余式

基本概念

同余式的定义:设mZ+,f(x)=anxn+an1xn1+....+a1x+a0,aiZm\in Z_{+},f(x)=a_nx_n+a_{n-1}x^{n-1}+....+a_{1}x+a_0,a_i\in Z,则f(x)0( mod m)f(x)\equiv0(~mod~m)叫做模m的同余式,若ai≢0 (mod m)a_i\not\equiv 0~(mod~m),则n叫做上述同余式的次数。

同余式解的定义:若a是使得f(a)0 (mod m)f(a)\equiv0~(mod~m)成立的一个整数,则称a所在的模m的剩余类kak_a为同余式的一个解(一类解),记为xa (mod m)x\equiv a~(mod~m)

一次同余式

一次同余式:形如axb (mod m)ax\equiv b~(mod~m),且满足a≢0 (mod m)a\not\equiv 0~(mod~m),这样形式的式子称为一次同余式。

定理1:一次同余式axb (mod m),a≢0 (mod m)  ax\equiv b~(mod~m),a\not\equiv 0~(mod~m)~~(a,m)=1(a,m)=1时有唯一解。即当未知数系数与模数互质的时候有唯一解。

定理2:设(a,m)=d(a,m)=d,则一次同余式axb (mod m),a≢0 (mod m)  ax\equiv b~(mod~m),a\not\equiv 0~(mod~m)~~,有解的充分必要条件是dbd\mid b,当该同余式有解时,它恰好有d个解

解一次同余方程

  • 熟悉完定理1和定理2的推导过程,其实可以了解到解一次同余方程其实就是解不定方程。下面俩个例题给出解同余方程的步骤。

求解同余方程步骤

  1. 先求解gcd(a,m)gcd(a,m),判断am是否互素。当gcd(a,m)gcd(a,m)互素可以得到该方程有唯一解,当不互素的时候令gcd(a,m)=dgcd(a,m)=d,判断d是否能整除b,能整除就说明有d个解,不能整除就说明该同余方程无解
  2. gcd(a,m)=1gcd(a,m)=1,直接去解二元不定方程as+mt=1as+mt=1,求出s,得到唯一解为xbs mod( m)x\equiv bs~mod(~m),就结束了。
  3. **当gcd(a,m)=dgcd(a,m)=d**并且同余方程有解时,我们就要将题目所给的同余方程中的a、b、m都除以d,得到一个唯一解的同余方程a1xb1 mod( m1)a_1x\equiv b_1~mod(~m_1),通过解二元不定方程as+mt=1as+mt=1,求出s,得到a1xb1 mod( m1)a_1x\equiv b_1~mod(~m_1)唯一解xb1s mod( m1)x\equiv b_1s~mod(~m_1)
  4. 求出原同余方程axb mod( m)ax\equiv b~mod(~m)的解,其解为xb1s+km1 mod( m),k=0,1,...,d1x\equiv b_1s + km_1~mod(~m),k=0,1,...,d-1

例题1:求一次同余方程2x3 (mod 5)2x\equiv 3~(mod~5)的所有解

例题2:求一次同余方程4x6(mod 10)4x\equiv 6(mod~10)的所有解

中国剩余定理(CRT)

中国剩余定理的个人认为算是基础部分最重要的一个定理,也是密码学中非常重要的一个定理。而二次互反律则是整个初等数论最重要的定理(听b站某位老师说的)。

  • 中国剩余定理是用来解同余方程组的,该同余方程组形式如下,特征有三个:
    • 未知数只有一个x
    • 每个同余式的模数不一样,但是一般要满足gcd(m1,m2,...,mk)=1才有解,否则需要使用拓展中国剩余定理
    • 每个同余式中b的值可以一样,也可以不一样

{xb1( mod m1)xb2( mod m2)............xbk( mod mk)                             (1)\begin{cases} x \equiv b_1(~mod~m_1)\\ x \equiv b_2(~mod~m_2)\\ ............\\ x \equiv b_k(~mod~m_k) \end{cases}~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(1)

  • 先来说解法再说其相关定理,这里默认gcd(m1,m2,...,mk)=1gcd(m_1,m_2,...,m_k)=1,按照从左到右的顺序计算出如下数据。
除数 余数 除数最小公倍数 衍数 乘率 各总
m1m_1 b1b_1 m1m2....mk=mm_1*m_2....m_k=m M1=mm1M_1=\frac{m}{m_1} y1=M11 (mod m1)y_1=M_1^{-1}~(mod~m_1) b1M1y1b_1*M_1*y_1
m2m_2 b2b_2 m1m2....mk=mm_1*m_2....m_k=m M2=mm2M_2=\frac{m}{m_2} y2=M21 (mod m2)y_2=M_2^{-1}~(mod~m_2) b2M2y2b_2*M_2*y_2
m1m2....mk=mm_1*m_2....m_k=m
mkm_k bkb_k m1m2....mk=mm_1*m_2....m_k=m Mk=mmkM_k=\frac{m}{m_k} yk=Mk1 (mod mk)y_k=M_k^{-1}~(mod~m_k) bkMky3b_k*M_k*y_3
  • 最后得到最终结果:

xb1M1y1+b2M2y2+....+bkMkyk mod( m)x\equiv b_1*M_1*y_1 + b_2*M_2*y_2+....+b_k*M_k*y_k~mod(~m)

定理1(中国剩余定理):

m1,m2,...,mkm_1,m_2,...,m_k是k个两两互素的正整数,m=m1m2...mk,m=miMi(i=1,2,...,k)m=m_1m_2...m_k,m=m_iM_i(i=1,2,...,k),则同余式组(1)的解为:

xM1M1b1+M2M2b2+....+MkMkbk      (2)x\equiv M_1'M_1b_1+M_2'M_2b_2+....+M_k'M_kb_k~~~~~~(2)

其中MiMi1( mod mi),i=1,2,...,kM_i'M_i\equiv1(~mod~m_i),i=1,2,...,k

定理2

b1,b2,...,bkb_1,b_2,...,b_k分别过模m1,m2,...,mkm_1,m_2,...,m_k的完全剩余系,则(2)过模m=m1m2...mkm=m_1m_2...m_k的完全剩余系。

高次同余式的解数及解法

相关定理

定理1

m1,m2,...,mkm_1,m_2,...,m_k是k个两两互素的正整数,m=m1m2...mkm=m_1*m_2*...*m_k,则同余式f(x)0( mod m)f(x)\equiv0(~mod~m)与下面同余式组等价(即任意适合f(x)0( mod m)f(x)\equiv0(~mod~m)都适合如下式子,反之也成立)。

f(x)0( mod mi),i=1,2,...,kf(x)\equiv0(~mod~m_i),i=1,2,...,k

并且若用TiT_i表示f(x)0(mod mi),i=1,2,...,kf(x)\equiv0(mod~m_i),i=1,2,...,k对模mim_i的解数,T表示f(x)0( mod m)f(x)\equiv0(~mod~m)的解数,则解数之间有如下关系

T=T1T2...TKT=T_1T_2...T_K

定理2——命题1

α\alpha是大于等于2的正整数,则适合f(x)0 (mod pα)f(x)\equiv0~(mod~p^{\alpha}),p为素数。适合上式的每个整数都适合同余式f(x)0 (pα1)f(x)\equiv 0~(p^{\alpha-1})

进而适合f(x)0 (mod pα)f(x)\equiv0~(mod~p^{\alpha})的每个正数都适合同余式f(x)0 (mod p)f(x)\equiv0~(mod~p)

因此解f(x)0 (mod pα)f(x)\equiv0~(mod~p^{\alpha}),可以转换为求f(x)0 (mod pα)f(x)\equiv0~(mod~p^{\alpha})的解。

定理2

xx1 (mod p)x\equiv x_1~(mod~p)x=x1+pt1,tZx=x_1+pt_1,t\in Zf(x)0 (mod p)f(x)\equiv0~(mod~p)的一解,并且pf(x1)p\nmid f'(x_1)(f(x)f'(x)f(x)f(x)的导函数),

x=x1+pt1x=x_1+pt_1刚好给出f(x)0 (mod pα)f(x)\equiv0~(mod~p^{\alpha})的一解,对于模pαp^{\alpha}来说:

x=xα+pαtα,tαZx=x_{\alpha}+p^{\alpha}t_{\alpha},t_{\alpha}\in Z

xxα( mod pα)x\equiv x_{\alpha}(~mod~p^{\alpha}),其中xαx1 mod( p)x_{\alpha}\equiv x_1~mod(~p)

解高次同余式

情形一:

f(x)0( mod m)m=m1m2...mkf(x)\equiv0(~mod~m)\\ m=m_1m_2...m_k

解法:

  1. 利用定理一将同余式转化为f(x)0( mod mi),i=1,2,...,kf(x)\equiv0(~mod~m_i),i=1,2,...,k
  2. 分别解出这些转化后的同余方程的解
  3. 列出同余方程组,使用CRT即可求得原方程的解

{xb1 (mod m1)xb2 (mod m2).....xbk (mod mk)\begin{cases} x\equiv b_1 ~(mod~m_1)\\ x\equiv b_2 ~(mod~m_2)\\ .....\\ x\equiv b_k ~(mod~m_k) \end{cases}

情形二:

f(x)0( mod pα)f(x)\equiv0(~mod~p^{\alpha})

解法:

  1. 转化为解同余方程f(x)0 (mod p)f(x)\equiv0~(mod~p)的解

  2. 求出f(x)0 (mod p)f(x)\equiv0~(mod~p)的解后可以得到x=x1+pt1,tZx=x_1+pt_1,t\in Z,将x1x_1带入f(x)+f(x)pt10(mod p2)f(x)+f'(x)*pt_1\equiv0(mod~p^{2}),并求出t1x2 (mod p2)t_1\equiv x_2~(mod~p^{2}),这样就有t1=x2+p2t2t_1=x_2+p^{2}*t_2

  3. t1=x2+p2t2t_1=x_2+p^{2}*t_2带回去x=x1+pt1x=x_1+pt_1,可以得到x=x1+p(x2+p3t2)=(x1+px2)+p3t2x=x_1+p*(x_2+p^{3}*t_2)=(x_1+p*x_2)+p^{3}t_2从而得到xx1+px2 mod( p2)x\equiv x_1+p*x_2~mod(~p^{2})

  4. 之后将x=(x1+px2)+p3t2x=(x_1+p*x_2)+p^{3}t_2,带入到f(x)+f(x)p2t20(mod p3)f(x)+f'(x)*p^{2}t_2\equiv0(mod~p^{3})并且重复2、3步骤,直到带回到f(x)0( mod pα)f(x)\equiv0(~mod~p^{\alpha})

情形三:

f(x)0( mod m)m=p1α1p2α2...pαnf(x)\equiv0(~mod~m)\\ m=p_1^{\alpha_1}p_2^{\alpha_2}...p^{\alpha_{n}}

  1. 先使用定理一,将同余方程转换为f(x)0( mod piαi)f(x)\equiv0(~mod~p_i^{\alpha_i}),即转换为情形二的情况

  2. 根据情形二的情况去解出f(x)0( mod piαi)f(x)\equiv0(~mod~p_i^{\alpha_i})的解,每组方程都需要求出解。

  3. 求出每个方程的所有解后,可以构造出如下同余方程组,使用CRT,即可求出情形三的所有解:

{xb1 (mod p1α1)xb2 (mod p2α2).....xbk (mod pkαk)\begin{cases} x\equiv b_1 ~(mod~p^{\alpha_1}_1)\\ x\equiv b_2 ~(mod~p^{\alpha_2}_2)\\ .....\\ x\equiv b_k ~(mod~p^{\alpha_k}_k) \end{cases}

素数模的同余式

主要是利用如下定理化简素数模的多项式同余式

定理1

设p是素数,则同余式f(x)=anxn+an1xn1+...+a00( mod p),an≢0 mod( p)f(x)=a_nx^{n}+a_{n-1}x^{n-1}+...+a_0 \equiv 0(~mod~p),a_n\not\equiv 0~mod(~p),与一个次数不超过p-1的素数模p同余式等价。

定理2

kn,xαi( mod p)(i=1,2,...,k)k≤n,x\equiv \alpha_i(~mod~p)(i=1,2,...,k)是同余式f(x)=anxn+an1xn1+...+a00( mod p),an≢0 mod( p)f(x)=a_nx^{n}+a_{n-1}x^{n-1}+...+a_0 \equiv 0(~mod~p),a_n\not\equiv 0~mod(~p)的k个不同解,则对任何整数x来说,有如下式子:

f(x)(xα1)(xα2)...(xαk)fk(x)( mod p)f(x)\equiv(x-\alpha_1)(x-\alpha_2)...(x-\alpha_k)f_k(x)(~mod~p)

其中fk(x)f_k(x)是首项系数ana_nnkn-k次多项式.

定理3——(1)

对于任何整数x来说,都有如下式子:

xp1(x1)(x2)...(x(p1)) (mod p)x^{p-1}\equiv(x-1)(x-2)...(x-(p-1))~(mod~p)

定理3——(2):又称威尔逊定理

(p1)!+10( mod p)(p-1)!+1\equiv0(~mod~p)

定理4

同余式f(x)=anxn+an1xn1+...+a00( mod p),an≢0 mod( p)f(x)=a_nx^{n}+a_{n-1}x^{n-1}+...+a_0 \equiv 0(~mod~p),a_n\not\equiv 0~mod(~p)的解数不超过它的次数。

定理5——命题1:

同余式f(x)=anxn+an1xn1+...+a00( mod p),an≢0 mod( p)f(x)=a_nx^{n}+a_{n-1}x^{n-1}+...+a_0 \equiv 0(~mod~p),a_n\not\equiv 0~mod(~p)总和某个次数的首一同余式等价。

定理5

npn≤p,则同余式f(x)=xn+an1xn1+...+a00( mod p)f(x)=x^{n}+a_{n-1}x^{n-1}+...+a_0 \equiv 0(~mod~p),有n个解的充分必要条件是(xpx)/f(x)(x^{p}-x)/f(x)所得余式的一切系数都是p的倍数。