数论与代数

前几天在b站上看到一本书 CINTA(a Concrete Introduction to Number Theory and Algebra),感觉内容着实不错,几乎把同余这一块的内容讲完了,本文将会参考这本书,可以当作这本书的一个阅读笔记。

Github

其实这一块的内容并不常考,需要用到的时候通常没有很强的技巧,所以记好模板是很重要的一步。本文并不局限于这一本书,将有更多OI结合的应用。

写于2023年1月30日,深圳罗湖家中。


在原书的内容上添加了一部分东西

2.4 类欧几里德算法

2.5 万能欧几里德算法

4.4 扩展欧拉定理

10.2 扩展中国剩余定理

11.5 Cipolla 算法

15.5 Index Calculus 算法

17 连分数

18 二次域


1 整数和二进制#

只用 01 来表示所有数值的系统,二进制是计算机的基础。

b 可以用 2 的次幂唯一表达:b=i=0nbi2i

加法与乘法,这里必须采用位运算。

定义 带余除法

对于给定的任意整数 a,b 其中 b>0,存在唯一的整数对 q,r 使得 a=bq+rr[0,b1]

证明用到良序原则(自然数的非空子集必然存在一个最小元素)。

先考虑存在性的证明。构造集合 S={abk|kZ,abk0},那么必然存在一个最小元 r,满足 aqb=r,a=qb+r。接下来证明 r<b,这很显然,因为如果 rb 则存在 r=rb0,也在集合 S 中,矛盾。接下来证明唯一性。由于 r0,所以只存在唯一的 r 满足 0r<b,那么也就只有唯一 q

定义 整除

a=bq 则称 b 整除 a,记作 ba。称 ba 的约数。整除具有传递性。

定义 最大公约数

对于两个数a,b,如果 dadb 则称其为公约数,其中最大的称为最大公约数记作 gcd(a,b)

定义 素数 只有 1 和其本身两个约数的数是素数,除了 1 和其本身还有其他的约数的数是合数,10 既不是素数也不是合数。

定理 唯一分解定理

任意正整数 n=i=1mpici,其中 pi 是质数。

习题#

1.6 任意整数的平方都形如 3k 或者 3k+1kZ

  • n=3q 时,n2=9q2 满足。
  • n=3q+1 时,n2=9q2+6q+1 满足。
  • n=3q1 时,n2=9q26q+1 满足。

1.8 证明 2n1,n>2 不是完全平方数。

由于 2n 一定是偶数,所以 x 一定是奇数 x=2k+1,所以 x2=4(k2+k)+1,显然没有解。

2 欧几里德算法及相关定理#

2.1 欧几里德算法#

定理 欧几里德算法

gcd(a,b)=gcd(b,amodb)。其中 amodb 代表 a 带余除法 b 所得的余数 r

证明之前先证一个引理,gcd(a,b)=gcd(b,ab)。假设 gcd(a,b)=d,设 a=ud,b=vd,gcd(u,v)=1,那么 ab=d(uv)。所以 gcd(b,ab)=dgcd(v,uv)=d。这里用反证法,如果 gcd(v,uv)=k>1,那么设 v=kx,(uv)=kyu=k(x+y),那么 gcd(u,v) 至少是 k。这和 gcd(u,v)=1 矛盾。

有了这个引理之后,根据带余除法的定义,就可以证明该定理了。

命题 对任意两个正整数 ab,有 gcd(2a1,2b1)=2gcd(a,b)1

扩展一下,变成 gcd(xa1,xb1)=xgcd(a,b)1

不妨设 ab,那么 xa1(xb1)=xb(xab1)=k(xr1)。故 gcd(xa1,xb1)=gcd(xb1,xb(xab1))。又因为相邻两个数互质,所以 gcd(xb1,xb(xab1))=gcd(xb1,xab1)。所以有了 gcd(xa1,xb1)=gcd(xb1,xamodb1)

终止的条件即为 ba,这时即为 xgcd(a,b)1

定理 欧几里德算法的时间复杂度是 O(logmax(a,b))

只需证明,ab,amodb<a2。分类讨论即可。

2.2 Stein 算法#

