hdu 6143 Killer Names（第二类斯特林数）

生活总是那么操蛋，(a+b)%mod写成了(a+b%mod)，找错找了一个多小时。。
思路：
m种字符，名字分两部分，都能放下n个字符。先取出名字的一部分，即n个位置。对于这n个位置，能放i种字符，i属于[1, min(m,n)]，就是n个元素放进k个盒子里，并且盒子不同，所以要乘阶乘，这对应第二类斯特林数，对于某个i，共$C^{i}_{m}*S[n][i]*fac[i]$种情况，剩下m-i个字符就放到名字的另一部分，方法同上。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL mod = 1e9+7;
const int MAXN = 2010;

LL S[MAXN][MAXN];
LL C[MAXN][MAXN];
LL Fac[MAXN];

void init()
{
    Fac[0] = 1;
    C[0][0] = 1;
    S[0][0] = 0;
    for(LL i = 1; i < MAXN; ++i)
    {
        Fac[i] = (Fac[i-1]*i)%mod;
        S[i][i] = 1;
        S[i][0] = 0;
        C[i][i] = 1;
        C[i][0] = 1;
        for(LL j = 1; j < i; ++j)
        {
            C[i][j] = (C[i-1][j] + C[i-1][j-1])%mod;
            S[i][j] = (S[i-1][j-1] + j*S[i-1][j])%mod;
        }
    }
}

int main()
{
    init();
    int T,n,m;
    scanf("%d",&T);
    while(T--)
    {
        LL res = 0,half;
        int h1,h2;
        scanf("%d %d",&n,&m);
        for(int i = 1; i < m && i <= n; ++i)
        {
            h1 = i;
            half = (C[m][i]*S[n][i])%mod*Fac[i]%mod;
            h2 = m-i;
            for(int j = 1; j <= h2 && j <= n; ++j)
                res = (res + ((C[h2][j]*S[n][j]%mod)*Fac[j]%mod)*half%mod)%mod;
        }
        printf("%lld\n",res);
    }
    return 0;
}