载谭 Binomial Sum：多项式复合、插值与泰勒展开

对于微分有限（D-Finite） 的生成函数 $F(x)$，又设一生成函数 $G(x)$，有数列 $a$，我们已经对每个 $0\le k\le n$ 知晓了

$$\sum _{j=0}^n{a}_j[x^j]G(x{)}^k$$

那么我们可以在 $\Theta (n)$ 时间内计算

$$\sum _{j=0}^n{a}_j[x^j]F(G(x))$$

设 $c=G(0)$，考虑 $F$ 在 $c$ 处的泰勒展开式、带入 $G(x)$：

$$F(G(x))=\sum _{k\ge 0}{F}^{(k)}(c)\frac{(G(x)-c{)}^k}{k!}$$

由于 $x\mid G(x)-c$，可知对于我们关心的 ${\sum }_{j=0}^n{a}_j[x^j]F(G(x))$ 之计算而言，只展开到 $n$ 次项足够了。

设 $F(x)$ 满足的微分方程为 ${\sum }_{j=0}^m{p}_j(x)F^{(j)}(x)=0$，可知

$$\sum _{j=0}^m{p}_j(x+c)F^{(j)}(x+c)=0$$

我们考虑对 $F(x+c)$ 进行截取，设 $\mathcal{F}(x+c)=F(x+c)\mathrm{m}\mathrm{o}\mathrm{d}{x}^{n+1}$，考察原微分方程的影响

$$\sum _{j=0}^m{p}_j(x+c){\mathcal{F}}^{(j)}(x+c)$$

对于仅涉及 $\le n$ 的项来说，方程依旧成立。而对于仅涉及 $>n$ 的项来说也成立，因此我们知道在这个区间内，有一个扰动因子 ${\sum }_{j=0}^m{p}_j(x+c){\mathcal{F}}^{(j)}(x+c)=\mathcal{D}(x)$。因此有

$$\sum _{j=0}^m{p}_j(x){\mathcal{F}}^{(j)}(x)=\mathcal{D}(x-c)$$

依此进行递推。由于此时 $\mathcal{F}$ 只有不超过 $n$ 次项，可以直接依据已有的信息得到答案。

---

目前能见到的问题多为 $G(x)={\mathrm{e}}^x$ 的情况。这是线性插值的基本推广。

根据这一直接推论，我们统一了以下问题的 $\Theta (n)$ 结果。

1. 线性插值：$G(x{)}^k$
2. 前缀和：$\frac{1-G^{k+1}}{1-G}$
3. “Binomial Sum”：$(1+G{)}^k$
4. ${\sum }_{x\ge 0}{x}^n{q}^x$：$\frac{1}{1-qG}$
5. Bell 数：${\mathrm{e}}^{G-1}$
6. Bernoulli 数：$\frac{\ln (1+(G-1))}{G-1}$
7. “如何优雅地求和”：$(aG+(1-a){)}^k$
8. “TJOI 求和”：$\frac{1}{3-2G}$，这也是我当年第一次做出了一点有意思的结果，现在我们将它收尾。