第七章 同余方程

1、Roots of congruences(同余方程的根)

满足同余方程

f(x)=c0xn+c1xn1++cn0(modm)f(x)=c_{0} x^{n}+c_{1} x^{n-1}+\ldots+c_{n} \equiv 0(\bmod m)

xx 称为同余方程的 root(根),如果 aa 是同余方程的根,那么与 aa 同余的数都是方程的根,同余的根我们视为等价的,因此,当我们说一个同余方程有 ll 个根时,指的是方程有 ll 个不同余的根

根据分析学的成果,我们知道 nn 次代数方程在复数范围内有且仅有 nn 个根,且 nn 次多项式可以分解为 nn 个线性因子的乘积,我们自然而然的想到去探究同余方程有没有类似的性质,然而下面的几个例子告诉我们同余方程的性质并非这么简单

我们知道根据欧拉定理,方程

xp110(modp)x^{p-1}-1 \equiv 0(\bmod p)

p1p-1 个根,分别为 1,2,,p11,2, \ldots, p-1

而方程

x410(mod16)x^{4}-1 \equiv 0(\bmod 16)

88 个根,分别为 1,3,5,7,9,11,13,151,3,5,7,9,11,13,15

下面的方程方程则没有根

x420(mod16)x^{4}-2 \equiv 0(\bmod 16)

由此可见,同余方程的性质比代数方程要复杂得多


2、Integral polynomials and identical congruences(整多项式及其同余)

c0,c1,,cnc_{0}, c_{1}, \ldots, c_{n} 是整数,那么

c0xn+c1xn1++cnc_{0} x^{n}+c_{1} x^{n-1}+\cdots+c_{n}

称为 integral polynomial(整多项式),若

f(x)=r=0ncrxnr,g(x)=r=0ncrxnrf(x)=\sum_{r=0}^{n} c_{r} x^{n-r}, \quad g(x)=\sum_{r=0}^{n} c_{r}^{\prime} x^{n-r}

满足 crcr(modm)c_{r} \equiv c_{r}^{\prime}(\bmod m),那么我们称 f(x),g(x)f(x),g(x) 同余 (modm)(\bmod m),记作

f(x)g(x)(modm)f(x) \equiv g(x)(\bmod m)

若另有一个整多项式 h(x)h(x),那么

f(x)g(x)f(x)h(x)g(x)h(x)f(x) \equiv g(x) \rightarrow f(x) h(x) \equiv g(x) h(x)

Theorem 104
(1)若 pp 是素数,且
f(x)g(x)0(modp)f(x) g(x) \equiv 0(\bmod p)
那么 f(x)0f(x) \equiv 0g(x)0(modp)g(x) \equiv 0(\bmod p)
(2)更一般的,若
f(x)g(x)0(modpa),f(x)≢0(modp)f(x) g(x) \equiv 0\left(\bmod p^{a}\right),\quad f(x) \not\equiv 0(\bmod p)
那么 g(x)0(modpa)g(x) \equiv 0(\bmod p^{a})

(1)我们设 f1(x),g1(x)f_{1}(x),g_{1}(x) 分别是 f(x),g(x)f(x),g(x) 删掉被 pp 整除的项之后的结果

假设 f(x)≢0,g(x)≢0f(x) \not\equiv 0,g(x) \not\equiv 0,那么 f1(x),g1(x)f_{1}(x),g_{1}(x) 中的第一项系数不能被 pp 整除,因此 f1(x)g1(x)f_{1}(x)g_{1}(x) 的第一项系数也不能被 pp 整除

f(x)g(x)f1(x)g1(x)≢0(modp)f(x) g(x) \equiv f_{1}(x) g_{1}(x) \not\equiv 0(\bmod p)

(2)我们把 f(x)f(x) 中系数是 pp 的倍数的项删掉,g(x)g(x) 中系数是 pap^{a} 的倍数的项删掉,那么用同样的方法就可以证明,定理的这一部分将在第八章用到

值得注意的是,若 f(x)g(x)f(x) \equiv g(x),那么 f(a)g(a)f(a)\equiv g(a) 对于所有的 aa 成立,但是反过来不一定成立,例如

