HDU 3625 Examining the Rooms：第一类stirling数

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

题意：

　　有n个房间，每个房间里放着一把钥匙，对应能开1到n号房间的门。

　　除了1号门，你可以踹开任意一扇门（不用钥匙），但你最多只能踹k次。

　　问你能将所有门打开的概率。

题解：

　　· P(打开所有门) = 能打开所有门的钥匙放置情况数 / 钥匙放置的总情况数

　　· 钥匙放置的总情况数 = n!

　　

　　那么考虑下能打开所有门的钥匙放置情况数。。。

　　由于每个房间里有且只有一把钥匙，所以如果将每个房间连向房间内钥匙对应的房间，会得到一个有向图，并且由若干个独立的环组成。

　　

　　所以，只要将一个环内的任意一扇门踹开，就能打开这个环上的所有房间。

　　如果不考虑1号门，那么要算的就是n个元素组成1~k个环排列的情况数。

　　第一类Stirling数的定义啊！

　　递推式：s(n,k) = s(n-1,k-1) + (n-1)*s(n-1,k)

　　

　　So...

　　· 能打开所有门的钥匙放置情况数 = ∑S(n,i)　(1<=i<=k)

　　

　　考虑到1号门不能踹，也就是1不能独立成环，它的情况数就等于用剩下n-1的元素组成i-1个环的情况数。

　　· 能打开所有门的钥匙放置情况数 = ∑( S(n,i)-S(n-1,i-1) )

　　所以答案为：∑( S(n,i)-S(n-1,i-1) ) / n!　　(1<=i<=k)

AC Code:

 // s(n,k) = s(n-1,k-1) + (n-1)*s(n-1,k)
 // s(n,0) = 0    s(n,n) = 1
 // P = sigma(s[n][i]-s[n-1][i-1])/fact(n)    1<=i<=k

 #include <iostream>
 #include <stdio.h>
 #include <string.h>
 #define MAX_N 25
 #define MAX_K 25

 using namespace std;

 int n,k,t;
 long long sum;
 long long s[MAX_N][MAX_K];
 long long fact[MAX_N];

 void stirling()
 {
     memset(s,,sizeof(s));
     for(int i=;i<MAX_N;i++)
     {
         s[i][i]=;
         for(int j=;j<i;j++)
         {
             s[i][j]=s[i-][j-]+(i-)*s[i-][j];
         }
     }
 }

 void cal_fact()
 {
     fact[]=;
     for(int i=;i<MAX_N;i++)
     {
         fact[i]=fact[i-]*i;
     }
 }

 int main()
 {
     stirling();
     cal_fact();
     cin>>t;
     for(int cas=;cas<=t;cas++)
     {
         cin>>n>>k;
         sum=;
         for(int i=;i<=k;i++)
         {
             sum+=s[n][i]-s[n-][i-];
         }
         printf("%.4f\n",(double)sum/fact[n]);
     }
 }