Processing math: 100%
Live2D

Solution Set - Stirling 数相关杂题

  《好多题的题解》


  「洛谷 P5408」第一类斯特林数·行

  根据结论

x¯n=i[ni]xi,

我们只需要求出 x¯n 的各项系数。显然的 O(nlog2n) 做法就足够过掉洛谷上的原题了,但是我们 OJ 比较卓越,所以得用 O(nlogn) 的做法。

  考虑倍增。令 Fn(x)=x¯n,那么

F2n(x)Fn(x)=Fn(x+n)=ni=0[ni](x+n)i=ni=0[ni]ij=0(ij)xjnij=nj=0xjj!ni=j[ni]i!nij(ij)!.

可见 Fn(x+n) 是差卷积形式,O(nlogn) 求出后再和 Fn(x) 卷起来,O(nlogn) 完成倍增。总复杂度 T(n)=T(n/2)+O(nlogn)=O(nlogn)


   「洛谷 P5409」第一类斯特林数·列

  你永远可以相信组合意义。

  轮换数量的 EGF 为 G(x)=i>0(i1)!xii!,而显然“构成 m 个圆排列”等价于随便拿出 m 个总长为 n 的轮换,再为轮换中的数安排具体的值,最后除掉不同轮换间的顺序。令 R(x)=i[ik]xii!,那么

R(x)=Gk(x)k!.

O(nlogn) 快速幂即可。注意 [x0]G(x)1,需要位移后进行 ln 运算。


  「洛谷 P5395」第二类斯特林数·行

  考虑 {nm} 的组合意义:n 个不同的球放入 m 个相同的盒子且盒子不空的方案数。而 m!{nm}n 个不同的球放入 m 个不同的盒子且盒子不空的方案数。所以可以二项式反演求 {nm}

mn=mi=0(mi)i!{ni}{nm}=1m!mi=0(1)mi(mi)in=mi=0(1)mi(mi)!ini!.

好看又好记的式子。而这个通项恰好就是关于 m 这维的卷积形式,所以 O(nlogn) 卷一卷即可。


  「洛谷 P5396」第二类斯特林数·列

  组合意义嘛。把若干求放入一个盒子的 EGF 为 G(x)=i>0xii!=ex1,令 R(x)=i{ik}xii!,那么

R(x)=Gk(x)k!.

还是 O(nlogn) 快速幂。于是你惊讶地发现整个算法与第一类求列的区别只有:初始 G(x) 的系数从 x 的逆元变成了 x! 的逆元。所以这两类数确实是近亲(?)


  「HEOI 2016 / TJOI 2016」「洛谷 P4091」求和

f(n)=ni=0ij=0{ij}2jj!=nj=02jj!ni=0{ij}=nj=02jj!ni=0jk=0(1)jkkik!(jk)!=nj=02jj!jk=0(1)ni=0kik!(jk)!

特殊处理 k=0k=1,其余情况最内层和式可以等差数列求和,得到卷积形式。O(nlogn) 求出卷积再每项带 2jj! 的系数求和即可。


  「TC 13444」CountTables

  如果只有行上限制就可以直接组合数算,考虑用 DP 化简列上的组合问题。令 f(i) 表示 m=i 时的答案,总数减去不合法方案数,得到

f(i)=(cin)i1j=1{ij}f(j).

O(n2) 大力算就能过。


  「CF 715E」Complete the Permutations

  不管怎么说,这个 3400 还是虚高了。

  显然的结论:对于已补全的排列 pq,它们的距离为 n 减去置换 σ:piqi (i) 中的轮换数量。讨论一对 pi,qi 目前的状态:

  • a,b (ab0),将 a,b 视作同一点(缩点),环的数量不改变;
  • 0,b (b0),称之为 A 类边;
  • a,0 (a0),称之为 B 类边;
  • 0,0,称之为 C 类边。

  同时又注意到,若 ab,bc (b0) 已经确定,那么无论 a,c 是否为 0b 都可以被忽略,得到 ac。考虑完成上述缩点后该操作的执行:

  • A 类边 + A 类边,两者顺序任意,得到新的 A 类边;
  • A 类边 + B 类边,不可能执行;
  • A 类边 + C 类边,C 放 A 后面,得到新的 C 类边;
  • B 类边 + B 类边,两者顺序任意,得到新的 B 类边;
  • B 类边 + C 类边,C 妨 B 前面,得到新的 C 类边;

  因此,对于 A 类和 B 类边,要不自己连成完整的环,要不直接或间接地被缩入某个 C 类边。注意到以上操作不改变 C 类边的数量,且 ABC 间不相互影响。所以“A 类边构成环,剩下的缩入 C 类边”“B 类边构成环,剩下的缩入 C 类边”以及“C 类边构成环”三个组合问题是独立的。记 A,B,C 分别为三类边在完成所有缩边操作后的数量,考虑它们的方案数关于环数量的 GF:

[xk]GA(x)=Ai=0(Ai)[ik](Ai+C1)Ai_;[xk]GB(x)=Bi=0(Bi)[ik](Bi+C1)Bi_;[xk]GC(x)=C![Ck].

解释一下 GA(x):枚举 i 条自己成环,乘上选边和成环方案,其余 (Ai) 条边依次选一条边合并(可能是 A+A 或者 A+C,和剩下的任意边都能合并)。

  记 R(x)=GA(x)GB(x)GC(x),设初始时已经有 l 个环确定了,那么 [xnk]xlR(x) 即为距离为 k 的方案。中途涉及的运算全部暴力,O(n2) 可过。


  「HDU 4625」JZPTREE

  设某条路径的长度为 d,注意到

dk=k!ki=0{ki}(dk),

其中 (dk) 能够被赋予组合意义,对它换根 DP 即可。复杂度 O(Tnk)


  「FJOI 2016」「洛谷 P4609」建筑师

  非常自然的想法是将原问题用最大值的位置划分为两个关于最长贪心上升子序列的问题,但不要被这个想法限制了。令 AA1,BB1,先用除最大值外的 (n1) 个数构成 (A+B) 个轮换,不放认为轮换中的最大值在第一个的位置,那么把轮换按第一个位置的大小依次放置,就能得到对应长度的最长贪心上升子序列。而所有轮换确定后,最大值 n 自然只有一个位置能放(注意轮换内部形态)。因此,答案为

(A+BA)[n1A+B].

复杂度 O(nA)O(1)

posted @   Rainybunny  阅读(118)  评论(0编辑  收藏  举报
编辑推荐:
· 智能桌面机器人:用.NET IoT库控制舵机并多方法播放表情
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
阅读排行:
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· 新年开篇:在本地部署DeepSeek大模型实现联网增强的AI应用
· 程序员常用高效实用工具推荐,办公效率提升利器!
· Janus Pro:DeepSeek 开源革新,多模态 AI 的未来
· 【译】WinForms:分析一下(我用 Visual Basic 写的)
点击右上角即可分享
微信分享提示