apa(modp)a^{p} \equiv a(\bmod p)

对于所有的 aa 成立,但是

xpx(modp)x^{p} \equiv x(\bmod p)

显然不成立


3、Divisibility of polynomials(多项式的整除性)

如果存在一个整多项式 h(x)h(x) 使得

f(x)g(x)h(x)(modm)f(x) \equiv g(x) h(x)(\bmod m)

那么我们称 f(x)f(x) 能被 g(x)g(x) 整除 (modm)(\bmod m),记作

g(x)f(x)(modm)g(x)|f(x)(\bmod m)

Theorem 105
使得
(xa)f(x)(modm)(x-a)|f(x)(\bmod m)
成立的一个充要条件是
f(a)0(modm)f(a)\equiv 0(\bmod m)

(xa)f(x)(modm)(x-a)|f(x)(\bmod m),那么存在 h(x)h(x) 使得

f(x)(xa)h(x)(modm)f(x) \equiv(x-a) h(x)(\bmod m)

因此 f(a)0(modm)f(a)\equiv 0(\bmod m),必要性得证

f(a)0(modm)f(a)\equiv 0(\bmod m),那么

f(x)f(x)f(a)(modm)f(x) \equiv f(x)-f(a)(\bmod m)

f(x)=crxnrf(x)=\sum c_{r} x^{n-r},那么

f(x)f(a)=(xa)h(x)f(x)-f(a)=(x-a) h(x)

其中

h(x)=f(x)f(a)xa=cr(xnr1+xnr2a++anr1)h(x)=\frac{f(x)-f(a)}{x-a}=\sum c_{r}\left(x^{n-r-1}+x^{n-r-2} a+\cdots+a^{n-r-1}\right)

也就是说,存在 h(x)h(x) 使得

f(x)(xa)h(x)(modm)f(x)\equiv (x-a)h(x)(\bmod m)

(xa)f(x)(modm)(x-a)|f(x)(\bmod m),证毕


4、Roots of congruences to a prime modulus(模素数意义下同余方程的根)

Theorem 106
pp 是素数,且
f(x)g(x)h(x)(modp)f(x) \equiv g(x) h(x)(\bmod p)
那么 f(x)f(x) 的根是 g(x)g(x)h(x)h(x) 的根 (modm)(\bmod m)

aaf(x)f(x) 的根,那么

f(a)g(a)h(a)0(modp)f(a) \equiv g(a)h(a)\equiv 0(\bmod p)

因此 g(a)0g(a)\equiv 0h(a)0h(a)\equiv 0,因此 aag(x)g(x)h(x)h(x) 的根 (modm)(\bmod m)

这里,模数 pp 是素数的条件是必要的,例如

x2x24(x2)(x+2)(mod4)x^{2} \equiv x^{2}-4 \equiv(x-2)(x+2)(\bmod 4)

44x20(mod4)x^2\equiv 0(\bmod 4) 的根,但不是 x20(mod4)x-2\equiv 0(\bmod 4)x+20(mod4)x+2\equiv 0(\bmod 4) 的根

Theorem 107
f(x)f(x)nn 次多项式,且有多于 nn 个根 (modp)(\bmod p),那么
f(x)0(modp)f(x)\equiv 0(\bmod p)

这个定理等价于 nn 次非零多项式至多有 nn 个根,我们用数学归纳法证明

n=1n=1 时,根据 Theorem 57 显然成立

假设当 f(x)f(x) 次数小于 nn 时成立,那么当次数等于 nn 时,设 f(a)0f(a)\equiv 0

f(x)(xa)g(x)(modp)f(x)\equiv (x-a)g(x)(\bmod p)

根据 Theorem 106f(x)f(x) 的根要么是 aa,要么是 g(x)g(x) 的根,若 f(x)f(x) 有多于 nn 个根,那么 g(x)g(x) 有多于 n1n-1 个根,因此

g(x)0f(x)0(modp)g(x)\equiv 0\Rightarrow f(x)\equiv 0(\bmod p)

这里,模数 pp 是素数的条件也是必要的,例如

