12+22+32+42+…+n²求和公式如何推导？

Question

在数学与算法分析中，求解前n个正整数平方和（即1² + 2² + 3² + … + n²）的公式为： $$ S = \frac{n(n+1)(2n+1)}{6} $$ 许多初学者在推导该公式时，常困惑于如何从直观累加过渡到闭合表达式。常见问题包括：为何不能直接类比等差数列求和方法？数学归纳法虽可验证公式正确性，但无法解释公式的来源。那么，**如何通过代数构造、组合恒等式或几何直观等方式，系统地推导出该平方和公式的原始来源？** 特别是利用立方差公式或递推关系构建辅助序列的方法，在实际推导中应如何操作并确保逻辑严密？

Answer 1

前n个正整数平方和公式的系统推导：从直觉到严密证明

在算法分析与数学建模中，求解前 $ n $ 个正整数的平方和：

$$ S_n = 1^2 + 2^2 + 3^2 + \cdots + n^2 $$

是一个经典问题。其闭合形式为：

$$ S_n = \frac{n(n+1)(2n+1)}{6} $$

尽管该公式可通过数学归纳法验证，但其“来源”常令人困惑。本文将从多个角度——代数构造、组合恒等式、几何直观以及递推关系——循序渐进地揭示该公式的原始推导路径。

1. 直观对比：为何不能类比等差数列求和？

• 等差数列求和公式 $ S = \frac{n(a_1 + a_n)}{2} $ 基于线性增长特性，项之间差值恒定。
• 而平方序列 $ a_k = k^2 $ 是非线性的，相邻项差值 $ (k+1)^2 - k^2 = 2k + 1 $ 随 $ k $ 变化，呈线性增长。
• 因此，无法通过首尾配对或平均值法直接推广。
• 这提示我们需要更高阶的方法：如构造辅助表达式、利用多项式拟合或恒等变形。

2. 方法一：代数构造法（立方差公式）

我们利用立方差恒等式构建递推关系。考虑以下恒等式：

$$ (k+1)^3 - k^3 = 3k^2 + 3k + 1 $$

对 $ k = 1 $ 到 $ n $ 求和：

$$ \sum_{k=1}^{n} \left[(k+1)^3 - k^3\right] = \sum_{k=1}^{n} (3k^2 + 3k + 1) $$

左边是望远镜求和：

$$ (n+1)^3 - 1^3 = 3\sum_{k=1}^{n}k^2 + 3\sum_{k=1}^{n}k + \sum_{k=1}^{n}1 $$

代入已知公式 $ \sum k = \frac{n(n+1)}{2}, \sum 1 = n $，得：

$$ (n+1)^3 - 1 = 3S_n + 3\cdot\frac{n(n+1)}{2} + n $$

整理后解出 $ S_n $：

$$ 3S_n = (n+1)^3 - 1 - \frac{3n(n+1)}{2} - n $$ $$ S_n = \frac{1}{3} \left[ (n+1)^3 - 1 - \frac{3n(n+1)}{2} - n \right] $$

经代数化简可得：

$$ S_n = \frac{n(n+1)(2n+1)}{6} $$

3. 方法二：递推关系与辅助序列构造

假设 $ S_n $ 是一个关于 $ n $ 的三次多项式（因平方和增长趋势接近积分 $ \int x^2 dx \sim n^3 $）：

$$ S_n = an^3 + bn^2 + cn + d $$

利用初始条件建立方程组：

n	S_n	方程
0	0	d = 0
1	1	a + b + c = 1
2	5	8a + 4b + 2c = 5
3	14	27a + 9b + 3c = 14
4	30	64a + 16b + 4c = 30
5	55	125a + 25b + 5c = 55
6	91	216a + 36b + 6c = 91
7	140	343a + 49b + 7c = 140
8	204	512a + 64b + 8c = 204
9	285	729a + 81b + 9c = 285

解此线性系统可得 $ a = \frac{1}{3}, b = \frac{1}{2}, c = \frac{1}{6}, d = 0 $，即：

