题目链接:
http://acm.hdu.edu.cn/showproblem.php?
pid=3625
题目大意:
有N个房间,每一个房间的要是随机放在某个房间内,概率同样。有K次炸门的机会。
求能打开全部房间门,进入全部房间的概率有多大。
解题思路:
门和钥匙的相应关系出现环。打开一个门后,环内的门都能够打开。
也就意味着:
N个房间的钥匙与门形成1~K个环的概率有多大。
也就是求N个元素。构成K个以内的环,且1不成自环的概率。
N个元素形成K个环的方法数是第一类stirling数 S(N。K)。
N个元素形成K个环。且1成自环的方法数是S(N-1,K-1)。
则N个元素形成K个环。且1不成自环的方法数是S(N。K) - S(N-1,K-1)。
要是随机放的总的方法数为N!。
则概率P(N,K)为( S(N,K) - S(N-1,K-1) + S(N。K-1) - S(N-1。K-2) + … +
S(N,1) - S(N。0) ) / N!
AC代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define LL long long
using namespace std; LL F[25],Stirling[25][25]; void Solve()
{
F[0] = 1;
for(int i = 1; i <= 20; ++i) //阶乘数组
F[i] = i*F[i-1]; for(int i = 1; i <= 20; ++i) //算出第一类stirling数
Stirling[i][0] = 0;
Stirling[1][1] = 1;
for(int i = 1; i <= 20; ++i)
{
for(int j = 1; j <= i; ++j)
{
if(i == j)
Stirling[i][j] = 1;
else
Stirling[i][j] = Stirling[i-1][j-1] + (i-1)*Stirling[i-1][j];
}
} for(int i = 1; i <= 20; ++i) //取绝对值
{
for(int j = 1; j <= 20; ++j)
{
Stirling[i][j] = abs(Stirling[i][j]);
}
}
} int main()
{
Solve();
int T,N,K;
scanf("%d",&T);
while(T--)
{
scanf("%d %d",&N,&K);
LL sum = 0;
for(int i = 1; i <= K; ++i)
sum += Stirling[N][i] - Stirling[N-1][i-1];
printf("%.4f\n",(double)sum / (double)F[N]); //.4lf输出不了正确结果
} return 0;
}