斯特林数和欧拉数(还没学完 但是学不动了于是开摆)

斯特林数和欧拉数(还没学完 但是学不动了于是开摆)

斯特林数

生成函数

同一列第二类斯特林数:一个盒子的 egf 是 ex1,有 k 个盒子但是盒子之间无区分,所以同一列斯特林数的 egf 是 (ex1)k/k!.还有一个东西是根据它的递推式解出来 ogf 是 xk(1ix)

如果不限定盒子个数,那么 egf 是 exp(ex1),这就是贝尔数的 egf。

同一行第二类斯特林数:利用 egf 列通项:{nk}=n!k![xn]i=0k(ki)ei(1)ki=i=0kin(1)kii!(ki)!,对所有 k 求这个式子的值于是卷积计算。

同一列第一类斯特林数:单个轮换 egf 是置换 egf 取 ln,而置换 egf 是 11x 那么轮换的 egf 是 ln(11x)=ln(1x),和上面同理计算其 k 次方除掉 k! 就行。

同一行第一类斯特林数:利用递推式解出来 ogf 是 xn¯,直接分治 ntt 是 2log。

套路是给定 F(x) 的系数,可以通过一次卷积求出 F(x+c) 的系数,列出式子之后交换求和号,发现是个差卷积的形式。利用这个将所有 x2i¯ 算出之后选取一些乘起来,复杂度是 T(n)=T(n/2)+O(nlogn)=O(nlogn)

幂之间的转换

xn=k{nk}xk_ xn¯=k[nk]xk

这个应该怎么记忆,xn¯xnxn_,所以应该是若干个下降幂乘上一个什么系数得到普通幂,普通幂得到上升幂。利用 xn_=(1)n(x)n¯ 再还元就能得到:

xn=k{nk}(1)nkxk¯xn_=k[nk](1)nkxk

反转公式

上一节中把其中一个往另一个式子的右侧代入,能够得到四个比较对称的式子,两边提取 [xn] 就得到了 [n=m] 的展开,为了更好记忆这里记作 [l=r],那么就有:

k=lr{rk}[kl](1)kl=[l=r]k=lr{rk}[kl](1)rk=[l=r]k=lr[rk]{kl}(1)kl=[l=r]k=lr[rk]{kl}(1)rk=[l=r]

这四个式子非常对称,仅需要记忆 {rk}[kl][rk]{kl} 可以互换;(1)kl(1)rk 可以互换就能把这些记住。

斯特林反演

f(n)=i=0n{ni}g(i)g(n)=i=0n(1)ni[ni]f(i)

证明(来自 cjl)

尝试利用 egf 解决问题,令 F,Gf,g 的 egf:

F(x)=kfkxkk!=k=0xkk!{ki}gi=igik=i{ki}xkk!=i=0gi(ex1)ii!=G(ex1)

ex1 的复合逆是 ln(1+x),于是有 F(ln(1+x))=G(x),那么对于另一方向:

G(x)=kfkln(1+x)kk!=kfk(1)k(ln(1(x)))kk!=kfk(1)ki=k[ik](x)ii!=ixii!k=i[ik](1)ikfk

类比二项式反演,尝试证另一个方向的斯特林反演:

f(m)=km{km}g(k)g(m)=km(1)km[km]f(k)

证明(感谢 qwaszx 指点)

来类比这里的证明,普通幂和斯特林数不是很搭,用下降幂来解决。

m0f(m)zm_=m0km{km}g(k)zm_=k0g(k)mk{km}zm_=k0g(k)zk

直接两边提取 [zk],那么左侧要计算 f(m) 的贡献系数是 [zk]zm_=[mk](1)mk(下降幂转普通幂),于是就完成了证明。每一步都是当且仅当所以逆命题也成立。

恒等式

具体数学 P221。

(6.15) 组合意义 (nk) 是把 n+1 所在的集合给选出来;(6.16) 归纳不是很难(这个怎么组合意义啊?);(6.17) 和 (6.18) 是 (6.15) 和 (6.16) 二项式反演得到的。

upd:感谢鱼鱼教我 (6.16) 的简洁证明:

k=0n[nk](km)=[xm]k=0n{nk}(1+x)k=[xm](1+x)n¯=[xm]1xxn+1¯=[n+1m+1]

(6.19) 就是同一行第二类斯特林数那个式子,另一个推法是对 mn=k(mk){nk}k! 二项式反演,这个式子是 n 个不同球放进 m 个可空盒子方案数,右侧是枚举非空盒子个数。

