在数学分析的广阔领域中，欧拉倒易公式（Euler's Reciprocal Formula）以其简洁的形式和深刻的内涵，成为连接离散与连续、对称与对偶关系的重要桥梁，该公式揭示了正整数幂次和与伯努利多项式之间的内在联系，不仅为级数求和提供了新工具，还在组合数学、数论和物理中有着广泛应用，本文将从欧拉倒易公式的基本形式出发，通过逐步推导，展现其从对称性观察到严格证明的全过程,并揭示其背后的数学思想。
欧拉倒易公式的基本形式
欧拉倒易公式通常表述为：对于正整数 ( n \geq 1 )，有
[
\sum{k=1}^{\infty} \frac{1}{k^n} + \sum{k=1}^{\infty} \frac{(-1)^{k+1}}{k^n} = 2 \sum{k=1}^{\lfloor n/2 \rfloor} \frac{1}{(2k)^n} + \begin{cases}
1 & \text{若 } n \text{ 为偶数}, \
0 & \text{若 } n \text{ 为奇数}.
\end{cases}
]
更简洁的版本可统一表示为：
[
\sum{k=1}^{\infty} \frac{1}{k^n} + \sum{k=1}^{\infty} \frac
( \zeta(n) = \sum{k=1}^{\infty} \frac{1}{k^n} ) 为黎曼 zeta 函数，( \delta{n \text{ 偶}} ) 当 ( n ) 为偶数时为 1，否则为 0。
这一公式的核心在于：将全正项级数与交错级数的和，转化为偶数项级数的加权和，体现了级数项的“对称性”与“对偶性”。
推导前的准备：伯努利多项式与生成函数
欧拉倒易公式的推导离不开伯努利多项式（Bernoulli Polynomials）这一工具，伯努利多项式 ( Bn(x) ) 是由雅各布·伯努利引入的一类特殊多项式，其生成函数为：
[
\frac{t e^{x t}}{e^t - 1} = \sum{n=0}^{\infty} B_n(x) \frac{t^n}{n!}, \quad |t| < 2\pi.
]
伯努利多项式具有以下重要性质（我们仅列出与推导相关的部分）：
- 积分关系：( \int_0^1 B_n(x) \, dx = 0 )（对 ( n \geq 1 )）；
- 递推关系：( Bn'(x) = n B{n-1}(x) )，且 ( B_n(0) = B_n(1) )（对 ( n \geq 2 )）；
- 对称性：( B_n(1-x) = (-1)^n B_n(x) )；
- 特殊值：伯努利数 ( B_n = Bn(0) )，且 ( B{2m+1} = 0 )（对 ( m \geq 1 )）。
黎曼 zeta 函数与伯努数的关系（欧拉公式）为：
[
\zeta(2m) = (-1)^{m+1} \frac{(2\pi)^{2m} B_{2m}}{2 (2m)!}, \quad m \geq 1.
]
推导步骤：从生成函数到级数分解
第一步：利用生成函数表达 ( \zeta(n) ) 与交错级数
考虑生成函数 ( \frac{t}{e^t - 1} = \sum_{n=0}^{\infty} Bn \frac{t^n}{n!} )，其展开式中 ( t^n ) 的系数与伯努利数直接相关，对于 ( \zeta(n) )，当 ( n \geq 2 ) 时，有：
[
\zeta(n) = \sum{k=1}^{\infty} \frac{1}{k^n} = \frac{1}{\Gamma(n)} \int_0^{\infty} \frac{t^{n-1}}{e^t - 1} \, dt,
]
这是通过积分变换将级数转化为积分的结果（利用 ( \frac{1}{k^n} = \frac{1}{\Gamma(n)} \int_0^{\infty} t^{n-1} e^{-k t} \, dt ) 并交换求和与积分顺序）。
类似地，交错级数 ( \eta(n) = \sum{k=1}^{\infty} \frac{(-1)^{k+1}}{k^n} )（称为 Dirichlet eta 函数）可表示为：
[
\eta(n) = \sum{k=1}^{\infty} \frac{(-1)^{k+1}}{k^n} = \frac{1}{\Gamma(n)} \int_0^{\infty} \frac{t^{n-1}}{e^t + 1} \, dt.
]
第二步：将全正项级数与交错级数相加
考虑 ( \zeta(n) + \eta(n) )：
[
\zeta(n) + \eta(n) = \sum{k=1}^{\infty} \frac{1}{k^n} + \sum{k=1}^{\infty} \frac{(-1)^{k+1}}{k^n} = \sum{k=1}^{\infty} \frac{1 + (-1)^{k+1}}{k^n}.
]
观察分子 ( 1 + (-1)^{k+1} )：当 ( k ) 为奇数时，( (-1)^{k+1} = 1 )，故分子为 2；当 ( k ) 为偶数时，( (-1)^{k+1} = -1 )，故分子为 0，级数仅保留奇数项：
[
\zeta(n) + \eta(n) = 2 \sum{m=1}^{\infty} \frac{1}{(2m-1)^n}.
]
第三步：将全正项级数拆分为偶数项与奇数项
全正项级数 ( \zeta(n) ) 可拆分为偶数项与奇数项：
[
\zeta(n) = \sum{k=1}^{\infty} \frac{1}{k^n} = \sum{m=1}^{\infty} \frac{1}{(2m)^n} + \sum{m=1}^{\infty} \frac{1}{(2m-1)^n} = 2^{-n} \zeta(n) + \sum{m=1}^{\infty} \frac{1}{(2m-1)^n}.
]
解得奇数项和为：
[
\sum_{m=1}^{\infty} \frac{1}{(2m-1)^n} = \zeta(n) - 2^{-n} \zeta(n) = (1 - 2^{-n}) \zeta(n).
]
第四步：代入第二步结果并整理
将第三步的奇数项和代入第二步的 ( \zeta(n) + \eta(n) ) 表达式：
[
\zeta(n) + \eta(n) = 2 \cdot (1 - 2^{-n}) \zeta(n) = 2 \zeta(n) - 2^{1-n} \zeta(n).
]
整理得：
[
\eta(n) = (2 - 2^{1-n} - 1) \zeta(n) = (1 - 2^{1-n}) \zeta(n).
]
但这只是中间结果，我们需要更精细的表达，以体现 ( n ) 的奇偶性差异。
第五步：利用伯努利多项式对称性修正奇偶性差异
回到伯努利多项式的对称性 ( B_n(1-x) = (-1)^n B_n(x) )，对于 ( x = 0 )，有 ( B_n(1
