系数为多项式的幂级数求和

最新推荐文章于 2021-04-29 15:00:27 发布

MILLOPE

最新推荐文章于 2021-04-29 15:00:27 发布

阅读量1.2k

点赞数 2

文章标签： ldy

求
$\sum_{n=0}^{ \infty }f(n)r^{n}$
其中 $\ r \in (0,1)$ ， $\ f(n)$ 是一个最高次为 $\ k$ 的多项式。

我们先求最简单的形式：

$s=\sum_{n=0}^{ \infty }r^{n}$

即

$s=1+r+r^{2}+r^{3}+\cdots$

$\cdot s=r+r^{2}+r^{3}+\cdots$

两式相减得

$\cdot s=1$

$\frac{1}{1-r}$

在考虑其他简单形式

$s=\sum_{n=0}^{ \infty }n \cdot r^{n}$

$s=r+2r^{2}+3r^{3}+4r^{4}+ \cdots$

$\cdot s=r^{2}+2r^{3}+3r^{4}+4r^{5}+ \cdots$

相减得

$\cdot s=r+r^{2}+r^{3}+r^{4}+ \cdots$

这样和第一种情况一样，答案为

$s=\frac{r}{(1-r)^{2}}$

这使我们意识到，一个高次多项式的问题可以转换为低次多项式的问题。
为简化算法，我们需要多项式转化为一坨相邻数字连乘。
由于公式（提醒一句， $\ k$ 是最高次）

$=a_{0}\binom{n}{1}+a_{1}\binom{n}{2}+ \cdots +a_{k}\binom{n}{k}\,\,\,\,\,\,\,\,\,(证明（1）)$

所以我们就一个组合数考虑

$s=\sum_{n=0}^{ \infty }\binom{n}{m}r^{n}$

$\binom{n}{m}=\frac{n(n-1)(n-2) \cdots (n-m+1)}{m!}$

提出 $\ m!$ 得

$\frac{1}{m!} \sum_{n=0}^{ \infty }n(n-1)(n-2) \cdots (n-m+1)r^{n}$

$\cdot s$

$n<m时，\binom{n}{m}=0$

$\times 2 \times \cdots m \cdot r^{m}+2 \times 3 \times \cdots (m+1) \cdot r^{m+1}+3 \times 4 \times \cdots (m+2)r^{m+2}+ \cdots$

$\cdot s0 =1 \times 2 \times \cdots m \cdot r^{m+1}+2 \times 3 \times \cdots (m+1) \cdot r^{m+2}+3 \times 4 \times \cdots (m+2)r^{m+3}+ \cdots$

相减得

$\times 2 \times \cdots m \cdot r^{m}+2 \times 3 \times \cdots m \cdot m \cdot r^{m+1}+3 \times 4 \times \cdots (m+1) \cdot m \cdot r^{m+2}+ \cdots$

可以看到我们将除了 $r^{m}$ 以外的系数多了m。我们继续用这样方法

$设 s 1 = (1 - r) s 0$

$\times 2 \times \cdots m \cdot r^{m}+2 \times 3 \times \cdots m \cdot m \cdot r^{m+1}+3 \times 4 \times \cdots (m+1) \cdot m \cdot r^{m+2}+ \cdots$

$\cdot s1=1 \times 2 \times \cdots m \cdot r^{m+1}+2 \times 3 \times \cdots m \cdot m \cdot r^{m+2}+3 \times 4 \times \cdots (m+1) \cdot m \cdot r^{m+3}+ \cdots$

两式相减得

$\times 2 \times \cdots m \cdot r^{m}+2 \times 3 \times \cdots m(m-1)r^{m+1}+3 \times 4 \times \cdots m \cdot m(m-1) r^{m+2} +\cdots$

$设 s 2 = (1 - r) s 1$

$\times 2 \times \cdots m \cdot r^{m}+2 \times 3 \times \cdots m(m-1)r^{m+1}+3 \times 4 \times \cdots m \cdot m(m-1) r^{m+2} +\cdots$

$\cdot s2=1 \times 2 \times \cdots m \cdot r^{m+1}+2 \times 3 \times \cdots m(m-1)r^{m+2}+3 \times 4 \times \cdots m \cdot m(m-1) r^{m+3} +\cdots$

两式相减得

$\times 2 \times \times \cdots m \cdot r^{m}+2 \times 3 \times \cdots m(m-2)r^{m+1}+3 \times 4 \times \cdots m(m-1)(m-2)r^{m+2}+ \cdots$

