🔢 mathematics

General formulas for a class of Euler sums

本文提出了一种易于实现的算法，利用偏微分函数和多项式函数将特定有理函数与调和数乘积的欧拉和转化为部分分式分解问题，并给出了分母含一至两项且幂次不超过 3 的显式闭式解公式。

原作者： David H Bailey, Ross McPhedran, Bruno Salvy

发布于 2026-04-06

📖 1 分钟阅读🧠 深度阅读

原作者： David H Bailey, Ross McPhedran, Bruno Salvy

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

这篇论文就像是一份**“数学无穷级数翻译指南”**。

想象一下，你面前有一堆极其复杂的数学算式，它们长得像是一串串无穷无尽的珠子（级数），中间还夹杂着一种叫“调和数”（ $H_k$ ）的累加项。数学家们（就像这篇论文的作者）发现，虽然这些算式看起来像是一团乱麻，但如果你用特定的“魔法钥匙”去打开它们，它们其实都能变成非常简洁、优雅的公式。

下面我用几个生活中的比喻来解释这篇论文的核心内容：

1. 什么是“欧拉和”？（那串乱麻）

想象你在玩一个无限长的多米诺骨牌游戏。

骨牌：每一个骨牌代表一个数字。
规则：第 $k$ 个骨牌的价值，不仅取决于它自己的编号，还取决于前面所有骨牌价值的总和（这就是“调和数” $H_k$ ）。
问题：如果你把这一整串无穷无尽的骨牌价值加起来，结果是多少？
- 有些简单的情况，结果是一个漂亮的数字（比如 $\pi$ 的几次方）。
- 但大多数情况，尤其是分母很复杂的时候（比如 $(2k+1)^3(3k+1)^2$ ），直接算简直是不可能的，就像试图数清大海里有多少滴水。

2. 作者做了什么？（制造“翻译机”）

作者 David H. Bailey、Ross McPhedran 和 Bruno Salvy 发明了一套通用的“翻译算法”。

以前的做法：遇到一个新的复杂算式，数学家得像个侦探一样，一个个去猜、去试，或者花几年时间推导，看看能不能凑出一个公式。
现在的做法：作者写了一个计算机程序（算法），只要你把那个复杂的算式（分母是多项式的形式）扔进去，程序就能自动帮你把它“翻译”成另一种语言。

3. 翻译成了什么？（从“乱麻”到“乐高积木”）

这个程序把复杂的无穷求和，翻译成了由**“特殊函数”**组成的公式。

比喻：
- 原来的算式像是一团乱糟糟的毛线球，你很难看清它的结构。
- 翻译后的公式像是乐高积木。虽然积木本身（也就是论文里提到的“双伽马函数” $\psi$ 和“多伽马函数” $\psi^{(n)}$ ）也是特殊的零件，但它们有固定的形状和规则。
- 一旦你有了这些标准积木，无论原来的毛线球多乱，你都能用这些积木拼出它的最终形状。

4. 他们是怎么做到的？（“留数”魔法）

论文的核心技术叫**“留数计算”**（Residue Computation）。

比喻：想象你在一个巨大的迷宫里找出口。
- 传统的求和方法是沿着迷宫的墙壁一步步走（直接计算每一项），这太慢了。
- 作者的方法是：他们发现迷宫里只有几个特定的“魔法点”（极点）。只要计算出这几个点的“魔法能量”（留数），就能直接算出整个迷宫的出口在哪里，完全不需要走完整个迷宫。
- 这个“魔法能量”的计算，正好对应了那些特殊的“乐高积木”（双伽马函数）。

5. 为什么这很重要？（不仅仅是为了好看）

你可能会问：“把算式变复杂了（从求和变成了函数公式），有什么意义？”

速度：直接算那个无穷级数，就像是用勺子一勺一勺地舀干大海，需要算几亿项才能得到几个小数位，慢得要死。
效率：用作者提供的公式，就像是用抽水机。计算机算那些“特殊函数”的速度极快，能瞬间给出几百位精度的答案。
验证：作者用超级计算机算出了成千上万个例子，发现这些公式在几百种情况下都完美吻合。这就像他们造了一台机器，证明了“只要输入这种类型的算式，输出一定符合这个规律”。

6. 总结

这篇论文就像是一本**“数学食谱大全”**。

以前，厨师（数学家）遇到一道新菜（复杂的欧拉和），得自己摸索做法。
现在，作者提供了一套标准化的烹饪流程（算法）和预制调料包（双伽马函数公式）。
只要你的菜符合“分母比分子复杂至少两阶”这个条件，你就可以直接套用这个流程，瞬间得到这道菜的“终极味道”（精确的解析解）。

这不仅让数学研究变得更快，也让那些原本看起来深不可测的无穷级数，变得像搭积木一样有章可循。

这是一份关于论文《一类欧拉和的通用公式》（General formulas for a class of Euler sums）的详细技术总结。该论文由 David H. Bailey、Ross McPhedran 和 Bruno Salvy 撰写，主要探讨了如何计算包含调和数（Harmonic Numbers）的无穷级数。

1. 研究问题 (Problem)

欧拉和（Euler Sums） 是一类重要的无穷级数，形式为 $\sum_{k=1}^{\infty} R(k)H_k$ ，其中：

$H_k = 1 + 1/2 + \dots + 1/k$ 是第 $k$ 个调和数。
$R(k)$ 是一个有理函数（两个多项式的商）。

核心挑战：
许多此类级数具有优雅的解析解（通常涉及黎曼 $\zeta$ 函数、 $\pi$ 的幂次等），但直接求和非常困难。特别是当分母包含多个线性因子（混合分母）或具有较高幂次时，手动推导闭式解极其繁琐且容易出错。

