hdu3625-第一类斯特林数

第一类斯特林数S(n, m)是将n个物品分成m个非空循环排列的方法数

对于第n个物品，可以单独构成一个非空循环排列，这样前n-1个物品构成m-1个非空循环排列，方法数是 S(n - 1, m - 1),

也可以前n - 1个物品构成m个非空循环排列，第n个物品插入第i个物品的左边，方法数是 (n - 1) * S(n - 1, m),

所以递推公式为S(n, m) = S(n - 1, m - 1) + (n - 1) * S(n - 1, m)，S(i, 0) = 0, S(i, i) = 1

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

题意：有n个房间都被锁上了，开门的钥匙被放在房间里面(每个房间一个)，第一个房间的门不能被破坏，求最多破坏k个门，有可能将所有房间都打开的概率。

因为第一个房间的门不能被破坏，所以第一个房间不能单独成一个非空循环排列，所以破坏i个门打开所有房间的方法数为S(n, i) - S(n - 1, i - 1),

求出破坏1,2,3,...,k个门的方法数之和，求出n的全排列，两者比值即为所求。

代码：

# include <iostream>
# include <algorithm>
# include <cstdio>
# include <cstring>
using namespace std;
typedef long long ll;
const int maxn = 20 + 1;
int n, k;
ll S[maxn][maxn];
ll A[maxn][maxn];

void init() {
    memset(S, 0, sizeof S);
    S[1][1] = 1;
    for (int i = 2; i < maxn; ++i) {
        for (int j = 1; j <= i; ++j) {
            S[i][j] = S[i - 1][j - 1] + (i - 1) * S[i - 1][j];
        }
    }
    for (int i = 1; i < maxn; ++i) {
        A[i][1] = i;
        for (int j = 2; j <= i; ++j) {
            A[i][j] = A[i][j - 1] * (i - j + 1);
        }
    }
}

int main(void)
{
    init();
    int T; scanf("%d", &T);
    while (T-- && scanf("%d %d", &n, &k)) {
        ll cnt2 = A[n][n];
        ll cnt1 = 0;
        for (int i = 1; i <= k; ++i) {
            cnt1 += S[n][i] - S[n - 1][i - 1];
        }
        printf("%.4f\n", cnt1 * 1.0 / cnt2);
    }

    return 0;
}