斯特林数 HDU 3625

/*
 斯特林数第一类
 n个元素的集合分成k个环排列的方法
 s[n][0] = 0;
 s[1][1] = 1;
 s[n][k] = s[n-1][k-1]+(n-1)*s[n-1][k];
 当放入第n个元素，n可以为单独的循环排列，当n-1个元素已经占据了k个环排列，n可以插入n-1个元素的左边
 该题要去除1号门的单独排列，所以。。。
*/

#include<cstdio>
#include<iostream>
using namespace std;
typedef long long LL;
const int maxn = 25;
LL f[maxn];
LL st[maxn][maxn];

void init() {
    f[1]=1;
    for (int i=2; i<maxn; i++)
        f[i] = i*f[i-1];

    for (int i=0; i<maxn; i++)
        st[i][0] = 0;
    st[1][1] = 1;
    for (int i=1; i<maxn; i++) {
        for (int j=1; j<=i; j++)
            if (i == j)
                st[i][j] = 1;
            else
                st[i][j] = st[i-1][j-1]+(i-1)*st[i-1][j];
    }

    for (int i=1; i<maxn; i++) {
        for (int j=1; j<maxn; j++)
            st[i][j] = abs(st[i][j]);
    }
}

int main() {
    int T;
    scanf("%d", &T);
    init();
    while (T--) {
        int n, k;
        scanf("%d%d", &n, &k);
        LL sum = 0;
        for (int i=1; i<=k; i++)
            sum += st[n][i]-st[n-1][i-1];
        printf("%.4lf\n",1.0*sum/f[n]);
    }
    return 0;
}