[论文解读] On Two Types of Harmonic Number Identities

M. J. Kronenburg|arXiv (Cornell University)|Feb 17, 2012

Scientific Research and Discoveries参考文献 6被引用 8

一句话总结

本文推导出两类有限调和数求和的闭式恒等式：$ H_k H_{n-k} $ 的乘积形式与 $ H_k^2 $ 的平方形式，且均以 $ k $ 的幂加权。通过生成函数、分部求和法以及涉及广义调和数和伯努利数的递归公式，作者提供了可计算的递推关系，从而实现对这些求和式在 $ p=5 $ 以内的高效计算与符号运算。

ABSTRACT

Two types of finite series of products of harmonic numbers involving nonnegative integer powers are evaluated, also yielding two other important harmonic number identities. The recursion formulas for these sums are derived, which are easily translated into a computer program.

研究动机与目标

• 推导有限求和 $ \sum_{k=0}^n k^p H_k H_{n-k} $ 与 $ \sum_{k=0}^n k^p H_k^2 $ 的显式闭式表达式，其中 $ H_k $ 为第 $ k $ 个调和数。
• 利用分部求和法及广义调和数与伯努利数的相关恒等式，建立这些求和式的递推公式。
• 提供一种系统化方法，用于计算高阶幂 $ p $ 下的这些求和，支持算法实现。
• 呈现 $ p = 5 $ 以内的符号表达式，并通过生成函数与系数提取方法加以验证。
• 将推导出的公式整合进 Mathematica 程序中，实现对这类调和数求和的自动求值。

提出的方法

• 使用生成函数：$ \sum_{k=0}^n H_k H_{n-k} = [z^n] \frac{1}{(1-z)^2} \left(\ln \frac{1}{1-z}\right)^2 $，其已知系数涉及第一类斯特林数。
• 通过恒等式 $ \sum_{k=a}^{b-1} x_k y_k = s_{b-1} y_b - \sum_{k=a}^{b-1} s_k (y_{k+1} - y_k) $ 应用分部求和法，其中 $ s_k = \sum_{i=a}^k x_i $，以推导递推关系。
• 采用恒等式 $ \sum_{k=0}^n k^p H_k^{(m)} = F(n,p,m) $，其中 $ F(n,p,m) $ 通过 $ H_n^{(-p)} $ 结合广义调和数与伯努利数构造而成。
• 利用符号运算与系数提取，分别推导出 $ \sum k^p H_k H_{n-k} $ 与 $ \sum k^p H_k^2 $ 的递推公式 (2.14) 与 (3.4)。
• 将递推关系转化为 Mathematica 程序，使用 `HarmSumR`、`HarmSumS` 与 `HarmSumT` 函数，并借助 `CoefficientArrays` 与 `Factor` 实现符号化简化。
• 通过计算 $ p = 0 $ 到 $ 5 $ 的显式表达式验证结果，其与已知恒等式一致，并给出新的闭式表达。

实验结果

研究问题

• RQ1 求和 $ \sum_{k=0}^n H_k H_{n-k} $ 的闭式表达式为何？如何将其推广至包含 $ k $ 的幂的形式？
• RQ2 如何利用调和数恒等式与生成函数，递归计算求和 $ \sum_{k=0}^n k^p H_k H_{n-k} $？
• RQ3 求和 $ \sum_{k=0}^n k^p H_k^2 $ 的结构为何？其与广义调和数及伯努利数有何关联？
• RQ4 能否为 $ \sum k^p H_k H_{n-k} $ 与 $ \sum k^p H_k^2 $ 推导出系统化的递推关系？如何实现算法化计算？
• RQ5 $ p = 5 $ 以内的这些求和的显式符号表达式为何？当 $ p = 0 $ 时，它们如何退化为已知恒等式？

主要发现

• 求和 $ \sum_{k=0}^n H_k H_{n-k} $ 的结果为 $ (n+1)[(H_{n+1} - 1)^2 - H_{n+1}^{(2)} + 1] $，该闭式恒等式通过生成函数与系数提取方法得到验证。
• 当 $ p = 1 $ 时，$ \sum_{k=0}^n k H_k H_{n-k} = H_n^{(-1)}[(H_{n+1} - 1)^2 - H_{n+1}^{(2)} + 1] $，其中 $ H_n^{(-1)} = \frac{1}{2}n(n+1) $，表明其直接依赖于负一阶调和数。
• 当 $ p = 2 $ 时，求和 $ \sum_{k=0}^n k^2 H_k H_{n-k} $ 的表达式为 $ H_n^{(-2)}[H_{n+1}^2 - H_{n+1}^{(2)}] - \frac{1}{18}n(n+1)(13n+5)H_{n+1} + \frac{1}{108}n(n+1)(71n+37) $，其中 $ H_n^{(-2)} = \frac{1}{6}n(n+1)(2n+1) $。
• 求和 $ \sum_{k=0}^n H_k^2 $ 的结果为 $ (n+1)H_{n+1}^2 - (2n+3)H_{n+1} + 2(n+1) $，该为已知恒等式，通过分部求和与下标平移推导得出。
• 当 $ p = 3 $ 时，$ \sum_{k=0}^n k^3 H_k^2 = H_n^{(-3)}H_{n+1}^2 - \frac{1}{24}n(n+1)(n+2)(3n+1)H_{n+1} + \frac{1}{288}n(n+1)(9n^2 + 13n + 14) $，其中 $ H_n^{(-3)} = \frac{1}{4}n^2(n+1)^2 $。
• 本文提供了一个可运行的 Mathematica 程序，可计算至 $ p = 5 $ 的这些求和，通过递归求值与 `CoefficientArrays` 和 `Factor` 实现的符号化简化。