x410(mod16)x^{4}-1\equiv 0(\bmod 16)

88 个根,事实上刚才的过程已经证明了下面的定理

Theorem 108
f(x)f(x) 所有的根为 a1,a2,,an(modp)a_{1}, a_{2}, \ldots, a_{n}(\bmod p),那么
f(x)c0(xa1)(xa2)(xan)(modp)f(x) \equiv c_{0}\left(x-a_{1}\right)\left(x-a_{2}\right) \ldots\left(x-a_{n}\right)(\bmod p)


5、Some applications of the general theorems(扩展理论的应用)

我们知道根据费马定理,当 d=p1d=p-1 时,同余方程

xd1(modp)x^{d} \equiv 1(\bmod p)

p1p-1 个根,我们接下来证明当 d(p1)d|(p-1) 时,结论也是正确的

Theorem 109
pp 为素数且 d(p1)d|(p-1),那么同余方程
xd1(modp)x^{d} \equiv 1(\bmod p)
dd 个根

我们有 xp11=(xd1)g(x)x^{p-1}-1=\left(x^{d}-1\right) g(x),其中

g(x)=xp1d+xp12d++xd+1g(x)=x^{p-1-d}+x^{p-1-2 d}+\cdots+x^{d}+1

由于 xp110x^{p-1}-1\equiv 0p1p-1 个根,根据 Theorem 106g(x)0g(x)\equiv 0 最多有 p1dp-1-d 个根,因此 xd10x^{d}-1\equiv 0 至少有 dd 个根

结合 xd10x^{d}-1\equiv 0 至多有 dd 个根,即可证明有且仅有 dd 个根

在这 dd 个根中,有些以 dd 为阶,另外的一些则是以更小的 p1p-1 的因子为阶,接下来的定理将给出以 dd 为阶的根的个数

Theorem 110
在同余方程
xd1(modp)x^{d}\equiv 1(\bmod p)
的所有根中,有 ϕ(d)\phi(d) 个根以 dd 为阶

ψ(d)\psi(d) 表示以 dd 为阶的根的个数,那么

d(p1)ψ(d)=p1\sum_{d | (p-1)} \psi(d)=p-1

根据 Theorem 63,有

d(p1)ϕ(d)=p1\sum_{d | (p-1)} \phi(d)=p-1

那么接下来只需要证明 ψ(d)ϕ(d)\psi(d)\leq \phi(d) 即可证明 ψ(d)=ϕ(d)\psi(d)=\phi(d)

ψ(d)>0\psi(d)>0,那么至少有一个根 ffdd 为阶,考虑下面 dd 个数

1,f,f2,,fd11,f,f^{2},\cdots,f^{d-1}

它们都是方程的根,且两两不同余(这是因为若存在 h<h<dh'<h<d 使得 fhfhf^{h}\equiv f^{h'},则 fk1f^{k}\equiv 1,其中 k=hh<dk=h-h'<d,与阶的定义不符)

根据 Theorem 109,它们是方程的所有根

fhf^{h}dd 为阶,那么必须满足 (h,d)=1(h,d)=1,这是因为若 k=(h,d)>1k=(h,d)>1,则

(fh)d/k=(fd)h/k1\left(f^{h}\right)^{d / k}=\left(f^{d}\right)^{h / k} \equiv 1

fhf^{h} 的阶必将小于 dd,因此 hh 必须在与 dd 互质的 ϕ(d)\phi(d) 个数中,即 ψ(d)ϕ(d)\psi(d)\leq\phi(d),证毕

在定理中取 d=p1d=p-1,我们得到推论:ppϕ(p1)\phi(p-1) 个原根

在证明过程中,我们已经证明了下面的定理

Theorem 111
pp 是奇素数,gg 是模 pp 的原根,那么
1,g,g2,,gp21,g,g^2,\cdots,g^{p-2}
两两不同余 (modp)(\bmod p)

我们知道

f(x)=xp11f(x)=x^{p-1}-1

的根为 1,2,,p11,2,\cdots,p-1,结合 Theorem 108

