其他分享
首页 > 其他分享> > 浅谈 Binomial Sums 相关

浅谈 Binomial Sums 相关

作者:互联网

from EI

基本就是复述一遍,仅供参考。

\[\]

考虑如下问题:

对于一生成函数 \(G(x)\) 和一数列 \(a\),\(\forall 0\le k\le n\) 已知

\[\sum_{i=0}^na_i[x^i]G(x)^k \]

给出另一生成函数 \(F(x)\),求

\[\sum_{i=0}^na_i[x^i]F(G(x)) \]

若 \(F(x)\) 微分有限,并将微分方程相关视为常数,则我们可以做到 \(O(n)\) 的复杂度。

设 \([x^0]G(x)=c\),\(F(x)\) 满足的微分方程阶数为 \(r\)。

我们考虑求出一生成函数 \(H(x)\),满足 \(H(x+c)\equiv F(x+c) \pmod{x^{n+1}}\)。

那么有:

\[\begin{align*} &\sum_{i=0}^na_i[x^i]F(G(x))\\ =&\sum_{i=0}^na_i[x^i]F(G(x)-c+c)\\ =&\sum_{i=0}^na_i[x^i]\sum_{j\ge 0}([x^j]F(x+c))(G(x)-c)^j\\ =&\sum_{i=0}^na_i[x^i]\sum_{j=0}^n([x^j]H(x+c))(G(x)-c)^j\\ =&\sum_{i=0}^na_i[x^i]H(G(x))\\ =&\sum_{i=0}^nh_i\sum_{j=0}^na_j[x^j]G(x)^i \end{align*} \]

于是问题转为求 \(H(x)\)。

考虑 \(F(x)\) 满足的微分方程:

\[\sum_{i=0}^rp_i(x)F^{(i)}(x)=0 \]

则:

\[\sum_{i=0}^rp_i(x+c)F^{(i)}(x+c)=0 \]

考虑将 \(F\) 替换成 \(H\) 后的影响,此时我们只失去了 \(>n\) 的项,所以若某一项只与 \(\le n\) 的项有关或只与 \(>n\) 的项有关我们就不必去考虑它,唯一影响的只有在一段区间,我们设这个扰动因子为 \(D(x)\),那么:

\[\sum_{i=0}^ip_i(x)H^{(i)}(x)=D(x-c) \]

此时我们就可以递推出 \(H(x)\) 了。

现在大多数为 \(G(x)=e^x\) 的情况。

例:

参考资料

标签:le,frac,浅谈,na,sum,Sums,nF,Binomial
来源: https://www.cnblogs.com/leukocyte/p/15914997.html