生活總是那麼操蛋，(a+b)%mod寫成了(a+b%mod)，找錯找了一個多小時。。
思路：
m種字元，名字分兩部分，都能放下n個字元。先取出名字的一部分，即n個位置。對於這n個位置，能放i種字元，i屬於[1, min(m,n)]，就是n個元素放進k個盒子裡，並且盒子不同，所以要乘階乘，這對應第二類斯特林數，對於某個i，共Cim∗S[n][i]∗fac[i]種情況，剩下m-i個字元就放到名字的另一部分，方法同上。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL mod = 1
 
e9+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;
}