Stirling数笔记
Updating....
这几个玩意儿要记的东西太多太乱所以写blog整理一下
虽然蒯的成分会比较多全部
我居然开始记得写blog了??
第一类
这里讨论的是无符号类型的。
OEIS编号A130534
表示方法
\(s(n,m)\)或者\(\begin{bmatrix}n \\ m\end{bmatrix}\)
注意前者是小写s
意义
\(n\)个元素的项目分作\(m\)个非空环排列的方法数目
求法
递归求解法
\[\begin{bmatrix}n\\m\end{bmatrix}=\begin{bmatrix}n-1\\m-1\end{bmatrix}+(n-1)\begin{bmatrix}n-1\\m\end{bmatrix}\]
这个就是说新建一个环排列或者插入已有的环排列
可怕 这很\(O(n^2)\)
各种性质
\(\begin{bmatrix}n\\1\end{bmatrix}=(n-1)!\)
\(\begin{bmatrix}n\\2\end{bmatrix}=(n-1)!\times\sum_{i=1}^{n-1}\frac 1 i\)
\(\sum_{i=0}^n \begin{bmatrix}n\\i\end{bmatrix}=n!\)
\(\begin{bmatrix}n\\n-1\end{bmatrix}=\binom{n}{2}\)
这里就不给出证明了
别的地方都有
也挺好记好想的
maybe
第二类
OEIS编号A008277
表示方法
\(S(n,m)\)或者\(\left\{\begin{matrix}n \\ m\end{matrix}\right\}\)
当然这里是大写S
意义
\(n\)个元素的集定义\(m\)个等价类的方法数目
。。。wiki害人
就是从环排列变成集合划分了
当然也要保证非空
求法
递归求解法
\[\begin{Bmatrix}n\\m\end{Bmatrix}=\begin{Bmatrix}n-1\\m-1\end{Bmatrix}+m\begin{Bmatrix}n-1\\m\end{Bmatrix}\]
同样也可以解释,新建or插入已有的
再次\(O(n^2)\)??别啊
幸好这玩意儿能搞容斥,通项就有了
\[\begin{Bmatrix}n\\m\end{Bmatrix}=\frac{1}{m!}\sum\limits_{k=0}^{m}(-1)^k\binom{m}{k}(m-k)^n\]
\(O(n)\)求解不是梦
好吧只求一个用这个会快
最重要的是这个能卷,也好搞些别的???
稍微整理一下
\[\left\{\begin{matrix}n\\m\end{matrix}\right\}=\sum_{k=0}^m\frac{(-1)^k}{k!}\frac{(m-k)^n}{(m-k)!}\]
就很舒服
Stirling 反演
两个柿子挺好记
但我暂时还搞不清具体是干嘛的。。。
\[f(x) = \sum_{i=0}^x \begin{Bmatrix}x\\i\end{Bmatrix} g(i) \Leftrightarrow g(x) = \sum_{i=0}^x (-1) ^ {x - i}\begin{bmatrix}x\\i\end{bmatrix} f(i)\]
\[f(x) = \sum_{i=0}^x \begin{bmatrix}x\\i\end{bmatrix} g(i) \Leftrightarrow g(x) = \sum_{i=0}^x (-1) ^ {x - i}\begin{Bmatrix}x\\i\end{Bmatrix} f(i)\]
Bell数
OEIS编号A000110
就是把第二类stirling数的集合划分个数限制去掉了
只限制了基数
也就是
\[B_n=\sum_{i=0}^n\begin{Bmatrix}n\\i\end{Bmatrix}\]
当然也可以直接\(O(n)\)递推
\[B_{n+1}=\sum_{k=0}^n\binom{n}{k}B_{k}\]
参考
%%%
https://www.cnblogs.com/NaVi-Awson/p/9242645.html
https://www.cnblogs.com/ezoiLZH/p/9424911.html
https://www.cnblogs.com/owenyu/p/6724661.html
https://blog.csdn.net/winycg/article/details/70233717
Stirling数笔记的更多相关文章
- 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数不同的是,集合 ...
- Bell(hdu4767+矩阵+中国剩余定理+bell数+Stirling数+欧几里德)
Bell Time Limit:3000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit Status ...
- 贝尔数(来自维基百科)& Stirling数
贝尔数 贝尔数以埃里克·坦普尔·贝尔(Eric Temple Bell)为名,是组合数学中的一组整数数列,开首是(OEIS的A000110数列): Bell Number Bn是基数为n的集合 ...
随机推荐
- HDU 3467 Song of the Siren(圆交)
Problem Description In the unimaginable popular DotA game, a hero Naga Siren, also known as Slithice ...
- POJ 1815 Friendship(最大流最小割の字典序割点集)
Description In modern society, each person has his own friends. Since all the people are very busy, ...
- H5页面 绝对定位元素被 软键盘弹出时顶起
H5页面 绝对定位元素被 软键盘弹出时顶起 在h5页面开发的过程中,我们可能会遇到下面这个问题,当页面中有输入框的时候,系统自带的软盘会把按钮挤出原来的位置.那么我们该怎么解决呢?下面列出一下的方法: ...
- “Hello World!”团队第三周召开的第一次会议
今天是我们团队“Hello World!”团队第三周召开的第一次会议.博客内容: 一.会议时间 二.会议地点 三.会议成员 四.会议内容 五.Todo List 六.会议照片 七.燃尽图 一.会议时间 ...
- 本周PSP图
本周共写博文5篇,共计4800字,知识点:知道了博客应当如何写,接触了博客园,阅读了构建之法 内容 开始时间 结束时间 中断时间 共计时间 9月8日博文 22:00 22:55 10min聊天 45m ...
- Swift-闭包理解
/* 闭包(Closures) * 闭包是自包含的功能代码块,可以在代码中使用或者用来作为参数传值. * 在Swift中的闭包与C.OC中的blocks和其它编程语言(如Python)中的lambda ...
- js图片转换为base64
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- python爬虫-使用xpath方法
#coding=utf-8 import re from lxml import etree import requests response = requests.get("http:// ...
- Cacti自定义脚本,监测Docker信息(Script/Command方式)
一 环境背景 监控主机A:192.168.24.231:被监控主机B:192.168.24.233 A/B主机,通过公私钥建立ssh连接 [操作B主机时不需要输入密码,详见笔记:http://app. ...
- 求csdn博客优良编辑方法
看见很多大牛的csdn博客编写的非常好,阅读体验也非常强.我就纳闷了,为啥我插公式也不行,插图片也不行呢... 插图片问题:图片不能复制招贴,否则在编辑的时候可以显示但是在发表之后就无法显示了.想要显 ...