本文目标：
提出一种易于实现的算法，用于计算形如 $\sum_{k=1}^{\infty} R(k)H_k$ 的显式闭式解。要求 $R(k)$ 的分母次数至少比分子次数大 2（确保级数收敛），且分母在整数处无极点（除可能的 $k=0$ 外）。

2. 方法论 (Methodology)

论文的核心方法基于复变函数中的留数定理（Residue Computation），结合Digamma（ $\psi$ ）函数和Polygamma（ $\psi^{(n)}$ ）函数的性质。

2.1 理论基础

Digamma 函数性质： 利用 $\psi(z) = \Gamma'(z)/\Gamma(z)$ 与调和数的关系： $\psi(n+1) = H_n - \gamma$ （ $\gamma$ 为欧拉常数）。
极点展开： 在负整数附近， $\psi(-s) + \gamma$ 具有特定的级数展开形式，涉及 $\zeta(k)$ 。
留数公式（命题 1）：
对于满足条件的有理函数 $R(s)$ ，级数和可以表示为 $R(s)$ 所有极点 $\alpha$ 处的留数之和：
$\sum_{k=1}^{\infty} R(k)H_k = \frac{1}{2} \sum_{\alpha \in \text{poles of } R} \text{Res}_{s=\alpha} \left[ R'(s)(\psi(-s) + \gamma) - R(s)(\psi(-s) + \gamma)^2 \right]$
该公式将无穷级数求和转化为有限个复平面上极点的留数计算。

2.2 算法实现 (Algorithm HarmonicSum)

作者设计了一个计算机代数系统（CAS）友好的算法：

部分分式分解： 将 $R(s)$ 分解为形如 $\frac{C_{i,j}}{(s+t_j)^i}$ 的项。
符号计算： 不需要显式求出极点的数值，而是通过多项式的因式分解（在代数数域上）进行符号处理。
泰勒/洛朗展开： 在极点 $\alpha$ 处，对 $R(s)$ 、 $R'(s)$ 以及 $\psi(-s)+\gamma$ 进行泰勒展开。
提取系数： 计算展开式中 $(s-\alpha)^{-1}$ 项的系数（即留数），并求和。
通用公式推导： 基于上述算法，推导出了针对单极点、高阶极点以及混合分母（如 $(m_1k+n_1)^{p_1}(m_2k+n_2)^{p_2}$ ）的通用解析公式。

2.3 数值验证

为了验证解析公式的正确性，作者采用了高精度的数值计算方法：

直接计算级数的前 $10^8$ 项。
利用欧拉 - 麦克劳林求和公式（Euler-Maclaurin）处理尾部。
精度达到约 290 位小数。
使用 PSLQ 整数关系算法 将高精度数值与假设的解析形式（包含 $\pi, \zeta, \gamma, \psi$ 等）进行匹配，从而确定有理系数。

3. 主要贡献与结果 (Key Contributions & Results)

3.1 通用解析公式

论文给出了基于 Digamma 和 Polygamma 函数的通用公式：

单极点情况： 对于 $\sum \frac{H_k}{(k+t)^p}$ ，给出了用 $T(t, p)$ 表示的闭式解（定理 2）。
混合分母情况： 对于 $\sum \frac{H_k}{\prod (m_i k + n_i)^{p_i}}$ ，通过部分分式分解，将其转化为上述单极点公式的线性组合（定理 3）。
带分子幂次的情况： 扩展到了 $\sum \frac{k^q H_k}{(mk+n)^p}$ 的形式（推论 2）。

3.2 具体公式库 (Catalog of Formulas)

作者利用 Maple 和 Mathematica 生成了一个包含数百个具体公式的目录（附录 B），涵盖了：

分母项数：1 项或 2 项。
幂次 $p_i$ ：最高到 3（部分公式扩展到 8）。
系数 $m_i$ ：最高到 4。
例如： $\sum \frac{H_k}{(2k+1)^3(3k+1)^2}$ 等复杂形式的显式解。

3.3 计算效率分析

论文指出，虽然解析公式将求和转化为了特殊函数的计算，但计算效率显著提升：

直接数值求和达到 280 位精度需要极长时间。
使用基于 Digamma/Polygamma 的解析公式计算，速度比直接求和快约 30,000 倍。
这使得高精度计算和验证成为可能。

4. 意义与影响 (Significance)

自动化与通用性： 提供了一种系统化的方法，将原本需要极高数学技巧才能手动推导的欧拉和求和问题，转化为可自动执行的算法步骤。
连接离散与连续： 建立了离散调和级数与连续特殊函数（ $\psi, \zeta$ ）之间的深刻联系，揭示了这些看似复杂的级数背后的代数结构。
高精度计算工具： 为数学物理、量子场论（微扰计算中常出现此类级数）以及数论研究提供了高效的计算工具。
验证与发现： 通过“数值实验 + 整数关系算法 + 解析证明”的闭环流程，不仅验证了已知结果，还发现并证明了数百个新的恒等式。
软件实现： 提供了 Maple 代码实现，使得其他研究者可以方便地计算特定形式的欧拉和。

总结

这篇文章不仅提供了一套强大的算法来计算一类广泛的欧拉和，还通过大量的数值实验和解析推导，建立了一个包含数百个具体公式的数据库。其核心创新在于利用留数定理将无穷级数求和转化为特殊函数的代数运算，极大地提高了计算效率和公式生成的自动化程度，对数学物理和实验数学领域具有重要价值。