又称作二进制欧几里德算法。

  • 如果 2a,2b,则 gcd(a,b)=2gcd(a/2,b/2)
  • 如果 2a,2b,则 gcd(a,b)=gcd(a/2,b)
  • 如果 2a,2b,则 gcd(a,b)=gcd(a,b/2)
  • 如果 2a,2b,则 gcd(a,b)=gcd(ab,b)

这个算法中没有乘法和除法以及取模运算,只用到了位运算,效率得到了提升。

然而实现表明这个算法是最慢的,因为它的理论常数会是最大的。即欧几里得算法准确复杂度应该是 Θ(log(max(a,b)gcd(a,b)))

2.3 扩展欧几里德算法#

定理 裴蜀定理 方程 ax+by=gcd(a,b) 存在整数解 x=s,y=t

纯数学证明可以构造集合并继续应用良序定理,这里给出扩展欧几里德算法的直接构造。

gcd(a,b)=d,欧几里得算法的最后一步会得到两个数 d,0,然后令s=1,t=0,就得到了一组解d+0=d。假设当前我们要求一组解满足 ax+by=d,且我们已经得到了一组解满足 bs+(amodb)t=d。拆开来得到 at+b(sabt)=d,所以我们可以令s=t,t=sabbt 而得到一组解s,t。一直这么推下去就可以得到最后的答案。

其实这是一个高斯消元的过程。

同样我们有二进制版本的扩展欧几里得算法。区别在于 a 是偶数的时候用 x=xb+12 而不是直接 x=x/2

Luogu5656 【模板】二元一次不定方程 (exgcd)

求解二元一次不定方程 ax+by=c

题解

2.4 类欧几里德算法#

看这里

2.5 万能欧几里德算法#

习题#

更多命题的证明在习题

2.5 将裴蜀定理扩展到 n 元。

首先证明 gcd(a1,a2,,an)=gcd(gcd(a1,a2,,an1),an)。这个好像比较显然,那么就可以归纳证明了。

2.6 假设 ga1(modm)gb1(modm),请证明 ggcd(a,b)1(modm)

用之前的证明的命题,ggcd(a,b)1=gcd(ga1,gb1),由于 mga1,mgb1,所以 mgcd(ga1,gb1)=ggcd(a,b)1

2.7 对任意不同的正整数 mnFmFn 互素。其中 F 是费马数。

不妨设 mnFm=22m1+222n+122n+1122m。于是有 Fm2(modF)n,又 2Fn,所以 gcd(Fm,Fn)=1

2.12 证明:对任意整数 a,b,dgcd(ad,bd)=gcd(a,b)d

k=gcd(a,b),a=ku,b=kv,则 gcd(kdud,kdvd)=kdgcd(ud,vd)=kd

2.14 证明:对任意整数 n,2n+1,9n+4 互素。

构造裴蜀定理 (2)n+2n+1=1。同理 9(2n+1)2(9n+4)=1

2.16 对非零整数 ab,如果整数 m 满足:ambm,则称 mab 的公倍数,称最小的 m 为最小公倍数,记作 lcm(a,b)。证明 lcm(a,b)=abgcd(a,b)

d=gcd(a,b),则 a=dm,b=dn,则最小公倍数得是 m,n,d 的倍数又 (m,n)=1 所以lcm(a,b)=dmn

最大公约数和最小公倍数的唯一分解定理

a=pici,b=pidi,则 gcd(a,b)=pimin(a,b),lcm(a,b)=pimax(a,b)

CF510D Fox And Jumping

根据裴蜀定理,ax+by 能表示出所有数当且仅当 (a,b)=1。所以这变成了一个背包问题,设 dpi 代表选出的数的最大公约数为 i 的最小价值。最后的答案为 dp1

我们可以通过 dpj 来更新 dpgcd(j,c[i])。但 gcd 可能很大,不能直接枚举。所以我们可以用平衡树把所有已经可能出现的 gcd 存下来,然后用其来更新。

3 同余与模运算#

3.1 同余与模算术运算#

定义 同余

如果整数 a,bm 满足 m(ab),则称 a,bm 同余。记作 ab(modm)

命题 如果 a1b1(modm) 而且 a2b2(modm),则 a1±a2b1±b2(modm)a1a2b1b2(modm)

定理 消去律

