HDU3625 Examining the Rooms（第一类斯特林数）

题目链接

圆排列：把n个元素排在一个圆周上，如果旋转之后两个圆周上的排列一样，那么这两个排列相同 n个元素的圆排列=n-1个元素的排列数
第一类斯特林数S(n,m)表示把n个不同元素构成m个圆排列的方案数
由定义得S(n,0)=0,S(n,n)=1
第一类斯特林数有递推公式S(n,m)=S(n−1,m−1)+(n−1)×S(n−1,m)
证明：我们考虑把第n个元素放在什么位置。如果单独形成一个圆排列，则答案就是S(n-1,m-1).
如果插入到原来的m个圆排列中，那本质上就是在n-1个元素中选一个插到它后面去（定义“后面”为顺时针方向下一个数，定义为逆时针同理)，有n-1种选法，答案就是(n-1)S(n-1,m)

n个房间都有钥匙，1号门不能破门而入，要打开其他门有一种环状关系 比方说第i号房间有第j号房间的钥匙，那么你就要去找第i号房间的钥匙的所在地或者破门而入 - - - - 总之理解一下这个连环关系

ll dp[50][50];
int main(){
ll t=read();


while(t--){
ll n=read();
ll k=read();
ll ans=1;
for(int i=0;i<=n;i++){
dp[i][0]=0;
dp[i][i]=1;
 for(int j=1;j<=i;j++){
  dp[i][j]=dp[i-1][j-1]+(i-1)*dp[i-1][j];
 }
}

for(int i=1;i<=n;i++){
    ans=ans*i;
}
ll add=0;
for(int i=1;i<=k;i++){
   add=add+dp[n][i]-dp[n-1][i-1];
}
double rnm=add*1.0/ans;
printf("%.4f\n",rnm);
}
return 0;
}