[论文解读] On Bernoulli Numbers and Its Properties

Lin Cong|ArXiv.org|Aug 6, 2004

History and Theory of Mathematics参考文献 2被引用 4

一句话总结

这篇综述论文全面概述了伯努利数的历史发展、通过生成函数的现代定义，以及它们在三角函数和双曲函数的级数展开中的出现。论文详细阐述了关键性质，如克莱森-冯·施陶特定理、它们通过正则素数在费马最后定理中的作用，以及高效计算方法，使得能够计算出极大伯努利数，包括具有超过三百万位数的 $ B_{750,000} $。

ABSTRACT

In this survey paper, I first review the history of Bernoulli numbers, then examine the modern definition of Bernoulli numbers and the appearance of Bernoulli numbers in expansion of functions. I revisit some properties of Bernoulli numbers and the history of the computation of big Bernoulli numbers.

研究动机与目标

• 从幂的和的古老公式到现代应用，呈现伯努利数的历史与数学综述。
• 解释伯努利数在余切、双曲余切、正切和双曲正切函数的级数展开中的作用。
• 探讨深刻的数论性质，包括克莱森-冯·施陶特定理及其与正则素数和奇异素数的联系。
• 回顾计算进展，从欧拉的计算到现代高精度算法，实现对具有数百万位数的伯努利数的计算。
• 强调伯努利数在数论中的重要性，包括其在库默尔研究费马最后定理中的关键作用。

提出的方法

• 通过生成函数 $ \frac{x}{e^x - 1} = \sum_{n=0}^\infty \frac{B_n x^n}{n!} $ 推导伯努利数，确立其定义与递推关系。
• 利用生成函数推导 $ \coth(x) $、$ \cot(x) $、$ \tanh(x) $ 和 $ \tan(x) $ 的级数展开，将其与偶数索引的伯努利数联系起来。
• 应用克莱森-冯·施陶特定理，将 $ B_{2k} $ 的分数部分计算为满足 $ p-1 \mid 2k $ 的素数 $ p $ 的倒数之和，证明 $ B_{2k} + \sum_{(p-1)\mid 2k} \frac{1}{p} \in \mathbb{Z} $。
• 利用 zeta 函数关系 $ \zeta(2k) = (-1)^{k-1} \frac{(2\pi)^{2k}}{2(2k)!} B_{2k} $，实现无需递归计算即可直接计算大伯努利数。
• 利用斯特林公式推导大 $ k $ 的渐近估计 $ B_{2k} \simeq (-1)^{k-1} 4 \left( \frac{k}{\pi e} \right)^{2k} \sqrt{\pi k} $。
• 回顾历史与现代计算里程碑，包括利用 zeta-伯努利关系在个人计算机上直接计算 $ B_{200,000} $ 和 $ B_{750,000} $。

实验结果

研究问题

• RQ1 伯努利数如何出现在 $ \sum_{k=1}^{n-1} k^r $ 的求和公式中？它们在将幂和表示为 $ n $ 的多项式时起什么作用？
• RQ2 伯努利数与 $ \cot(x) $、$ \coth(x) $、$ \tan(x) $ 和 $ \tanh(x) $ 的级数展开之间有何联系？
• RQ3 克莱森-冯·施陶特定理如何确定 $ B_{2k} $ 的分数部分？这对伯努利数的分母有何影响？
• RQ4 伯努利数通过库默尔的正则素数准则在费马最后定理的证明中起什么作用？
• RQ5 哪些计算技术使得能够直接且高效地计算出极大伯努利数，如 $ B_{750,000} $ ？

主要发现

• 生成函数 $ \frac{x}{e^x - 1} = \sum_{n=0}^\infty \frac{B_n x^n}{n!} $ 提供了伯努利数的现代定义，其中 $ B_0 = 1 $，$ B_1 = -\frac{1}{2} $，$ B_2 = \frac{1}{6} $，且对 $ n \geq 1 $ 有 $ B_{2n+1} = 0 $。
• $ x \cot(x) = \sum_{n=0}^\infty (-1)^n \frac{2B_{2n} (2x)^{2n-1}}{(2n)!} $ 的级数展开将伯努利数与三角函数联系起来。
• 克莱森-冯·施陶特定理指出 $ B_{2k} + \sum_{(p-1)\mid 2k} \frac{1}{p} $ 是整数，从而可精确确定 $ B_{2k} $ 的分数部分；例如，$ B_{16} \equiv \frac{463}{510} \pmod{1} $。
• 对于素数 $ p = 3n+1 $，$ B_{2p} $ 的分数部分恰好为 $ \frac{1}{6} $，如 $ B_{14} $、$ B_{26} $、$ B_{38} $ 等，这是由于符合条件的素数集合有限所致。
• 库默尔证明：若素数 $ p $ 是正则素数，即 $ p $ 不整除任何 $ B_{2k} $（$ k = 1, 2, \dots, \frac{p-3}{2} $）的分子，则费马最后定理对指数 $ p $ 成立；37 是第一个奇异素数，它整除 $ B_{32} $ 的分子。
• 利用 zeta 函数关系与直接计算，普鲁夫与菲计算出 $ B_{750,000} $，该数具有 3,391,993 位数，仅用 21 小时就在个人计算机上完成，创下最大显式计算的伯努利数纪录。