a,b,cZm 是一个正整数。如果 gcd(c,m)=1,且 acbc(modm),则 ab(modm)

m(acbc)=c(ab),由于 gcd(c,m)=1 所以 mab

定义 乘法逆元

a1 满足 aa11(modm)。显然乘法逆元不一定存在。

NOIP2012 同余方程

求关于 x 的同余方程 ax1(modb) 的最小正整数解。

也就是 ax+by=1 的最小正整数解,运用扩展欧几里德算法即可。

1,2,np 的逆元。

p=ki+j,0j<i

ki+j0(modp)(ki+j)i1j10(modp)kj1+i10(modp)i1kj1(modp)i1(p[p/i])×(pmodi)1(modp)

定义 等价类

显然,同余关系是一种等价关系。如果有一种等价关系定义在集合 S 上,则该等价关系会把集合 S 划分为不相交的子
集,这些子集我们称为等价类。

定义 同余类

考虑整数集合 Z,设 m 为正整数,则模 m 的同余关系把 Z 划分为 m 个等价类。任意整数 a 所落在的等价类记为 [a]m,即所有模 m 同余的整数形成的集合。

m 同余类的集合记为 Z/mZ,该集合恰好有 m 个元素,即:Z/mZ=[0]m,[1]m,,[m1]m

定义 最小剩余系

Z/mZ 的所有最小非负代表元的集合 {0,1,2,,m1}

3.2 二进制补码#

考虑 1x2n11,那么对于负数 x,我们找到他其实就相当于找到一个 x+x0(mod2)n,而这个值就应该是 2nx。取反加一的操作,其实就相当于变成了 2nx

3.3 模指数运算#

定理 快速幂

abai=0nbi2i(modm),于是就可以在 O(logb) 的时间求出 a2iab

Luogu1593 因子和

我们知道对于一个整数 n,可以分解成 p1c1×p2c2××pmcm。其因子和应为 Πi=1mj=0cipj

答案为 Πi=1mj=0cibpj。我们发现 里面即为等比数列求和。直接套公式为s=pcib+11p1

习题#

3.4 给定整数 pN,且 pN。证明:对任意整数 x(xmodN)modp=xmodp

使用剩余系即可证明。

3.11 证明 acbc(modm)gcd(c,m)=d,则 ab(modmd)

用 3.4 的结论即可说明。

4 费马小定理和欧拉定理#

4.1 费马小定理#

引理 a×imodp,1a<i1,2,,p1 的一个置换。

反证法,若存在 aiaj(modp),则 ij(modp) 矛盾。

定理 费马小定理

ap11(modp)

i=1p1ii=1p1ai=ap1i=1p1i(modp)

所以显然有 ap11(modp)

4.2 欧拉函数#

1,2,,n 中与 n 互质的数的个数为 φ(n),即为欧拉函数。

命题p 是素数,φ(p)=p1,φ(pc)=pc1(p1)

显然 p 的剩余系里都与 p 互质,那么每 p 个数里会有 p1 个互质,于是 φ(pc)=pcp1p

定理 欧拉函数是积性函数,即 gcd(a,b)=1φ(ab)=φ(a)φ(b)

这其实可以用剩余系得思想来推,即证明 Z/aZ×Z/bZ 双射 Z/abZ。当然还可以用容斥的角度,类似上面的证明。还有一个比较巧妙的用了狄利克雷卷积,这里就不过论述。

命题 d|nφ(d)=n

{1,2,,n}=dn{x[1,n]Zgcd(x,n)=d},而后者刚好就是 φ(nd)

欧拉函数的应用

4.3 欧拉定理#

定理 欧拉定理

aφ(m)1(modm),gcd(a,m)=1

定义 简化剩余系

所有和 m 互质的数构成的剩余类的集合。

证明欧拉定理等价于证明简化剩余系 ×a 是简化剩余系的置换,其实和上面的反证法一模一样。

4.4 扩展欧拉定理#

定理 扩展欧拉定理

好复杂,但是可以根据欧拉定理记忆,先丢一个证明

习题#

六省联考2017 相逢是问候

给定 n,m,p,c 区间操作将每个数变成 cai,区间求和 modp

CF906D Power Tower

