【Stirling Number I】

hdu 4372 Count the Buildings

推荐这位小哥的，我觉得人家说的灰常的好。

注意数据范围，n,f,b均在(0,2000]范围内，而第一斯特林数的数组范围却是s[2000+5][2000+5]的（你要是开4000会内存超限），所以要加一个对（f+b-2）的判断，否则C++过不了，G++是能过的。亲身实践，虽然不知道G++为什么这么鲁棒...

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
const ll MOD = 1000000007;
const int maxn = 2005;
ll s[maxn][maxn];
ll C[maxn][maxn];

void get_Combination_Number()
{
    for(int i = 1; i <= 2000; i++)
    {
        C[i][0] = C[i][i] = 1;
        for(int j = 1; j < i; j++)
            C[i][j] = (C[i-1][j-1]+C[i-1][j])%MOD;
    }
}

void get_Stirling_Number_I()
{
    s[0][0] = 1;
    for(int i = 1; i <= 2000; i++)
    {
        s[i][0] = 0;
        s[i][i] = 1;
        for(int j = 1; j < i; j++)
        {
            s[i][j] = (s[i-1][j-1] + (i-1)*s[i-1][j]) % MOD;
        }
    }
}

int main()
{
    get_Combination_Number();
    get_Stirling_Number_I();
    int T;
    scanf("%d", &T);
    while(T--)
    {
        int n, f, b;
        scanf("%d%d%d", &n, &f, &b);
        ll ans = 0;
        if(f+b-2 <= n-1)
            ans = (s[n-1][f+b-2] * C[f+b-2][f-1]) % MOD;
        else
            ans = 0;
        printf("%lld
", ans);
    }
    return 0;
}