ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

浅谈 Binomial Sums 相关

2022-02-20 11:00:37  阅读:220  来源: 互联网

标签:le frac 浅谈 na sum Sums nF Binomial


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\) 的情况。

例:

  • \(\sum_{i=0}^ma_i[\frac{x^i}{i!}]e^{nx}\)

    \(F(x)=x^n\)

    显然有:

    \[F(x)'x=nF(x) \]

    考察 \(x^m\) 次项,有:

    \[(m+1)\binom{n}{m+1}(x-1)^m+H(x)'x=nH(x) \]

  • \(\sum_{i=0}^ma_i[\frac{x^m}{k!}](1+e^x)^n\)

    \(F(x)=(1+x)^n\)。

    显然有:

    \[F(x)'(1+x)=nF(x)\\ \]

    考察 \(x^m\) 项,可以得到:

    \[2(m+1)\binom{n}{m+1}2^{n-m-1}(x-1)^k+H(x)'(1+x)=nH(x) \]

  • \(\sum_{i=0}^m a_i[\frac{x^i}{i!}](qe^x+(1-q))^n\)

    \(F(x)=(qx+(1-q))^n\)。

    显然有:

    \[F(x)'(x+\frac{1-q}{q})=nF(x) \]

    考察 \(x^m\) 次项,可以得到:

    \[q^m(m+1)\binom{n}{m+1}(x-1)^m+H(x)'(x+\frac{1-q}{q})=nH(x) \]

参考资料

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

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有