给定一个数列 w 和模数 p,每次询问一个区间 [l,r],求 wlwl+1wl+2...modp 的值。

这两个题的思路是类似的,就是运用 p=φ(p) 这个迭代过程是 O(logp) 次的这个结论,进行暴力迭代。

5 素数与素性测试#

5.1 素数#

本章重点介绍素数的性质。

定理 素数的整除性

a,bZp 为素数。如果 pab,则 papb

推论,n>0,nab,gcd(n,a)=1nb

定理 欧几里得定理

素数有无穷多。

使用反证法。设素数有穷,p1,p2,pn,构造 N=1+i=1npi,这个数显然不能被任意一个素数整除,那么也就不能被任意一个数整除,所以它也是一个素数。

定理 狄利克雷定理

存在无穷多的素数满足 pa(modm),其中 gcd(a,m)=1

定理 素数定理

π(x) 为素数个数,则 limxπ(x)x/lnx=1

威尔逊定理

p 是素数的充要条件是 (p1)!1(modp)

p=1 时显然不行。

先证充分性,对于2a<bp2a1b1,我们只需将a[2,p2]a1 配对即可,充分性得证。

必要性显然,若 p 为合数则取一个质因子即可。

5.2 素性测试#

本节讨论素性测试问题,即如何判断一个大的奇整数 n 是否素数。

5.2.1 试除法#

枚举 2n1,看看能不能整除 n。这是最暴力的方法,时间复杂度O(n)

对上述算法优化,约数是成对出现的所以枚举到 n 即可,时间复杂度 O(n)

5.2.2 Miller-Rabin 素性测试#

费马小定理的逆定理虽然不成立,但是可以勉强用来判断。而满足费马小定理的合数称作卡麦尔数。

定理 二次探测定理

如果 p 为质数,则 x21modpx[1,p1] 的解为 x=1,p1

我们可以在费马小定理后以此判断 a(p1)/2,a(p1)/4 是不是 1/1

已经证明若随机选择 k 个底数则这样做出现伪素数的概率不大于 14k

code

习题#

6 群与子群#

定义

对于一个集合 S 和一个二元运算 ,如果其满足:

  • 封闭性:a,bS,abS
  • 结合律:a(bc)=(ab)c
  • 存在单位元 e,使得 aS,ae=ea=a
  • 存在逆元,使得 aS,bS,ab=e。、

定义 阿贝尔群

如果群 G 满足 a,bG,ab=ba 则称 G 是交换群或阿贝尔群。

img

6.1 显然乘法逆元在 Z 不封闭,而在 Q,R 内封闭。

6.2 没啥好解释的。

6.3 0 不存在乘法逆元。

6.4 根据费马小定理的证明可知。

6.5 根据欧拉定理的证明可知。

6.6 显然矩阵乘法不满足交换律

6.7 n 阶置换群,单位元是 (1,2,,n)

6.8 即对于 a2+b2=1,c2+d2=1,(acbd)2+(ad+bc)2=1,这是显然的。

6.2 群的基本性质#

命题 6.1 群的单位元唯一

如果存在两个不同的单位元 e,eG,则有 ee=e=e,矛盾。

命题 6.2 逆元唯一

若存在 m,n,使得 gm=gn=e,则 n=ne=ngm=em=m

命题 6.3G 为群,任取 a,bG,则 (ab)1=b1a1

由逆元唯一和 b1a1 满足逆元定义可得。

命题 6.4G 是一个群,对所有群元 gG,(g1)1=g

g1 带入 gg1=eg1(g1)1=e(g1)1=e(g1)1=gg1(g1)1=ge=g 所以命题成立。

对于群来说消去律一定存在。

定义 单射、满射和双射

给定映射 ϕ:AB,如果对任意 ϕ(x1)=ϕ(x2),则有 x1=x2,则称映射 ϕ 为单射(Injection),也称为一一映射。如果对任意 yB,都存在 xA,使得 ϕ(x)=y,则称映射 ϕ 为满射(Surjection)。如果映射 ϕ 同时满足单射性和满射性,则称映射 ϕ 为双射(Bijection)。

定义 置换群

定义映射 λa:gag,显然映射后还是群 G。容易发现这种映射是双射。所有映射组成一个群 G 称为置换群。

