实作正 17 边形
一、分组
按照高斯的步骤,我们首先选择一个 17 的原根 g=3,按照原根的幂次给除了 1 以外的 16 个本原单位根排序和分组,相同的颜色代表分在同一组中:
ζg0 ζg1 ζg2 ζg3 ζg4 ζg5 ζg6 ζg7 ζg8 ζg9 ζg10 ζg11 ζg12 ζg13 ζg14 ζg15
e=1
ζ1 ζ3 ζ9 ζ10 ζ13 ζ5 ζ15 ζ11 ζ16 ζ14 ζ8 ζ7 ζ4 ζ12 ζ2 ζ6
e=2 (P0,P1)
ζ1 ζ3 ζ9 ζ10 ζ13 ζ5 ζ15 ζ11 ζ16 ζ14 ζ8 ζ7 ζ4 ζ12 ζ2 ζ6
e=4 (Q0,Q1,Q2,Q3) P0=Q0+Q2,P1=Q1+Q3
ζ1 ζ3 ζ9 ζ10 ζ13 ζ5 ζ15 ζ11 ζ16 ζ14 ζ8 ζ7 ζ4 ζ12 ζ2 ζ6
e=8
ζ1 ζ3 ζ9 ζ10 ζ13 ζ5 ζ15 ζ11 ζ16 ζ14 ζ8 ζ7 ζ4 ζ12 ζ2 ζ6
我们为 e=8 这一层(每组 2 个元素)补充上组名,按照相同的拆分规则,它们被命名为 R0 到 R7。其中我们最关心的是包含 ζ1 的 R0 组:
R0=ζ1+ζ16
R4=ζ13+ζ4
上一层的 Q0 正好是由它们合并而成的:Q0=R0+R4。
下面我们自顶向下,依次建立二次方程并求解。
二、求解
1. 第一层拆分:求解 P0 和 P1
在 e=2 这一层,所有的 16 个本原单位根被分为了 P0 和 P1 两组,每组 8 个元素。
根据上一章的结论,所有本原根之和为 −1,所以:
P0+P1=−1
接下来计算它们的乘积 P0×P1。
根据多项式乘法,这两组相乘会产生 8×8=64 个项。根据我们上一章证明的“代数对称性”(性质2、3),这 64 个项必然能完美还原成上层已知组的正整数倍线性组合,而上层只有一组,所以必然是所有根出现相同的次数,因为 64÷16=4,所以 64 项之和一定是所有本原单位根之和的 4 倍,即:
P0×P1=4(ζ1+ζ2+⋯+ζ16)=4(−1)=−4
现在,我们得到了一个完美的二次方程:
X2+X−4=0
根据求根公式,解得 X=2−1±17。因为 P0 中包含了 2cos(172π) 这样较大的正数,通过简单的数值估算可知 P0>P1,所以我们取正根:
P0=2−1+17
P1=2−1−17
2. 第二层拆分:求解 Q0 和 Q1
在 e=4 这一层,我们要把 P0 拆分为 Q0 和 Q2。
已知:
Q0+Q2=P0
计算它们的乘积:Q0=(ζ1+ζ13+ζ16+ζ4), Q2=(ζ9+ζ15+ζ8+ζ2)。它们相乘产生 16 项,指数相加后(模 17),恰好不重不漏地构成了 1∼16 各一次!因此:
Q0×Q2=ζ1+ζ2+⋯+ζ16=−1
构造二次方程:
X2−P0X−1=0
由于 Q0 包含 ζ1,估算 Q0>Q2,取正根:
Q0=2P0+P02+4
同理,对于 P1 拆分出来的 Q1 和 Q3,也有 Q1+Q3=P1 且 Q1×Q3=−1。我们同样解得:
Q1=2P1+P12+4
代入第一步求出的 P0 和 P1:
Q0=4−1+17+34−217
Q1=4−1−17+34+217
3. 第三层拆分:求解 R0
这是最激动人心的一步。在 e=8 这一层,Q0 被拆分为了 R0 和 R4。已知:
R0+R4=Q0
我们来计算它们的乘积 R0×R4:
(ζ1+ζ16)×(ζ13+ζ4)=ζ14+ζ5+ζ12+ζ3
这 4 个项加起来是什么?请回头看一眼前面的分组结果,它恰好就是 Q1!因此:
R0×R4=Q1
构造出最后一个二次方程:
X2−Q0X+Q1=0
因为 R0=ζ1+ζ16=2cos(172π),这是正十七边形中心角对应的长度,估算 R0>R4。我们取正根:
R0=2Q0+Q02−4Q1
4. 走向终点的最简表达式
我们离真相只差最后一步了。因为
R0=ζ1+ζ−1=2cos(172π)
所以,正十七边形作图的核心坐标值,就是:
cos(172π)=2R0=4Q0+Q02−4Q1
我们将前面得到的 Q0 和 Q1 的精确表达式代入其中。经过一些代数上的整理(将分母统一提取为 16),就得到了最终的表达式:
cos(172π)=−161+1617+1634−217+8117+317−34−217−234+217
三、实作正 5 边形
5=22+1 也是一个费马素数,它的求解过程相对短小精悍,可以再演示一次。我们的目标是求解 cos(52π)=cos(72∘)
1. 寻找原根和分组
对于正五边形,费马素数 p=5
除 1 以外的 4 个本原单位根为:ζ,ζ2,ζ3,ζ4。
我们寻找模 5 的原根,发现 g=2 满足条件:
20≡1(mod5)
21≡2(mod5)
22≡4(mod5)
23≡3(mod5)
我们将本原单位根按照原根 g=2 的幂次顺序排列,并进行分组:
e=1
ζ1 ζ2 ζ4 ζ3
e=2(P0,P1)
ζ1 ζ2 ζ4 ζ3
可以看到,4 个根被完美分成了两组:
P0=ζ1+ζ4
P1=ζ2+ζ3
2. 计算和与积
x5−1=(x−1)(x4+x3+x2+x+1)=(x−1)(x−ζ1)(x−ζ2)(x−ζ3)(x−ζ4)
两边约去 (x−1):
x4+x3+x2+x+1=(x−ζ1)(x−ζ2)(x−ζ3)(x−ζ4)
比较 x3 项的系数知:
ζ1+ζ2+ζ3+ζ4=−1
即所有本原单位根之和为 −1,所以:
P0+P1=−1
接下来计算这两组的乘积 P0×P1:
P0×P1=(ζ1+ζ4)(ζ2+ζ3)
展开得到:
P0×P1=ζ3+ζ4+ζ6+ζ7
因为 ζ5=1,所以 ζ6=ζ1,ζ7=ζ2。代入后得到:
P0×P1=ζ3+ζ4+ζ1+ζ2=−1
3. 构造方程与求解
现在我们有了和与积,直接构造二次方程:
X2−(P0+P1)X+(P0×P1)=0
X2+X−1=0
使用求根公式解得:
X=2−1±5
因为 P0=ζ1+ζ4=ζ1+ζ−1=2cos(52π)=2cos(72∘)>0。我们取正根:
P0=2−1+5
由于 cos(52π)=2P0,我们直接得出:
cos(52π)=4−1+5
仅仅只用了一次拆分、解了一个二次方程,我们就严谨地得出了 cos(72∘) 的精确根式。