Theorem 112
pp 是素数,那么
xp11(x1)(x2)(xp+1)(modp)x^{p-1}-1 \equiv(x-1)(x-2) \ldots(x-p+1)(\bmod p)

比较两边的常数项,我们就能得到威尔逊定理

更一般的,比较 x,x2,,xp2x,x^{2},\cdots,x^{p-2} 项系数,可以得到

Theorem 113
pp 是奇素数,且 1l<p11\leq l<p-1,设 AlA_{l} 表示从集合
{1,2,,p1}\{1,2,\cdots,p-1\}
中选出 ll 个不同的数的乘积之和,则 Al0(modp)A_{l}\equiv 0(\bmod p)

我们可以用 Theorem 112 证明 Theorem 76,设

n=rps(r1,0s<p)n=r p-s(r \geqslant 1,0 \leqslant s<p)

那么有

(p+n1n)=(rps+p1)!(rps)!(p1)!=(rps+1)(rps+2)(rps+p1)(p1)!\begin{aligned}\binom{p+n-1}{n} &=\frac{(r p-s+p-1) !}{(r p-s) !(p-1) !} \\ &=\frac{(r p-s+1)(r p-s+2) \cdots(r p-s+p-1)}{(p-1) !} \end{aligned}

设结果为 ii,根据威尔逊定理,

(rps+1)(rps+2)(rps+p1)=(p1)!ii(modp)(r p-s+1)(r p-s+2) \dots(r p-s+p-1)=(p-1) ! i \equiv-i(\bmod p)

根据 Theorem 112,左边和下式同余

(s1)(s2)(sp+1)sp11(modp)(s-1)(s-2) \cdots(s-p+1) \equiv s^{p-1}-1(\bmod p)

因此

(p+n1n)(sp11)(modp)\binom{p+n-1}{n}\equiv -(s^{p-1}-1)(\bmod p)

pnp|n 时,s=0s=0,对应项系数为 11,否则,对应项系数为 00


6、Lagrange's proof of Fermat's and Wilson's theorem(拉格朗日对费马定理与威尔逊定理的证明)

我们对 Theorem 112 的证明是基于费马定理和 Theorem 108

而定理的发现者拉格朗日本人的证明更加简单粗暴,并且他的证明过程中包含了对费马定理及威尔逊定理的证明

我们设 pp 是奇素数,则

(x1)(x2)(xp+1)=xp1A1xp2++Ap1(x-1)(x-2) \cdots(x-p+1)=x^{p-1}-A_{1} x^{p-2}+\cdots+A_{p-1}

其中 A1,A2,A_{1},A_{2},\ldots 的定义与 Theorem 113 中相同,我们两边同乘 xx,然后把 xx 换成 x1x-1 得到

(x1)pA1(x1)p1++Ap1(x1)=(x1)(x2)(xp)=(xp)(xp1A1xp2++Ap1)\begin{aligned}(x-1)^{p}-A_{1}(x-1)^{p-1}+\cdots+A_{p-1}(x-1)&=(x-1)(x-2) \cdots(x-p) \\&=(x-p)\left(x^{p-1}-A_{1} x^{p-2}+\cdots+A_{p-1}\right) \end{aligned}

对比各项系数得到

(p1)+A1=p+A1(p2)+(p11)A1+A2=pA1+A2(p3)+(p12)A1+(p21)A2+A3=pA2+A3\begin{aligned}\binom{p}{1}+A_{1}&=p+A_{1} \\ \binom{p}{2}+\binom{p-1}{1}A_{1}+A_{2}&=pA_{1}+A_{2} \\ \binom{p}{3}+\binom{p-1}{2}A_{1}+\binom{p-2}{1}A_{2}+A{3}&=pA_{2}+A_{3}\\ &\cdots\end{aligned}

第一个式子是恒等式,其余的式子可以化成

A1=(p2)2A2=(p3)+(p12)A13A3=(p4)+(p13)A1+(p22)A2(p1)Ap1=1+A1+A2++Ap2\begin{aligned}A_{1}&=\binom{p}{2} \\ 2A_{2}&=\binom{p}{3}+\binom{p-1}{2}A_{1}\\ 3A_{3}&=\binom{p}{4}+\binom{p-1}{3}A_{1}+\binom{p-2}{2}A_{2}\\ &\cdots\\(p-1)A_{p-1}&=1+A_{1}+A_{2}+\cdots+A_{p-2} \end{aligned}