G 是一个群且 gG,则 gn 代表进行 n 次群运算后的结果。

由于群的计算是有结合律的所以 gmgn=gm+n(gm)n=gmn(gh)n=(h1g1)n

定义 阶(order)

有限群的阶定义为群的大小,无限群的阶定义为无穷大。对于群内的一个元素 gG,有 gm=e,则称最小的 mg 的阶写作 ord(g),假设不存在则该群是无限群。

6.3 子群#

G 为群,HGH 在群 G 的操作上也满足群公理,则称 HG 的子群,记作 HG

子群判定定理

G 的非空子集 HG 的子群,当且仅当 H 同时满足以下条件:

  • a,bH,则 abH
  • aH,则 a1H

已知 a,a1 均在群内,则 aa1=e 在群内,证毕。

G 的非空子集 HG 的子群,当且仅当对任意 a,bHab1H

aa1=eH

有限群 G 的非空子集 HG 的子群,当且仅当 H 中的元素在群 G 的操作下满足封闭性。

在之前的引理下只用证明 aHa1H

由于 H 封闭且有限,所以构造序列 a,a2,a3, 一定存在两个数 ai=aj,则 aij=e,aij1=a1

中心化子

对于任意群元素 aG,定义 C(a)={gG:ag=ga}。集合 C(a) 被称为 a 的中心化子。

中心

Z(G)={gG,xG,gx=xg}

显然 Z(G)=aGC(a)

习题#

6.5 证明对任意偶数阶群 G,都存在 gGgeg2=e

显然存在 gG,g|G|=e,则存在 g=g|G|2 满足条件。

6.8 证明:任意群 \mathbb G 的两个子群的交集也是群 G 的子群。

对于群 G 的两个子群 H1,H2,其满足 a,bH1H2 都有 abH1&abH2abH1H2

6.9 证明或证伪:任意群 G 的两个子群的并集也是群 G 的子群。

随便举例一个,2ZZ,3ZZ,2+3=52Z3Z

6.10 G 是阿贝尔群,HKG 的子群。请证明 HK={hk:hH,kK} 是群 G 的子群。如果 G 不是阿贝尔群,结论是否依然成立?

6.11 设 G 是阿贝尔群,m 是任意整数,记 Gm={gm:gG}。请证明 GmG 的一个子群。

gm,gm 均在群 Gm 里,

6.12 设 G 是阿贝尔群,m 是任意整数,记 G[m]={gG:gm=e}。请证明 G[m]G 的一个子群。

6.13 证明:阿贝尔群 G 中的有限阶元素形成子群。该子群有一个特殊的名字,称为挠子群(Torsion Subgroup)。

6.14 给定一个正整数 n(可以是素数也可以是合数),将 Zn 中所有元素相加并模 n 会等于什么?请证明。

会等于 0s=dd[gcd(d,n)=1]=nφ(n)2。在 n=2 时显然成立,在 n>2φ(n) 会是偶数。

7 循环群#

7.1 循环群#

G 是群,gG,集合 g={gx:xZ}G 的子群。称 gg 的生成循环群,且这是最小的包含 g 的子群。

循环群的生成元就是阶等于群的阶的元素。

命题 所有的循环群都是阿贝尔群。

G=g 是循环群,生成元是 g,则所有元素 a,b 可表示成 gr,gs 则显然 ab=gr+s=gs+r=ba

命题 循环群 G=g 的每一个子群都是循环群。

7.2 原根#

对于群 Zn 是循环群,则称其生成元是 n 的原根。即满足阶是 φ(n) 的元素 g

满足 n=2,4,pa,2pa 的群 Zn 都是循环群。

7.3 循环群的属性#

设群 G=g 是阶为 n 的循环群。如果 h=gk,则 h 的阶为 nd,其中 d=gcd(k,n)

h 的阶是 m,则 gkm=e 于是 nkmndkdmndm,所以被 nd 整除的最小的数就是 nd

设群 G=g 是阶为 n 的循环群,则群 G 中恰有 φ(n) 个生成元。

据此我们可以推出模 n 乘法循环群有 φ(φ(n)) 个生成元。

7.4 找原根#

原根判定定理