(6.20) 组合意义,考虑哪个盒子中最小编号球最大,枚举这个球是哪个,对于前面的直接 {km} 分出来 m 个盒子,对于后面 nk 个球每个球都有 m+1 个盒子可以选择放进去。(6.21) 也一样,后面 nk 个球放入的方案数从左到右分别是 k+1,k+2n,乘起来就是 nnk_

(6.22) (6.23):回想对二项式系数斜线求和是怎么处理的,那么这里也一样,将 {n+00} 看成 {n+10} 然后一层层合并,(6.23) 也一样。

/fn 后面的不想看了开摆

互不相等容斥

老题新做一下 TopCoder 13444 CountTables。广为人知的做法是 fi 表示 n×i 矩阵每行互不相同的个数,gi每行每列互不相同的个数。将每一列看成一个元素,那么将 fi 中的相等元素合并起来就得到了一个 gj,所以 fi=ji{ij}gj,而 fi=(ci)n_ 从而可以利用斯特林反演计算出 gm 的值。

还是想用常规的容斥思路来算。现在问题是有 n 个元素 a1,a2,an,每个元素都有若干取值可以取,计算满足 a 两两不同的方案数。那么限制就是 i<j[aiaj],写成 i<j(1[ai=aj]) 后将其展开,每个括号取 1 或者 [ai=aj],那么相等关系就会连成若干个连通块,一条边连起来表示这个括号取了 [ai=aj],没连表示取 1

把每个连通块看成一个点的集合,现在对于 1n 的一个集合划分计算容斥系数和,其等于各个连通块容斥系数和的乘积。那么现在要计算的就是对于这个连通块内部的边,有多少边集 E 能够将所有点连通,并且带来 (1)|E| 容斥系数的贡献。

用总的贡献减去不合法的贡献。假设现在有 k 个点,取边 (1,2),对于任意一个方案而言将其存在状态取反即得到另一个边个数奇偶性不同的方案,这是一个双射,两个方案权值相反可以抵消。仅有 k=1 的时候选不出一条边,所以总的贡献和是 [k=1]

假设 fi 表示 i 个点的答案,计算不合法贡献就枚举 i 所在连通块是哪些点,那么就有 fi=[i=1]1ji1(i1j1)fj[ij=1],所以 i>1 时就有 fi=(1)(i1)fi1,那么 fk=(1)k1(k1)!

假如枚举了一个集合划分 {S1,S2,,Sk},那么所有能得到这个集合划分的钦定方案的容斥系数的和就是 (1)|Si|1(|Si|1)!

注意到 (|Si|1)! 就是环排列方案数,而 (1)|S|1 看成初始有 (1)n 每多一个 S 就乘一个 (1),所以它其实和斯特林反演是等价的。

欧拉数

欧拉数:恰好满足有 k 个升高的排列 π1π2πn 的个数。其中升高定义为 πi<πi+1i 的个数,记做 nk,考虑将 n(n1) 的排列中插入可得递推式:nk=(k+1)n1k+(nk)n1k1

回想起不等关系,尝试用容斥去算欧拉数 nm。假定钦定了若干个间隔是 <,它们一共形成了 k 个极长连续上升段,那么一共钦定了 nk 个升高,根据二项式反演其容斥系数为 (1)nkm(nkm),而考虑每个极长连续段看成一个盒子,那么方案数就是将 n 个球放进 k 个盒子中,并且盒子之间有顺序(由于每个段必须要求上升所以不需要考虑球之间的顺序),那么方案数是 {nk}k!,至此得到了 (6.40):

(6.40)nm=k{nk}(nkm)(1)nkmk!

而 (6.39) 即为 (6.40) 逆用二项式反演得到的结果,在 (6.40) 中令 knk 换元即可看出 f(i)=ni,g(i)={nni}(ni)!f(m)=km(km)(1)kmg(k) 从而导出 g(m)=km(km)f(k),再换元 mnm 就得到了 (6.39).

posted @   do_while_true  阅读(157)  评论(2编辑  收藏  举报
相关博文:
阅读排行:
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· Manus的开源复刻OpenManus初探
· 三行代码完成国际化适配,妙~啊~
· .NET Core 中如何实现缓存的预热?
· 阿里巴巴 QwQ-32B真的超越了 DeepSeek R-1吗?

This blog has running: 1845 days 1 hours 33 minutes 19 seconds

点击右上角即可分享
微信分享提示