好久没有写过数学方面的博客了。
方程(1)
\[P(x)=\sum_{k=0}^{n}a_kx^k=0,a_n=1\]对应递推式
\[\sum_{k=0}^{n}a_kF_{k}=0\]即
\[F_n=-\sum_{k=0}^{n-1}a_kF{_{k}}\]令(1)的解为$x_1,x_2\ldots x_n$。
定义$S_m:=\sum_{k=1}^{n}x_k^m$,$m\le n$。显然当$m\gt n$时$S$符合$F$的递推式。
将\(S_m\)带入\(F\)的递推式中,不能完美符合。观察误差项的规律,得到
\[S_m=a_{n-m}(n-m)-\sum_{k=n-m}^{n-1}a_kS_{k-n+m}\] \[=a_{n-m}(n-m)-\sum_{k=1}^{m}a_{n-k}S_{m-k}\](这里坐标变换把我绕晕了很久)
这个式子几乎是一个递推式,相当于真正的\(S_0=n\),然而递推中把\(S_0\)看成了\(m\)
移一下项变成
\[a_{n-m}(n-m)=\sum_{k=0}^{m}a_{n-k}S_{m-k}\]令\(a'_k=a_{n-k}\)
\[a'_m(n-m)=\sum_{k=0}^{m}a'_{k}S_{m-k}\]多项式求逆板子。
一个多项式的逆序和另一个多项式逆序的逆做卷积,众所周知这是多项式除法。只不过保留的项数不只\(n_1-n_2\)。
写成相当简洁的封闭形式是
\[S=xP'(x)\div P(x)\]这背后到底是什么奥秘?
然后今天上Google搜了一下,搜索关键字是本文的标题Sum of Powers of Roots
。搜索中文内容是永远都啥都搜不到的(无论是谷歌还是百度)。
于是看到了这篇文章
惊讶的发现这是我已经学过的东西,我居然忘了。上次看zzq的博客学到了任意一个x数组的k次幂求和,然而那篇的重点是初等对称多项式,所以之前没有联想到。
那篇文章里面有一个简单的证明。
\[P(x)=\prod_{k=1}^{n}(x-x_k)\]两侧取对数得到
\[\log P(x)=\sum_{k=1}^{n}\log(x-x_k)\]求导得到
\[\begin{aligned}\frac{P'(x)}{P(x)}& =\sum_{k=1}^{n}\frac{1}{x-x_k} \\ & =\frac{1}{x}\sum_{k=1}^{n}\frac{1}{1-x_k/x} \\ & =\frac{1}{x}\sum_{k=1}^{n}\sum_{i=0}^{\infty}(\frac{1}{x})^ix_k^i \\ & =\frac{1}{x}\sum_{i=0}^{\infty}S_ix^{-i}\end{aligned}\]两侧逆序一下再乘个x就得到了之前找规律的式子
这个憨逼洛谷不支持\begin{aligned}
,然而我自己的博客支持。