Examining the Rooms 组合数学

题目：Examining the Rooms 题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3625

题目大意：给出N个房间，每个房间的钥匙随机放在某个房间内，概率相同。有K次炸门的机会，求能进入所有房间的可能性为多大。但是1号房间不能打开

我们可以先打开k个房间，也就是获得k把钥匙，我们要用这些钥匙打开所有房间。因为要求概率，所以可以用古典概型。
用可行方案数/总方案数。
先考虑总方案数即n的全排列，为n！
然后是可行方案数，那么我们先考虑不可行的原因：就是因为k个钥匙开房间会形成一个循环，最多k个循环，如果这些循环都没有包含n把钥匙，那这就不可行，反之，如果包含了，就可行。

那么问题就变成了n个物体组成1-k个非空环的方案数，这让我们很自然地想到第一类斯特林数。

什么是第一类斯特林数？请看我的博客：https://blog.csdn.net/luyan0408/article/details/98442689

知道这个了就好办了，答案就是从1到k枚举环的个数（因为有可能打开这几个房间中的钥匙处在一个循环），但是注意一号房间不能打开，就要排除单独成环的情况。

令s[i][j]表示i个物品成j个环的方案数
那么就要用s[i][j]-s[i-1][j-1],这样就好了，注意要枚举环的个数

要注意的是，具体是哪号门不能打开都可以，只要保证有一个自环即可，所以减去一个至少一个环的可能方案数。
下附AC代码：

#include<iostream>
#include<iomanip>
using namespace std;
int t; 
double s[300][300];//第一类斯特林数 
double fac[300];
int main(){
	fac[0]=1;
	for(int i=1;i<=200;i++){
		fac[i]=i*fac[i-1];
	}
	for(int i=1;i<=200;i++){
		s[i][i]=1;
	}
	for(int i=1;i<=200;i++){
		for(int j=1;j<i;j++){
			s[i][j]=s[i-1][j-1]+s[i-1][j]*(i-1);
		}
	} 
	cin>>t;
	while(t--){
		int k;
		int n;
		cin>>n>>k;
		double ans=0;
		for(int i=1;i<=k;i++){
			ans+=s[n][i]-s[n-1][i-1];//强制让一自成环 
		}
		cout<<fixed<<setprecision(4)<<ans/fac[n]<<endl;//注意保留四位小数 
	}
	
}