[hdu3625] Examining the Rooms (第一类斯特林数)

传送门
感觉这题的思想还是挺巧妙的
反正我想出来的时候有种被自己机智到了的感觉233（其实是我平时太蠢，所以偶尔想出一道题就觉得自己棒棒的。。）

用$key[i]$表示第$i$个房间放的钥匙编号，那么$key$数组就是$1$到$n$的一个排列，所以总方案数就是$n!$。
这就是最后概率的分母咯

再求分子。
我们可以这样考虑，将$i$和$key[i]$连一条边，可以发现整张图变成了若干个环（包括自环）。
为什么会这样呢？
因为每个房间只有一把钥匙，每把钥匙对应唯一一个房间，所以每个点的入度与出度都为$1$，那这样一张图最后只能是若干个环。
对于每个环，肯定要暴力破开其中一个房间的门，所以最后破门数量就是环的数量。
那我们很容易能想到第一类斯特林，于是好像就是$\sum_{i=1}^{m}S_{n,i}$啊？
才怪。
由于房间$1$的门是不能被暴力破开的，所以要减去$1$成自环的情况
非常好，那分子就是$\sum_{i=1}^{m}S_{n,i}-S_{n-1,i-1}$

最后

$$p=\frac{\sum_{i=1}^{m}S_{n,i}-S_{n-1,i-1}}{n!}$$

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define LL long long
const int N=21;
const int p=1e9+7;
LL T,n,m,s[N][N];

void read(LL &x){
    char ch=getchar();x=0;
    for(;ch<'0'||ch>'9';ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
}

int main(){
    read(T);
    s[0][0]=1;
    for(int i=1;i<N;i++)
     for(int j=1;j<=i;j++)
      s[i][j]=s[i-1][j-1]+s[i-1][j]*(i-1);
    while(T--){
        read(n);read(m);
        LL sum=1;
        for(int i=1;i<=n;i++) sum*=i;
        LL fz=0;
        for(int i=1;i<=m;i++) fz+=s[n][i]-s[n-1][i-1];
        double ans=fz*1.0/sum;
        printf("%.4lf\n",ans);
    }
    return 0;
}