g 是原根,当且仅当 i,gφ(n)pi1(modn),φ(n)=pici

最小原根的数量级是 n0.25 的。

习题#

7.4 如果 p 是素数,则 Zp 没有非平凡子群。

所有非单位元的元素的阶都是 p

7.5 设 n 为任意正整数,群 Zn 的生成元 r 满足以下条件 gcd(r,n)=1

设某个元素 g 的阶为 m,则 ngm,设 gcd(g,n)=ddngdmndm,这时 d1m 一定小于 n。又群 Zn 的阶是 n 所以一定只有 d=1 的元素才是生成元。

7.6 证明:如果群 G 没有非平凡子群,则群 G 是循环群。

任取 gG,ge,则 g 是循环群,又 G 没有非平凡子群所以 g=G 所以 G 是循环群。

7.8 证明:设 G 为任意群,且 gG。如果存在 m,nZ 使得 gm=1gn=1,则 ggcd(n,m)=1

gnm=1

8 陪集和拉格朗日定理#

8.1 陪集#

G 是群,H 是群 G 的子群。定义 H 的左陪集为 gH={gh:hH},其中,元素 gG 被称为代表元(Representative)。右陪集的定义类似:Hg={hg:hH}

证明费马小定理和欧拉定理。

ZpZ 的子群,那么 aZp 是一个陪集,而 aZp,这时 ah 会被吸收掉,于是 aZp=Zp

命题 8.1

G 是群,HG 的子群。任取 g1,g2G,则 g1H=g2H 当且仅当 g2g1H

g1h=g2hg2=g1hh1g1H

命题 8.2G 是群,HG 的子群。对任意的 gG,群 H 的阶与集合 gH 的阶相同。

h1h2gh1gh2

命题 8.3G 为群,HG 的子群。任取不同的两个群元素 g1,g2G,则 g1H=g2Hg1Hg2H=

g1h1=g2h2g2=g1h1h21g2h3=g1h1h21h3g2H=g1H

命题 8.4G 是群,H 是群 G 的子群。则子群 H 的所有左陪集将划分群 G

显然所有左陪集覆盖整个群 G

定义 指标

G 是群,H 是群 G 的子群。H 在群 G 上的指标(Index)定义为 HG 上不相同的左陪集的个数,并记为 [G:H]

8.2 拉格朗日定理#

G 是有限群,HG 的子群,则 |G||H|=[G:H]

证明由命题8.4显然得到。

推论 8.1 任意元素的阶 ord(g)|G|.

推论 8.2 素数阶群任意元素均为生成元。

由推论 8.1 可得。

推论 8.3HK 是有限群 G 的子群,并且 K<H<G,则 [G:K]=[G:H][H:K]

[G:K]=|G||H||H||K|=[G:H][H:K]

由拉格朗日定理可得,任意有限群 G 及其中元素 gG,有 g|G|=e。这样即可证明费马小定理和欧拉定理。

9 同构,同态与商群*#

9.1 同构#

给定群 (G,) 和群 (H,),如果存在一种双射 ϕ:GH,使得群操作得以保持,即:任取 a,bGϕ(ab)=ϕ(a)ϕ(b)。称群 G 与群 H 同构,并记为 GH。映射 ϕ 被称为一种同构映射。

同构是等价关系,将所有群分成若干等价类。

定理 9.1 凯莱定理

每一个群都同构于一个置换群。

10 中国剩余定理#

10.1 数论角度下的中国剩余定理#

物不知数问题:

有物不知其数,三三数之剩二,五五数之剩三,七七数之剩二。问物几何?
化解为线性方程组来表示即为:

