EGF

F^(x)=nanxnn!

(全文都是以 EGF 为基础)

封闭形式

n1xnn!=ex

  • 这个有关于麦克劳林级数(泰勒展开的一种特殊情况)

泰勒公式

  • xx0 处可导,那么当 xx0 时,函数的展开式近似于:

f(x)=f(x0)+f(x0)1!x+f(x0)2!x2...

麦克劳林级数是我们在 x00 时的展开,也就是泰勒公式的一个特殊情况

ln(1x)=n=11nxn

ln(1+x)=n=1(1)n+1nxn

排列和圆排列

  • 长度为 n 的排列数的指数生成函数是

P(x)=n0n!xnn!=11x

  • 圆排列的定义是排成一个环的方案

  • n 个数的圆排列方案为 (n1)!

Q(x)=n1(n1)!xnn!=n1xnn=ln(1x)=ln(11x)

  • 那么 exp(Q(x))=P(x)

exp 的直观理解

  • 上述中一个排列的每种方案其实就是将原排列分成若干集合,每种集合组成一种置换环

  • 那么方案数就是每种划分方案后的置换环的方案数的乘积

  • 这个置换环的方案数实际上就是圆排列的方案数

  • 形式化的,假设当前有 n 个元素,求 A 的方案数,设生成函数为 F(x)

  • 每种方案可以对应将 n 个元素先划分成若干个集合,每个集合 B 的方案数的乘积,设 B 的生成函数为 G(x)

  • 那么 exp(G(x))=F(x)

推广

  • 如果 n 个点 带标号 生成树的 EGF 是 G(x)n 个点 带标号 森林的 EGF 是 F(x) ,那么 F(x)=exp(G(x))

  • 如果 n 个点带标号无向连通图的 EGF 是 G(x) ,那么 n 个点带标号无向图的 EGF 是 exp(G(x))

应用

长度为 n 的一个错排指 pii

求错排的指数生成函数

  • 从置换环的角度来看,实际上就是不存在长度为 1 的置换环,生成函数 G(x)=n2xnn=ln(11x)x

  • 那么错排的生成函数是 exp(G(x))

求有多少个映射 f:{1,2,...,n}{1,2,...,n} ,使得

任意一个 i 映射 k 次后等于映射 k1 次后映射的数是同一个数

  • 考虑将 ifi 连边。相当于我们从任意一个 ik 步和走 k1 步到达的是同一个点

  • 那么整个图就是一些带一个自环的基环树组成的,且环的深度不超过 k ,假设根的深度为 1

  • 我们将环当成根,把边的方向反一下,那么整个图可以当成一些深度不超过 k 的生成树组成的

  • 假设深度不超过 k 的生成树的方案数的生成函数为 Fk(x)

  • 实际上深度不超过 k 的生成树就是用一个节点将深度不超过 k1 的森林连起来

  • 那么 Fk(x)=xexp(Fk1(x))

  • 答案是 [xn]exp(Fk(x))

给一个 n 个数的序列 a1,a2,...,an ,和初值为 0 的变量 s ,要求你重复一下操作 k

  • 1,2,...,n 中等概率选择一个 x
  • s 加上 Πixai
  • ax 减一

k 次操作后 s 的期望

1n5000,1k109,0ai109

  • 首先随便减两个数,发现最终的结果和减的顺序没有关系

  • 那么我们将一个数所有的减操作放在一起,再归纳一下,设每个数减少 bi,可以得出

s=i=1naii=1n(aibi)

  • 那么我们的问题就转换成了 i=1n(aibi) 的期望

  • 考虑计算每种方案的 i=1n(aibi) 的和,最后除 nk

  • k 次序列中,要使得 i 出现 bi 次的方案数是

k!b1!b2!...bn!

  • 这和指数生成函数的系数类似

  • aj 的指数生成函数是

Fj(x)=i0(aji)xii!

  • 那么答案就是 [xk]j=1nFj(x)

  • 为了快速计算答案,我们需要将 Fj(x) 转换成封闭形式

  • Fj(x)=ajexxex=(ajx)ex

  • j=1nFj(x)=enxj=1n(ajx)

  • j=1n(ajx) 是个 n 次多项式,可以暴力算,假设展开式是 i=0ncixi ,那么

j=1nFj(x)=(i0nixii!)(i=0ncixi)=i0j=0icjxjnijxij(ij)!=i0xii!j=0inijijcj

  • 答案就是这个多项式 xk 的系数
posted @   Kzos_017  阅读(81)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 地球OL攻略 —— 某应届生求职总结
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 提示词工程——AI应用必不可少的技术
· .NET周刊【3月第1期 2025-03-02】
点击右上角即可分享
微信分享提示