以此类推，最后得出（提醒一句， $\ k$ 是最高次）

$\frac{r^{m}}{(1-r)^{m+1}}$

$s=\frac{s0}{m!}=\frac{r^{m}}{(1-r)^{m+1}}$

这样我们就得到了一个通项公式
现在我们求的就是 $\ a_{0},a_{1},a_{2} \cdots a_{k}$
由于

$f(n)=a_{0}\binom{n}{0}+a_{1}\binom{n}{1}+ \cdots +a_{k}\binom{n}{k}$

我们可以证明出

$\sum_{n=0}^{ \infty }f(n)r^{n}=\sum_{n=0}^{k}a_{n} \frac{r^{n}}{(1-r)^{n+1}}$

求得 $\ a_{0},a_{1}, \cdots a_{k}$ 即可求得正解（提醒一句， $\ k$ 是最高次）
又由于

$f(0)=a_{0}\binom{0}{0}$

$f(1)=a_{0}\binom{1}{0}+a_{1}\binom{1}{1}$

$f(2)=a_{0}\binom{2}{0}+a_{1}\binom{2}{1}+a_{2}\binom{2}{2}$

$\cdots \, \cdots \, \cdots \, \cdots$

$f(k)=a_{0}\binom{k}{0}+a_{1}\binom{k}{1}+ \cdots +a_{k}\binom{k}{k}$

这样就变成了一个多元方程组的问题，整理得

$a_{0}=\binom{0}{0} f{0}$

$a_{1}=-\binom{1}{0} f(0)+\binom{1}{1} f(1)$

$a_{2}=\binom{2}{0} f(0)-\binom{2}{1} f(1)+\binom{2}{2} f(2)$

$\cdots \, \cdots \, \cdots \, \cdots$

$a_{k}=(-1)^{k}\binom{k}{0} f(0)+(-1)^{k-1}\binom{k}{1} f(1)+ \cdots (-1)^{0} \binom{k}{k} f(k)$

对于计算机而言，这些问题都是可以在 O(n*n) 之内求得。优化后，可以使用多项式多点求值和NTT求得。O(n logn logn)

证明（1）

$f(n)=a_{0}\binom{n}{0}+a_{1}\binom{n}{1}+ \cdots +a_{k}\binom{n}{k}$

有于每个组合数展开后形成一个多项式，最高次递增，其各项系数可以满足多元方程组。

$由于\binom{n}{m}是一个关于n的m次多项式$

设

$\binom{n}{m}=b_{m,0}+b_{m,1}n+b_{m,2}n^{2}+ \cdots b_{m,m}n^{m}$

我们得到方程组

$\binom{n}{0}=b_{0,0}$

$\binom{n}{1}=b_{1,0}+b_{1,1}n$

$\binom{n}{2}=b_{2,0}+b_{2,1}n+b_{2,2}n^{2}$

$\cdots \, \cdots \, \cdots \, \cdots$

易证 $\ b_{ij} = \neq 0 \,\,\,\,\,\,\,\,\,\,\,\, i \in [0,m],j \in [0,m]$

设原方程

$f(n)=c_{0}+c_{1}n^{1}+c_{2}n^{2}+ \cdots c_{k}n^{k}$

有方程组

$\underbrace{ a_{0}b_{0,0}+a_{1}b_{1,0}+a_{2}b_{2,0}+ \cdots +a_{k}b_{k,0} }_{k+1项}=c_{0}$

$\underbrace{ a_{1}b_{0,1}+a_{2}b_{1,1}+a_{3}b_{2,1}+ \cdots +a_{k}b_{k,1} }_{k项}=c_{1}$

$\underbrace{ a_{2}b_{0,2}+a_{3}b_{1,2}+a_{4}b_{2,2}+ \cdots +a_{k}b_{k,2} }_{k-1项}=c_{2}$

$\cdots \, \cdots \, \cdots \, \cdots$

$a_{k}b_{k,k}==c_{k}$

由于项数不断减少，此方程显然有解

MILLOPE

关注

2
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
系数为多项式的幂级数求和

求∑n=0∞f(n)rn\sum_{n=0}^{ \infty }f(n)r^{n}n=0∑∞f(n)rn其中 r∈(0,1)\ r \in (0,1) r∈(0,1)， f(n)\ f(n) f(n)是一个最高次为 k\ k k的多项式。我们先求最简单的形式：s=∑n=0∞rns=\sum_{n=0}^{ \infty ...
复制链接

扫一扫