因此我们得到

pA1,pA2,,pAp2p\left|A_{1}, \quad p\right| A_{2}, \quad \ldots, \quad p | A_{p-2}

(p1)Ap11(modp)Ap11(modp)(p-1) A_{p-1} \equiv 1(\bmod p) \Rightarrow A_{p-1}\equiv -1(\bmod p)

这就证明了威尔逊定理以及 Theorem 112,费马定理随之得证


7、The residue of {12(p1)}!\{\frac{1}{2}(p-1)\}!

pp 是奇素数,且

ϖ=12(p1)\varpi=\frac{1}{2}(p-1)

根据威尔逊定理

(p1)!=1.212(p1){p12(p1)}{p12(p3)}(p1)(1)ϖ(ϖ!)2(modp)\begin{aligned}(p-1) ! &=1.2 \ldots \frac{1}{2}(p-1)\left\{p-\frac{1}{2}(p-1)\right\}\left\{p-\frac{1}{2}(p-3)\right\} \ldots(p-1) \\ & \equiv(-1)^{\varpi}(\varpi !)^{2}(\bmod p) \end{aligned}

因此

(ϖ!)2(1)ϖ1(modp)(\varpi !)^{2} \equiv(-1)^{\varpi-1}(\bmod p)

p=4n+1p=4n+1 时,

(ϖ!)21(modp)(\varpi !)^{2} \equiv -1(\bmod p)

由于 1-1 是模 pp 的二次剩余,ϖ!\varpi! 是方程 x2=1(modp)x^2=-1(\bmod p) 的某个根

p=4n+3p=4n+3 时,

(ϖ!)21ϖ!±1(modp)(\varpi !)^{2} \equiv 1 \Rightarrow \varpi ! \equiv \pm 1(\bmod p)

由于 1-1 是模 pp 的二次非剩余,上式中 ±1\pm 1 的取值取决于 ϖ!\varpi ! 是否为模 pp 的二次剩余

ϖ!\varpi ! 是小于 12p\frac{1}{2}p 的正整数的乘积,根据 Theorem 85 可得

Theorem 114
pp4n+34n+3 型素数,那么
{12(p1)}!(1)ν(modp)\left\{\frac{1}{2}(p-1)\right\} ! \equiv(-1)^{\nu}(\bmod p)
其中 ν\nu 表示小于 12p\frac{1}{2}p 的二次非剩余的个数


8、A theorem of Wolstenholme(沃斯滕霍姆定理)

根据 Theorem 113,我们知道

1+12+13++1p11+\frac{1}{2}+\frac{1}{3}+\cdots+\frac{1}{p-1}

的分子能被 pp 整除,因为这个分子就是 Ap2A_{p-2},更一般的,

Theorem 115
pp 是大于 33 的素数,那么
1+12+13++1p11+\frac{1}{2}+\frac{1}{3}+\cdots+\frac{1}{p-1}
的分子能被 p2p^2 整除

我们可以用另一种方式来表述这个定理,若 (i,m)=1(i,m)=1,方程

ix1(modm)i x \equiv 1(\bmod m)

仅有一个根 xx,且 i,xi,x 互为数论共轭,这对数论共轭比较特殊,它们的乘积是 11,我们称这样的数论共轭为乘法逆元,记作

1i=i\frac{1}{i}=\overline{i}

更一般的,我们用 ab\frac{a}{b} 表示方程 axb(modm)ax\equiv b(\bmod m) 的根

引入这样的记号后,我们就得到了 Wolstenholme's theorem