{x2(mod3)x3(mod5)x2(mod7)

对于这个方程组我们可以对于每个方程构造一个数使得其模其它模数为0,模自己的模数为答案。

即对于上面的方程可以构造:

{a1=35a2=63a3=30

全部加起来即可得到一组特解 x0=128。然后得到通解 x=x0+Mt,tZM 是模数的最小公倍数。

可以得到最小正整数解:x=23

通过这个例子我们可以得到下面的结论:

如果模数两两互素,设 M=Πi=1nmi,Mi=Mmi。根据古人的智慧我们可以构造出一个特解:x0=i=1nbiMiMi1

容易证明这个特解是对的,然后我们就可以得到通解:x=x0+M

多数情况下模数不是质数,所以要用exgcd求逆元。

10.2 扩展中国剩余定理#

这个时候,我们就要每个方程轮流按顺序求解。

假设我们已经求出一个解x=x0+Mt,tZ。这个解满足前i1个方程。现在我们将其和第i个方程合并。

具体操作就是将其代入第i个方程来解出得到未知数t的一个通解。

设第i个方程是形如xa(modb)的。

解方程的过程:

x0+Mta(modb)

Mtax0(modb)

如果这个方程无解则整个方程组无解。

g=gcd(M,b)

解得:t=t0+bgt,tZ

带回 x 的解,得:x=x0+M(t0+bgt)=x0+Mt0+lcm(M,b)t,tZ

x0 更新为 x0+Mt0M更新为 lcm(M,b) 即可。

第一个方程显然可以得到一组解x=a+bt,tZ

10.3 群论角度下的中国剩余定理#

n=pqp,q 是质数。

ZnZp×ZqZnZp×Zq

11 二次剩余#

11.1 定义#

p 是奇素数,a 是与 p 互质的整数,询问是否存在 x 使得 x2a(modp)

11.2 勒让德符号和欧拉准则#

定义 amodp 的勒让德符号为:

(ap)={1(exist)1(not)0(pa)

欧拉准则:

(ap)ap12(modp)

(ap)=1 时,x,x2a,(x2)(p1)/21ap121

否则,a(p1)/210,又 0ap11=(a(p1)/21)(a(p1)/2+1)a(p1)/2+10

推论

p 是奇素数,则

(1p)={1(pmod4=1)1(pmod4=1)

二次互反律

(pq)(qp)=(1)(p1)(q1)4

11.3 Cipolla算法#

我们随机一个 n,令 ω=n2a。将数系扩充成 Fω,相当于复数里面的 i。我们反复随机 n,直到 n2a 不是一个二次剩余。然后大概可以证明这样的 n 的大概有一半,所以期望随机二次既可以随机出一个,然后答案就是 (n+ω)p+12

证明先考虑两个引理:

ωp=ωap+bp(a+b)p(modp)

第一个考虑拆成 (ω2)p12×ω,第二个考虑二项式定理展开。

然后这个证明就很显然了,直接暴力平方展开然后替换,最后是一个平方差的形式,相信大家都会。

显然 x,x 都是二次剩余。

提交记录

12 环和域#

12.1 环#

定义

R 是一种包含两种运算的代数结构,分别称为加法和乘法。则加法是一个交换群,乘法满足交换律且存在单位元。乘法满足加法分配律。如果乘法也满足交换律,则称环 R 是交换环。

如果乘法逆元存在且满足交换律,则称其为域 F

13 有限域#

咕咕咕了

14 大数分解算法#

质因数分解
唯一分解定理
所有的合数都可以分解成下面的形式:

p1c1×p2c2××pmcm,一般p1<p2<<pm

试除法
枚举所有质数,然后判断是不是其因数,如果是再一直除下去求有多少个这个质因数。

注意枚举的因数只用到根号n,因为不可能有两个大于根号n的素因子。除到最后如果n不是1就说明n存在一个大于根号n的质因子,加到n的质因子序列中即可。

所以筛质数一般也只用筛到根号n。

Pollard-Rho
一个随机化算法,比较迷。

具体就是找到n的一个约数d,然后分治去求d的质因数分解和n/d的质因数分解。再用MR判断一下当前n是不是质数,如果是直接加到质因数序列即可。

至于怎么找到这个d,这就是随机化的思想了。

先搞出来一个x,y在[1,n-1]的范围内随机找,然后每次令x=f(x),y=f(f(y)),直到1<gcd(|x-y|,n)<n就退出循环gcd(|x-y|,n)就是n的一个约数。

f(x)=(x+d)%n, d=random[1,n-1]。

这样做的时间复杂度大约是O(n^{\frac{1}{4}})

15 离散对数#

posted @   Semsue  阅读(138)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· AI与.NET技术实操系列(五):向量存储与相似性搜索在 .NET 中的实现
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
Title
点击右上角即可分享
微信分享提示
主题色彩