基础DP-第一类斯特林数-递推-hdu-3625-挑战程序设计2-q

开心，第一次学到了这么牛逼的东西。
好实用，好开心哦
第一类斯特林数用来解决这类问题
n个物品（可以区分）组成k个非空循环排列。
（注意哦，这些排列是有顺序的，比方说4个物品组成一个环
就有 1234 1243 1324 1342 1423 1432这么多种，画图的时候要注意啊。
）
斯特林数的递推公式如下
dp[0][0]=0;//废话
dp[n][n]=1;//废话
dp[n][0]=0;//废话，根本就不递推这种情况。
dp[n][k]=dp[n][k]+dp[n-1][k]*(n-1)//代码头里有讲解。。

https://vjudge.net/problem/HDU-3625
这个题有例外，要求1不能破门，所以要删掉1自环的情况，
其实就是减去一种自环情况就好了。
那他会不会减掉所有元素自环的情况呢，

#include<stdio.h>
#include<string.h>
#define N 21
typedef long long ll;
using namespace std;
/*第一类斯特林数。
 n个不相同的东西组成k个环（这tm还不是一般的环，）
 比方说4个东西 构成一个环
 有 1234 1243 1324 1342 1423 1432 6种情况。
 可以有以下递推方程。
 dp[0][0]=0;//
 dp[i][i]=1;
 dp[i][0]=1;//直接不递推这里就好，从1开始。
dp[n][k]=dp[n-1][k-1]+dp[n-1][k]*(k)
如果构成了k-1个环，那么最后一个要自环
如果构成
*/
ll dis[N][N];
ll gl[N];
void init()
{   gl[1]=1;
    for(int i=2;i<N;i++)
       gl[i]=gl[i-1]*i;//计算全排列；
     dis[1][1]=1;
     dis[0][0]=0;
     for(int i=2;i<N;i++)
        for(int j=1;j<=i;j++)
        dis[i][j]=dis[i-1][j-1]+(i-1)*dis[i-1][j];
    //不同代码不同实现但是思想是相同的哦；
  /*这里就是不同的初始化了，因为开始没有dis[i][i],所以在这里初始化了，并且
   发现是可以的，当i和j相等时，后一项为0（i比j小就为0了）
   所以这个也可以实现把 dis[i][i]设置成1
  */
}
int main()
{ int m,n,k;
    init();
    int t;
    scanf("%d",&t);
    /*for(int i=0;i<5;i++)
    {for(int j=0;j<5;j++)
        cout<<dis[i][j]<<" ";
         cout<<endl;
    }*/
   while(t--)
   {  scanf("%d%d",&m,&n);
       ll sum=0;
       for(int i=1;i<=n;i++)
          sum=dis[m][i]-dis[m-1][i-1]+sum;
          //去除1自环的情况。
          //cout<<sum<<endl;
       printf("%.4lf\n",1.0*sum/gl[m]);
   }
    return 0;
}