$$ S_n = \frac{1}{3}n^3 + \frac{1}{2}n^2 + \frac{1}{6}n = \frac{n(n+1)(2n+1)}{6} $$

4. 方法三：组合恒等式视角

利用组合数恒等式：

$$ k^2 = 2\binom{k}{2} + \binom{k}{1} $$

因为：

$$ 2\binom{k}{2} = k(k-1),\quad \binom{k}{1} = k \Rightarrow 2\binom{k}{2} + \binom{k}{1} = k^2 - k + k = k^2 $$

于是：

$$ \sum_{k=1}^{n}k^2 = \sum_{k=1}^{n}\left(2\binom{k}{2} + \binom{k}{1}\right) = 2\sum_{k=1}^{n}\binom{k}{2} + \sum_{k=1}^{n}\binom{k}{1} $$

利用组合恒等式 $ \sum_{k=1}^{n}\binom{k}{r} = \binom{n+1}{r+1} $：

$$ \sum_{k=1}^{n}k^2 = 2\binom{n+1}{3} + \binom{n+1}{2} $$

展开计算：

$$ 2 \cdot \frac{(n+1)n(n-1)}{6} + \frac{(n+1)n}{2} = \frac{n(n+1)(n-1)}{3} + \frac{n(n+1)}{2} $$ $$ = n(n+1)\left(\frac{n-1}{3} + \frac{1}{2}\right) = n(n+1)\left(\frac{2n - 2 + 3}{6}\right) = \frac{n(n+1)(2n+1)}{6} $$

5. 方法四：几何直观与面积累加模型

考虑将 $ k^2 $ 视为边长为 $ k $ 的正方形面积。将这些正方形按层堆叠，形成阶梯状结构。

虽然无法完美拼接成矩形，但可通过三维类比：将 $ k^2 $ 看作第 $ k $ 层的“底面积”，构建一个近似金字塔体。

体积积分启发我们：$ \sum k^2 \approx \int_0^n x^2 dx = \frac{n^3}{3} $，说明结果应为三次函数。

进一步结合离散修正项，可引导出精确表达式。

6. 方法五：生成函数法（拓展视角）

定义生成函数 $ G(x) = \sum_{n=1}^{\infty} n^2 x^n $，可通过微分几何级数得到：

$$ \sum_{n=1}^{\infty} x^n = \frac{x}{1-x} \Rightarrow \sum_{n=1}^{\infty} n x^{n-1} = \frac{1}{(1-x)^2} \Rightarrow \sum_{n=1}^{\infty} n x^n = \frac{x}{(1-x)^2} $$

再次微分：

$$ \sum_{n=1}^{\infty} n^2 x^{n-1} = \frac{1+x}{(1-x)^3} \Rightarrow \sum_{n=1}^{\infty} n^2 x^n = \frac{x(1+x)}{(1-x)^3} $$

部分和可通过系数提取获得，但更适用于分析渐进行为。

7. 流程图：系统推导路径总览

```mermaid
graph TD
    A[目标: 推导 Σk² 公式] --> B[方法选择]
    B --> C1[代数构造: 立方差]
    B --> C2[递推假设: 多项式拟合]
    B --> C3[组合恒等式]
    B --> C4[几何模型]
    C1 --> D1[望远镜求和]
    C2 --> D2[解线性方程组]
    C3 --> D3[组合求和恒等式]
    D1 --> E[得出 S_n = n(n+1)(2n+1)/6]
    D2 --> E
    D3 --> E
    C4 --> D4[积分类比与修正]
    D4 --> E
```

8. 实际应用与算法意义

在算法复杂度分析中，嵌套循环的时间复杂度常涉及平方和：

for i in range(1, n+1):
    for j in range(1, i+1):
        # O(1) operation

内层执行次数为 $ 1 + 2 + \cdots + n = \frac{n(n+1)}{2} $，若再嵌一层，则可能涉及 $ \sum i^2 $。

掌握闭合形式可快速估算运行时间，避免仿真或数值逼近。