【学习笔记】自然数幂和

温馨提示: 本文文档大小约\(11KB\).

引入

自然数幂和是一个我们从小就耳熟能详的经典问题。定义\(S(n,k)=\sum^{n}_{i=0} i^k\), 显然\(S(n,k)\)为关于\(n\)的不超过\((k+1)\)次多项式，那么给定\(k\), 如何快速求出这个多项式的系数？或者给定\(n\)和\(k\), 如何快速求出\(S(n,k)\)? 这就是本文的讨论内容。

本文介绍三种不同的方法——拉格朗日插值法、第二类斯特林数法和伯努利数法，这三种方法在解决不同的问题上各有优劣。

符号约定

\(x^{\underline n}\)表示\(x\)的\(n\)次下降幂，即\(x^{\underline n}=\prod^{n}_{i=1}(x-i+1)\).

本文中提到的所有\(i\), 与虚数均无关系。

本文中所有运用数学归纳法进行证明的定理，归纳奠基均省略，请读者自行处理。

拉格朗日(Lagrange)插值法

概述

Lagrange插值法的应用范围很广，给定任意\((n+1)\)个互不相同的点值\((x_i,y_i)\) (\(i=0,1,...,n\))，其可以在\(O(n^2)\)时间内求出一个不超过\(n\)次的多项式\(f(x)\)，满足对于任意\(i\), \(f(x_i)=y_i\). (“互不相同”指对于任意\(0\le i<j\le n\)有\(x_i\ne x_j\))

暴力做法

设\(f(x_i)=\sum^{n}_{i=0} a_ix^i\). 根据\((n+1)\)个点值，我们可以列出\((n+1)\)个方程，第\(j\)个形如\(\sum^{n}_{i=0}x_j^ia_i=y_j\).

直接高斯消元求解，时间复杂度\(O(n^3)\).

插值多项式的唯一性

定理1 给定任意\((n+1)\)对互不相同的点值\((x_i,y_i)\)，恰好存在一个不超过\(n\)次的多项式\(f(x)\)，满足对于任意\(i\), \(f(x_i)=y_i\).

证明 根据暴力做法，我们就是要证明方程左边的\((n+1)\times (n+1)\)的系数矩阵可逆。

实际上系数矩阵就是著名的范德蒙德(Vandermonde)矩阵: \(\begin{bmatrix} 1&1&1&...&1\\ a_0&a_1&a_2&...&a_n\\ a_0^2&a_1^2&a_2^2&...&a_n^2\\ & & & ... & \\ a_0^n&a_1^n&a_2^n&...&a_n^n\end{bmatrix}\)

可以证明，该矩阵行列式