第一类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$

  1. (A,B)(C,D)
  2. (A,C)(B,D)
  3. (A,D)(B,C)
  4. (A)(B,C,D)
  5. (A)(B,D,C)
  6. (B)(A,C,D)
  7. (B)(A,D,C)
  8. (C)(A,B,D)
  9. (C)(A,D,B)
  10. (D)(A,B,C)
  11. (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$(自行前前面对比),

  1. (A,B)(C,D)
  2. (A,C)(B,D)
  3. (A,D)(B,C)
  4. (A)(B,C,D)
  5. (B)(A,C,D)
  6. (C)(A,B,D)
  7. (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数入门的更多相关文章

  1. Bell(hdu4767+矩阵+中国剩余定理+bell数+Stirling数+欧几里德)

    Bell Time Limit:3000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u Submit Status  ...

  2. lightOJ 1326 Race(第二类Stirling数)

    题目链接:http://lightoj.com/volume_showproblem.php?problem=1326 题意:有n匹马赛跑.问有多少种不同的排名结果.可以有多匹马的排名相同. 思路:排 ...

  3. 斯特灵数 (Stirling数)

    @维基百科 在组合数学,Stirling数可指两类数,都是由18世纪数学家James Stirling提出的. 第一类 s(4,2)=11 第一类Stirling数是有正负的,其绝对值是个元素的项目分 ...

  4. hdu 4372 第一类stirling数的应用/。。。好题

    /** 大意: 给定一系列楼房,都在一条水平线上,高度从1到n,从左侧看能看到f个, 从右侧看,能看到b个,问有多少种这样的序列.. 思路: 因为肯定能看到最高的,,那我们先假定最高的楼房位置确定,那 ...

  5. HDU 3625 Examining the Rooms:第一类stirling数

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3625 题意: 有n个房间,每个房间里放着一把钥匙,对应能开1到n号房间的门. 除了1号门,你可以踹开任 ...

  6. HDU 4372 Count the Buildings:第一类Stirling数

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4372 题意: 有n栋高楼横着排成一排,各自的高度为1到n的一个排列. 从左边看可以看到f栋楼,从右边看 ...

  7. 整理一点与排列组合有关的问题[组合数 Stirling数 Catalan数]

    都是数学题 思维最重要,什么什么数都没用,DP直接乱搞(雾.. 参考LH课件,以及资料:http://daybreakcx.is-programmer.com/posts/17315.html 做到有 ...

  8. [总结] 第二类Stirling数

    上一道例题 我们来介绍第二类Stirling数 定义 第二类Stirling数实际上是集合的一个拆分,表示将n个不同的元素拆分成m个集合的方案数,记为 或者 .和第一类Stirling数不同的是,集合 ...

  9. 贝尔数(来自维基百科)& Stirling数

    贝尔数   贝尔数以埃里克·坦普尔·贝尔(Eric Temple Bell)为名,是组合数学中的一组整数数列,开首是(OEIS的A000110数列):   Bell Number Bn是基数为n的集合 ...

随机推荐

  1. python 之 面向对象 (异常处理)

    7.15 异常处理 1.什么是异常 异常是错误发生的信号,程序一旦出错,如果程序中还没有相应的处理机制,那么该错误就会产生一个异常抛出来,程序的运行也随之终止 2.一个异常分为三部分: 异常的追踪信息 ...

  2. unity---为什么用Time.deltaTime * speed 表示每秒移动的距离的理解

    Time.deltaTime:代表时间增量,即从上一帧到当前帧消耗的时间, 这个值是动态变化的. dt 表示 deltaTime. 假如 1s渲染10帧,沿X轴方向的移动速度 speed = 10m/ ...

  3. arc079

    D. Decrease (Contestant ver.) 大意: 每次操作选一个最大数$-n$,其余数全$+1$. 要求构造一个序列$a$, 使得恰好$k$次操作后最大值不超过$n-1$. 只要让$ ...

  4. golang ---JSON-ITERATOR 使用

    jsoniter ( json-iterator )是一款快且灵活的 JSON 解析器 Jsoniter 是最快的 JSON 解析器.它最多能比普通的解析器快 10 倍之多, 独特的 iterator ...

  5. java之spring mvc之文件上传

    目录结构如下: 注意,下面说的配置文件,一般都是值的src下的配置文件,即mvc.xml.如果是web.xml,则直接说 web.xml 1. 文件上传的注意点 表单必须是post提交,必须将 enc ...

  6. python中字符串的常用(部分)处理方法

    myStr = "hello world itcast and hahaitcastcpp" 方法的查询方法: help(myStr.replace) myStr.find(&qu ...

  7. Grafana+prometheus+AlertManager+钉钉机器人

    一.Grafana (1)安装Grafana的Linux环境 在官网下载windows的Grafana的压缩包到指定目录,解压缩Grafana压缩文件到包含当前Grafana版本的文件夹.将该文件夹解 ...

  8. iptables限制访问

    iptables限制访问 常用命令 # 查看规则 iptables -L INPUT --line-numbers # 开放指定的端口 iptables -A INPUT -p tcp --dport ...

  9. 【开发工具】- Idea.2018.02注册码激活

    1.从下面地址下载一个jar包,名称是  JetbrainsCrack-3.1-release-enc.jar 下载地址: 链接: https://pan.baidu.com/s/1VZjklI3qh ...

  10. python day 18: thinking in UML与FTP作业重写

    目录 python day 18 1. thinking in UML读书小感 2. FTP作业重写 2.1 软件目录结构 2.2 FTPClient端脚本 2.3 FTPServer端脚本 pyth ...