Theorem 116(Wolstenholme's Theorem)
p>3p>3,且 1i\frac{1}{i} 表示 ii 的乘法逆元 (modp2)(\bmod p^{2}),那么
1+12+13++1p10(modp2)1+\frac{1}{2}+\frac{1}{3}+\cdots+\frac{1}{p-1} \equiv 0\left(\bmod p^{2}\right)

为了解释乘法逆元符号的应用,我们先来证明

1+12+13++1p10(modp)1+\frac{1}{2}+\frac{1}{3}+\cdots+\frac{1}{p-1} \equiv 0\left(\bmod p\right)

我们知道,对于 0<i<p0<i<p,有

i1i1,(pi)1pi1(modp)i \cdot \frac{1}{i} \equiv 1, \quad(p-i) \frac{1}{p-i} \equiv 1(\bmod p)

因此

i(1i+1pi)i1i(pi)1pi0(modp)i\left(\frac{1}{i}+\frac{1}{p-i}\right) \equiv i \cdot \frac{1}{i}-(p-i) \frac{1}{p-i} \equiv 0(\bmod p)

1i+1pi0(modp)\frac{1}{i}+\frac{1}{p-i} \equiv 0(\bmod p)

然后我们求个和就能证明这个问题了

接下来我们来说明一下 Theorem 115Theorem 116 是等价的

对于 0<x<p0<x<p,设 x\overline{x} 表示 xx 的乘法逆元 (modp2)(\bmod p^{2}),那么

x(p1)!=xx(p1)!x(p1)x(modp2)\overline{x}(p-1) !=x \overline{x} \frac{(p-1) !}{x} \equiv \frac{(p-1)}{x}\left(\bmod p^{2}\right)

因此

(p1)!(1+2++p1)(p1)!(1+12++1p1)(modp2)\begin{array}{l}{(p-1) !(\overline{1}+\overline{2}+\cdots+\overline{p-1})} \\ {\quad \equiv(p-1) !\left(1+\frac{1}{2}+\cdots+\frac{1}{p-1}\right)\left(\bmod p^{2}\right)}\end{array}

等式右边就是分子,等价性即证,接下来证明这个定理

(x1)(x2)(xp+1)=xp1A1xp2++Ap1(x-1)(x-2) \ldots(x-p+1)=x^{p-1}-A_{1} x^{p-2}+\ldots+A_{p-1}

x=px=p,得到

(p1)!=pp1A1pp2+Ap2p+Ap1(p-1) !=p^{p-1}-A_{1} p^{p-2}+\ldots-A_{p-2} p+A_{p-1}

我们知道 (p1)!=Ap1(p-1)!=A_{p-1},因此

pp2A1pp3++Ap3pAp2=0p^{p-2}-A_{1} p^{p-3}+\ldots+A_{p-3 p}-A_{p-2}=0

由于 p>3p>3pA1,pA2,,pAp3p\left|A_{1}, p\right| A_{2}, \ldots, p | A_{p-3},得到 p2Ap2p^{2} | A_{p-2},即

p2(p1)!(1+12++1p1)p^{2} |(p-1) !\left(1+\frac{1}{2}+\ldots+\frac{1}{p-1}\right)

证毕,另外还有一个定理

Theorem 117
p>3p>3,那么
Cp=1+122++1(p1)20(modp)C_{p}=1+\frac{1}{2^{2}}+\ldots+\frac{1}{(p-1)^{2}}\equiv 0(\bmod p)

事实上,CpC_{p} 的分子就是 Ap222Ap1Ap3A_{p-2}^{2}-2 A_{p-1} A_{p-3},显然成立


9、The theorem of von Staudt(施陶特定理)

我们最后介绍施陶特定理,这个定理与 Bernoulli's number(伯努利数) 有关

伯努利数通常定义为

xex1=112x+B12!x2B24!x4+B36!x6\frac{x}{e^{x}-1}=1-\frac{1}{2} x+\frac{B_{1}}{2 !} x^{2}-\frac{B_{2}}{4 !} x^{4}+\frac{B_{3}}{6 !} x^{6}-\dots

该展开式在 x<2π|x|<2\pi 时收敛,我们发现,另外一种写法更为简洁

xex1=β0+β11!x+β22!x2+β33!x3+\frac{x}{e^{x}-1}=\beta_{0}+\frac{\beta_{1}}{1 !} x+\frac{\beta_{2}}{2 !} x^{2}+\frac{\beta_{3}}{3 !} x^{3}+\ldots

其中 β0=1,β1=12\beta_{0}=1, \beta_{1}=-\frac{1}{2},且

β2k=(1)k1Bk,β2k+1=0(k1)\beta_{2 k}=(-1)^{k-1} B_{k}, \quad \beta_{2 k+1}=0 \quad(k \geqslant 1)

伯努利数与 mk\sum m^kEuler-Maclaurin sum-formula(欧拉-麦克劳林求和公式) 密切相关,事实上,

1k+2k++(n1)k=r=0k1k+1r(kr)nk+1rβr1^{k}+2^{k}+\ldots+(n-1)^{k}=\sum_{r=0}^{k} \frac{1}{k+1-r}\binom{k}{r} n^{k+1-r} \beta_{r}

这是因为左边的式子是

k!x(1+ex+e2x++e(n1)x)k ! x\left(1+e^{x}+e^{2 x}+\ldots+e^{(n-1) x}\right)

xk+1x^{k+1} 项系数,而

k!x(1+ex+e2x++e(n1)x)=k!x1enx1ex=k!xex1(enx1)=k!(1+β11!x+β22!x2+)(nx+n2x22!+)\begin{array}{l}{k! x\left(1+e^{x}+e^{2 x}+\ldots+e^{(n-1) x}\right)} \\ {\quad=k ! x \frac{1-e^{n x}}{1-e^{x}}=k ! \frac{x}{e^{x}-1}\left(e^{n x}-1\right)} \\ {\quad=k !\left(1+\frac{\beta_{1}}{1 !} x+\frac{\beta_{2}}{2 !} x^{2}+\cdots\right)\left(nx+\frac{n^{2} x^{2}}{2 !}+\ldots\right)}\end{array}

计算对应的 xk+1x^{k+1} 系数就得到了上面的结论

施陶特定理基于 BkB_{k} 的小数部分

Theorem 118(Staudt's Theorem)
k1k\geq 1,那么
(1)kBk(p1)2k1p(mod1)(-1)^{k} B_{k} \equiv \sum_{(p-1)|2k} \frac{1}{p}(\bmod 1)

例如,若 k=1k=1,那么和式中 p=2,3p=2,3,而 B1=56B_{1}=-\frac{5}{6},显然成立

我们把 BkB_{k} 换成 βk\beta_{k},即

βk+(p1)k1p=i\beta_{k}+\sum_{(p-1) | k} \frac{1}{p}=i

其中 k=1,2,4,6,k=1,2,4,6,\ldots,且 ii 是一个整数,如果定义 ϵk(p)\epsilon_{k}(p)

ϵk(p)={1,(p1)k0,(p1)k\epsilon_{k}(p)=\left\{\begin{matrix}1,&(p-1) \mid k \\ 0,&(p-1) \nmid k\end{matrix}\right.

那么式子就可以写成

βk+ϵk(p)p=i\beta_{k}+\sum \frac{\epsilon_{k}(p)}{p}=i

事实上,施陶特定理说明了伯努利数的分母中没有平方因子


10、Proof of von Staudt's theorem(施陶特定理的证明)

我们首先来证明一个引理

Theorem 119
m=1p1mkϵk(p)(modp)\sum_{m=1}^{p-1} m^{k} \equiv-\epsilon_{k}(p)(\bmod p)

(p1)k(p-1)\mid k,那么根据费马定理,mk1(modp)m^k\equiv 1(\bmod p),因此

mkp11ϵk(p)(modp)\sum m^{k} \equiv p-1 \equiv-1 \equiv-\epsilon_{k}(p)(\bmod p)

(p1)k(p-1)\nmid k,设 gg 是模 pp 的原根,那么根据 Theorem 88

gk≢1(modp)g^{k} \not\equiv 1(\bmod p)

由于集合 {g,2g,,(p1)g}\{g,2g,\ldots,(p-1)g\}{1,2,,p1}\{1,2,\ldots,p-1\} 同余 (modp)(\bmod p)

(mg)kmk(gk1)mk0(modp)\sum(m g)^{k} \equiv \sum m^{k}\Rightarrow \left(g^{k}-1\right) \sum m^{k} \equiv 0(\bmod p)

这样就得到了

mk0=ϵk(p)(modp)\sum m^{k} \equiv 0=-\epsilon_{k}(p)(\bmod p)

现在我们用数学归纳法来证明 Theorem 118

我们知道定理中 k=1,2,4,6,k=1,2,4,6,\ldots,显然当 k=1,2k=1,2 时定理成立

假设定理在 L={1,2,4,,k2}L=\{1,2,4,\ldots,k-2\} 的范围内成立,我们需要证明定理对于 kk 成立

ϖ\varpi 是任意素数,根据 Theorem 119

ϵk(ϖ)+r=0k1k+1r(kr)ϖk+1rβr0(modϖ)\epsilon_{k}(\varpi)+\sum_{r=0}^{k} \frac{1}{k+1-r}\binom{k}{r} \varpi^{k+1-r} \beta_{r} \equiv 0(\bmod \varpi)

由于 βk1=0\beta_{k-1}=0,得到

βk+ϵk(ϖ)ϖ+r=0k21k+1r(kr)ϖk1r(ϖβr)0(mod1)\beta_{k}+\frac{\epsilon_{k}(\varpi)}{\varpi}+\sum_{r=0}^{k-2} \frac{1}{k+1-r}\binom{k}{r} \varpi^{k-1-r}\left(\varpi \beta_{r}\right) \equiv 0(\bmod 1)

接下来我们考虑

uk,r=1k+1r(kr)ϖk1r(ϖβr)u_{k, r}=\frac{1}{k+1-r}\binom{k}{r} \varpi^{k-1-r}\left(\varpi \beta_{r}\right)

的分母能否被 ϖ\varpi 整除

r∉Lr\not\in L,则 βr=10\beta_{r}=1或0,否则根据归纳假设,βr\beta_{r} 的分母中不含平方因子,即 (ϖβr)(\varpi \beta_{r}) 的分母不能被 ϖ\varpi 整除,而 (kr)\binom{k}{r} 是整数,因此 uk,ru_{k,r} 的分母能被 ϖ\varpi 整除当且仅当

ϖk1rk+1r=ϖs1s+1\frac{\varpi^{k-1-r}}{k+1-r}=\frac{\varpi^{s-1}}{s+1}

的分母能被 ϖ\varpi 整除,那就需要满足

s+1ϖss+1 \geqslant \varpi^{s}

但是我们知道 s=kr2s=k-r\geq 2,因此

s+1<2sϖss+1<2^{s} \leqslant \varpi^{s}

所以 uk,ru_{k,r} 的分母不能被 ϖ\varpi 整除,故

βk+ϵk(ϖ)ϖ=akbk\beta_{k}+\frac{\epsilon_{k}(\varpi)}{\varpi}=\frac{a_{k}}{b_{k}}

其中 ϖbk\varpi\nmid b_{k},且

ϵk(p)p(pϖ)\frac{\epsilon_{k}(p)}{p}(p \neq \varpi)

显然有相同的形式,故

βk+ϵk(p)p=AkBk\beta_{k}+\sum \frac{\epsilon_{k}(p)}{p}=\frac{A_{k}}{B_{k}}

其中 BkϖB_{k}\nmid \varpi,由于 ϖ\varpi 是任意素数,所以 Bk=1B_{k}=1,证毕

特别的,若 kk3n+13n+1 型素数,那么在 (p1)2k(p-1)|2k 条件下

p=2,3,k+1,2k+1p=2,3,k+1,2k+1

但是 k+1k+1 是偶数,2k+1=6n+32k+1=6n+3 不是素数,因此 p=2,3p=2,3,故

Theorem 120
kk3n+13n+1 型素数,那么
Bk16(mod1)B_{k}\equiv \frac{1}{6}(\bmod 1)

事实上,这个定理可以推广到对于给定的 kk,有无穷多个 BlB_{l} 使得

BlBk(mod1)B_{l}\equiv B_{k}(\bmod 1)

我们需要 DirichletsTheorem 15 来证明它