Live2D

Note -「广义二项级数」浅赏

  上回 说到拉反和扩展拉反,那么这里先给几个小小变形或推广。

  基础 ver:

(1)[xn]G(x)=1n[xn1](F(x)/x)n,(2)[xn](HG)(x)=1n[xn1]H(x)(F(x)/x)n.

  令 H(x)=xk 带入 (2)

(3)[xn]Gk(x)=kn[xk]Fn(x).

  另一个变形,由于

kn[xk]Fn(x)=1n[xk]x(Fn(x))=[xk1]F(x)Fn1(x),

那么:

(4)[xn]Gk(x)=[xk1]F(x)Fn1(x).

既然都会求 Gk(x) 了,我们立马把 H(x)(4) 头上一套,得到:

(5)[xn](HG)(x)=[xn]H(x)F(x)(F(x)/x)n1(6)=[x0]xH(x)F(x)Fn1(x).

  做了一道例题之后不负责任的感言:(6) 挺有用的(事实上,后文只用了 (6) 这一种形式),或许是因为最多变的 H(x) 不用导了。


  关于 Bt(x),它是一个满足如下方程的 GF:

Bt(x)=xBtt(x)+1.

  当 t=2,这是 Catalan 数的 GF,于是乎依照 Catalan 的多样组合意义,你可以对 Bt(x) 进行包括但不限于多边形划分、走折线计数、 多叉树计数等组合意义解释。这个 GF 结合拉反,能够对相关 (tn+rn) 的求和进行有效变形。所以——

  学拉格朗日反演不学广义二项级数, 就像读四大名著不读红楼梦. 说明这个人文学造诣和自我修养不足, 他理解不了这种内在的阳春白雪的高雅艺术, 他只能看到外表的辞藻堆砌, 参不透其中深奥的精神内核, 他整个人的层次就卡在这里了, 只能度过一个相对失败的人生.

  整活完毕,我们来推一下 Bt(x) 的系数。令 F(x)=Bt(x)1 以去掉常数项,观察发现 F 的复合逆为 G(x)=x(1+x)t,那么

[xn]Btr(x)=[xn](1+F(x))r=[x0]x(1+x)r(1+x)txt(1+x)t1(1+x)2t(1+x)tn+txn+1=[x0]x(1+x)r1(t1)x(1+x)t+1(1+x)tn+txn+1=[xn](1+x)tn+r1(1(t1)x)=(tn+r1n)(t1)(tn+r1n1)=rtn+r(tn+rn).

其中第二步用公式 (6),代入 H(x)=(1+x)r 反演。

  然后你就会说了,它这组合数带个 rtn+r 的系数,格局就上不去,根本做不了题。我们来尝试弄一下 n0(tn+rn)xn 的表达式

(tn+rn)=[xn](1+x)tn+r=[x0]((1+x)t/x)n(1+x)r=[x0]Gn(x)(1+F(G(x)))r=[x0]xxn(1+F(x))rF(x)/F(x)=[xn](1+F(x))rxF(x)F(x)=[xn](1+F(x))rx(G(F(x)))1F(x)=[xn](1+F(x))rG(F(x))G(F(x))F(x)=[xn](1+F(x))r1+F(x)1(t1)F(x).

  神他妈的,第三步时,令 H(x)=xn(1+F(x))r 代公式 (6),居然能这么反演,我兔麻了。最后一步展开 G(x)G(x) 整理即可,不多讲。


  然后是喜闻乐见的例题时间。

 

  初步转化,枚举 x,求合法方案数。当 x 确定时,方案数为

ansx=b=0Br=bxR(Bb)x(r+br)(Rr+BbRr)=b=0Bd=0RBx(b(x+1)+db)((Bb)(x+1)+RBxdBb)=d=0RBxb=0B(b(x+1)+db)((Bb)(x+1)+RBxdBb).

t=x+1,你发现这个组合数 (bt+db) 我们好像挺会的。令 F(z)=Bt(z)1,带入结论:

ansx==d=0RBxb=0B([zb](1+F(z))d1+F(z)1(t1)F(z))([zBb](1+F(z))RBxd1+F(z)1(t1)F(z))=[zB]d=0RBx(1+F(z))RBx(1+F(z)1(t1)F(z))2=(RBx+1)[zB](1+F(z))RBx(1+F(z)1(t1)F(z))2.

然后我们叕来一次拉反,令 H(z)=(1+z)RBx(1+z1(t1)z)2,反正已经有 FG=z,丢进去反演:

ansxRBx+1==[z0]zH(z)G(z)GB1(z)=[z0]z(1+z)RBx(1+z1(t1)z)21(t1)z(1+z)t+1(1+z)Bt+Bzn+1=[zB](1+z)R+1t+B1(t1)z=[zB](1+z)R+B+11xz=i=0B(R+B+1i)xBi.

  最终,答案就是

ans=x=1R/B(RBx+1)i=0B(R+B+1i)xBi=i=0B(R+B+1i)((R+1)x=1R/BxBiBx=1R/BxBi+1).

于是乎求一堆范围固定,上指标可变的自然数幂和即可。

  为了水篇幅我们来推一下这题自然数幂和的求法。令 fk=i=0nik,并令其 EGF 为 F(z),那么

F(z)=i0zii!j=0nji=j=0ni0jizii!=j=0nejz=1e(n+1)z1ez.

分子分母约一个 z,多项式求逆求出 F(z)modzB+2 即可。复杂度 O(BlogB)

posted @   Rainybunny  阅读(328)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现
点击右上角即可分享
微信分享提示