Stirling数入门
第一类Stirling数
定义
$$
\begin{aligned}
(x)_n & =x(x-1)...(x-n+1)\\
&= s(n, 0) + s(n,1)x +..+s(n,n)x^n\\
\end{aligned}$$
例如,$n=3$ 时,
$(x)3 = x(x-1)(x-2)$
$(x)3 = x^0 + 2x -3x^2 + x^3$
于是 $s(3,0)=0,s(3,1)=2,s(3,2)=-3,s(3,3)=1$
有符号斯特林数和无符号斯特林数有如下关系:
$$s(n, k) = (-1)^{n-k}\begin{bmatrix} n\\ k \end{bmatrix}$$
下文的 $s(n, k)$ 都是指的无符号的了。
理解
$n$ 个人围着 $k$ 个相同圆桌,每个桌子非空的方案数就是 $s(n, k)$。
也就是将 $n$ 个不同元素分成 $k$ 组,每组中的元素再进行圆排列的方法数。
例如,$s(4, 2) = 11$
- (A,B)(C,D)
- (A,C)(B,D)
- (A,D)(B,C)
- (A)(B,C,D)
- (A)(B,D,C)
- (B)(A,C,D)
- (B)(A,D,C)
- (C)(A,B,D)
- (C)(A,D,B)
- (D)(A,B,C)
- (D)(A,C,B)
易得一个递推式,
人a独占一桌:$s(n-1, k-1)$
人a不独占一桌:先将 $n-1$ 个人安排好,再将a安排到任一人的右边,$(n-1)*s(n-1, k)$
所以,$s(n, k) = s(n-1. k-1) + (n-1)*s(n-1, k)$
第二类Stirling数
定义
与第一类Stirling数类似,可以用下降阶乘幂定义:
$$x^n = \sum_{k=0}^ns(n, k)(x)_k$$
例如,$n=3$ 时
$$x^3 = s(3, 0)(x)_0 + s(3,1)(x)_1 + s(3,2)(x)_2+s(3,3)(x)_3$$
即 $x^3 = s(3, 0) + s(3,1)x + s(3,2)x(x-1)+s(3,3)x(x-1)(x-2)$
开并合并同类项,得
$$x^3 = s(3,0) + [(3,1)-s(3,2)+2s(3,3)]x + [s(3,2)-3s(3,3)]x^2 + s(3,3)x^3$$
对比系数,解得
$s(3,0)=0, s(3,1)=1,s(3,2)=3,s(3,3)=1$
理解
将含有 $n$ 个元素的集合拆分成 $k$ 个非空子集的方法数就是第二类Stirling数。
也就是将 $n$ 个有区别的球放到 $k$ 个盒子里的方案数。
例如,$s(4,2)=7$(自行前前面对比),
- (A,B)(C,D)
- (A,C)(B,D)
- (A,D)(B,C)
- (A)(B,C,D)
- (B)(A,C,D)
- (C)(A,B,D)
- (D)(A,B,C)
也与第一类Stirling数有类似的递推式(初始条件都相同):
$$s(n, m) = s(n-1, m-1) + m*s(n-1, m)$$
证:
等价于将 $n$ 个有区别的球放到 $k$ 个盒子里的方案数,
若球a独占一盒,$s(n-1, m-1)$
若球a不独占一盒,先将剩下的 $n-1$ 个放入 $m$ 个盒子中且不允许有空盒,再将球a放入其中一盒,$ms(n-1, m)$.
补充:
参考链接:
1.https://blog.csdn.net/doyouseeman/article/details/50876786
2. https://zh.wikipedia.org/wiki/斯特林数
Stirling数入门的更多相关文章
- Bell(hdu4767+矩阵+中国剩余定理+bell数+Stirling数+欧几里德)
Bell Time Limit:3000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit Status ...
- lightOJ 1326 Race(第二类Stirling数)
题目链接:http://lightoj.com/volume_showproblem.php?problem=1326 题意:有n匹马赛跑.问有多少种不同的排名结果.可以有多匹马的排名相同. 思路:排 ...
- 斯特灵数 (Stirling数)
@维基百科 在组合数学,Stirling数可指两类数,都是由18世纪数学家James Stirling提出的. 第一类 s(4,2)=11 第一类Stirling数是有正负的,其绝对值是个元素的项目分 ...
- hdu 4372 第一类stirling数的应用/。。。好题
/** 大意: 给定一系列楼房,都在一条水平线上,高度从1到n,从左侧看能看到f个, 从右侧看,能看到b个,问有多少种这样的序列.. 思路: 因为肯定能看到最高的,,那我们先假定最高的楼房位置确定,那 ...
- HDU 3625 Examining the Rooms:第一类stirling数
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3625 题意: 有n个房间,每个房间里放着一把钥匙,对应能开1到n号房间的门. 除了1号门,你可以踹开任 ...
- HDU 4372 Count the Buildings:第一类Stirling数
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4372 题意: 有n栋高楼横着排成一排,各自的高度为1到n的一个排列. 从左边看可以看到f栋楼,从右边看 ...
- 整理一点与排列组合有关的问题[组合数 Stirling数 Catalan数]
都是数学题 思维最重要,什么什么数都没用,DP直接乱搞(雾.. 参考LH课件,以及资料:http://daybreakcx.is-programmer.com/posts/17315.html 做到有 ...
- [总结] 第二类Stirling数
上一道例题 我们来介绍第二类Stirling数 定义 第二类Stirling数实际上是集合的一个拆分,表示将n个不同的元素拆分成m个集合的方案数,记为 或者 .和第一类Stirling数不同的是,集合 ...
- 贝尔数(来自维基百科)& Stirling数
贝尔数 贝尔数以埃里克·坦普尔·贝尔(Eric Temple Bell)为名,是组合数学中的一组整数数列,开首是(OEIS的A000110数列): Bell Number Bn是基数为n的集合 ...
随机推荐
- 最新超简单解读torchvision
torchvision https://pytorch.org/docs/stable/torchvision/index.html#module-torchvision The torchvisio ...
- 前端向后端获取数据的三种方法:ajax、axios、fetch
1.jQuery中的ajax get方法: $.ajax({ url:"v4/api/film/now-playing?t=1539401039415&page=1&coun ...
- [转帖]浅谈分布式一致性与CAP/BASE/ACID理论
浅谈分布式一致性与CAP/BASE/ACID理论 https://www.cnblogs.com/zhang-qc/p/6783657.html ##转载请注明 CAP理论(98年秋提出,99年正式发 ...
- 微信JS从1.0.0升级到1.1.2的一个坑
因为1.0.0不支持电脑端日期Picker滚动,升级成了1.1.2,结果发现日期选择不起作用了经过跟踪发现 通过控制台查看 resut[1].toString()居然是number类型,修改代码为() ...
- Mysql 中 int(3) 和 int(11) 的区别
[1]int(3) 和 int(11)的区别(思维惯性认知错误) 这里的3或11代表的是存储在数据库中的具体的长度,总以为int(3)只能存储3个长度的数字,int(11)只会存储11个长度的数字. ...
- Spark学习(4) Spark Streaming
什么是Spark Streaming Spark Streaming类似于Apache Storm,用于流式数据的处理 Spark Streaming有高吞吐量和容错能力强等特点.Spark Stre ...
- 微信公众号开发 token 验证程序
<?php traceHttp(); define("TOKEN", "gmll001"); $wechatObj = new wechatCallbac ...
- 如何获取文件夹下所有文件名称(windows)
1. win+R -> cmd 打开dos命令窗口 2. 打开需获取文件名的位置 3. 获取名称 命令格式:dir /b >> 文件目标盘符:\文件夹位置(可省略)\目标名称.目标后 ...
- CLRS10.1-6练习 - 用双栈实现队列
双栈实现队列算法: 分别考虑队列两种操作入队和出队,我们假设使用栈s1 s2, s1用来模拟入队,s2用来模拟出队 入队: 入队操作直接执行s1.push即可 出队: 代码实现 package hel ...
- SQL Server 2012启动时提示:无效的许可证数据,需要重新安装
因为手咸,觉得电脑没有VS 2010版本的软件,就把Microsoft Visual C++ 2010某个组件给卸载了. 然后打开Sql Server 2012,就开始报错. 重装之后,也还是报错,将 ...