Skip to content

算法Algorithm

简单代数扩张的精确运算

Arithmetic in a simple algebraic extension · Polynomial quotient field arithmetic

用唯一的低次余式表示代数数,并以多项式 Euclid 回代给出可逐项核验的逆元证书及非单位诊断。

形式陈述 ​

不用小数也能除以一个代数数 ​

设 α3=α+1。乘法中出现 α4 时,可以换成 α2+α;但如何求 (α2+α+1)−1?近似算出 α 的小数再做除法,只会得到近似值。这里把问题化为两条多项式恒等式,最后得到精确答案 2−α2。

本页的输入是域 F、首一不可约多项式 m∈F[x],以及以多项式表示的扩域元素。取商环

E=F[x]/(m),α=x+(m),d=deg⁡m.

若 F 的四则运算与相等判断可执行,下述过程就是算法;对抽象域,它仍是正确的代数构造,但不能凭域公理声称已有程序实现。

先统一表示,再做运算 ​

多项式带余除法给出唯一表达式 a=qm+r,其中 deg⁡r<d,或 r=0。因此 a(α)=r(α)。两个低次余式若表示同一元素,其差是 m 的倍数,却次数小于 d,只能为零。每个元素因而具有唯一坐标

a0+a1α+⋯+ad−1αd−1.

这是极小多项式所给幂基的计算版本。加法逐坐标进行,乘法先相乘再除以 m 取余;判断相等只需比较约简后的系数。约简须按同一个 m 进行,不可在运算中随意换成另一个以 α 为根、但次数更高的多项式。

求逆需要一个额外证书。给定 g,先取其模 m 的余式。余式为零时不能除;否则不可约性保证 gcd(g,m)=1。我们要实际找出 s,t∈F[x],使

sg+tm=1.

代入 α 后得到 s(α)g(α)=1。所以逆元就是 s 模 m 的余式;最后再乘回取余为 1,便能独立核验输出。

直觉

Euclid 为什么同时产出逆元 ​

沿欧几里得整环中的带余除法,从 r0=m,r1=g 出发反复除法:

ri−1=qiri+ri+1,ri+1=0 或 deg⁡ri+1<deg⁡ri.

非零余式的次数严格下降,所以有限步后停止。相邻两项的公因子与下一对相同,因此最后非零项是 gcd 的一个非零常数倍。

为了把它写回输入,给每项保留 ri=Aim+Big。初值为

(A0,B0)=(1,0),(A1,B1)=(0,1),

每次更新

Ai+1=Ai−1−qiAi,Bi+1=Bi−1−qiBi.

这条关系由余式公式直接推出,所以每一步的表示都正确。如果末个非零余式为常数 c≠0,把全部系数除以 c,得到 1=(Ai/c)m+(Bi/c)g,逆元取 Bi/c。这里归一化不可漏掉:末项若是 2,原回代式只证明乘积为 2,还不是逆元。

例子与边界

一个完整的三次除法 ​

取 F=Q、m=x3−x−1。它的有理根只能是 ±1,而两处取值都是 −1,因此三次式不可约。令 g=x2+x+1,做两轮除法:

m=(x−1)g−x,g=(−x−1)(−x)+1.

第二式给 1=g+(x+1)(−x);再用第一式中的 −x=m−(x−1)g 回代,得到

1=(2−x2)g+(x+1)m.

所以

(α2+α+1)−1=2−α2.

乘回也可以不用信任 Euclid 的过程。展开

(2−x2)(x2+x+1)=−x4−x3+x2+2x+2.

用 x3≡x+1、x4≡x2+x 模 m 约简,结果恰为 1。求逆与验逆是不同任务;验逆只需一次乘法和一次取余。

可约模多项式不会让所有元素失去逆元 ​

若把不可约的 m 换成任意首一正次数多项式 h,唯一低次余式、加法和乘法仍然成立;但商环未必是域。此时 g+(h) 可逆的准确条件是

gcd(g,h)=1.

充分性仍由 Bézout 证书给出。必要性则来自:若 sg≡1(modh),便有 sg+th=1,任何公共因子都必须整除 1。

例如在 Q[x]/(x2−1) 中,x−1 非零,却与同样非零的 x+1 相乘为零;其 gcd 为 x−1,不能求逆。另一方面,x 是单位,因为 x2≡1。因此检测到模多项式可约,意味着“不能保证每个非零元素可逆”,并不意味着每个非零元素都不可逆。

同一扩域中的逆元也不一定留在整数系数生成的子环中。在 Q(2) 内,1/2 是 2 的逆元,但不属于 Z[2]。本页运算的系数域是 F;不要将 F[α] 与某个整系数子环混为一谈。

推论与应用

运算成本与一次迁移检查 ​

把一次 F 中的加、减、乘、除计为一个域运算,且输入已经约简到次数小于 d。朴素加法需要 O(d) 个域运算;朴素乘法加上长除法取余需要 O(d2) 个域运算。Euclid 至多进行 d 轮非零余式下降,系数回代多项式的次数也不超过 d,因此按每轮至多 O(d2) 估算,得到保守的 O(d3) 求逆上界。这里没有声称这是最锐界。

这个计数把域运算当作原子操作。在 F=Q 中,分子分母的位数会增长,真实位复杂度还需计入整数运算与约分成本;若原输入次数远大于 d,最开始的约简成本也要另算。

自测:在同一个可约商环 Q[x]/(x2−1) 中,x+2 能求逆吗?能,因为

(x+2)(2−x)=4−x2≡3(modx2−1).

故逆元为 (2−x)/3。这里既展示了可约商环中的单位,也说明回代所得非零常数必须先归一化。

回到四次根式的运算证书 ​

已知 θ=2+3 的极小多项式为 m=x4−10x2+1,则

x(10x−x3)+m=1.

因此 θ−1=10θ−θ3。原有根式恢复公式遂可以完全写成幂基坐标:

2=θ3−9θ2,3=11θ−θ32.

这些不是重新证明极小多项式;不可约性与次数已在原例中完成。这里获得的是统一的算法接口:给定任意 g(θ)≠0,都能用同一种 Euclid 过程求逆,而不必为每个分母单独猜一个根式共轭。

参考资料
  • J. S. Milne,Fields and Galois Theory,v5.10,2022,Chapter 1,“Construction of some extensions”,条目 1.25 之前的 (a)–(e),印刷页16–17;Example 1.27 给另一组三次扩张逆元计算。
  • Keith Conrad,The Division Algorithm in Z and F[T],Theorem 1.2 与 §3,多项式带余除法与唯一性;本页使用的扩展系数递推由所列恒等式直接推出。
关系图谱15 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系