hdu 6143 Killer Names （组合数，递推）

链接

今天下午多校赛的痛苦经历……

下午的时候在选出一些数字放入n有多少种方式计算时，只想着怎么用排列组合的方式计算出来，然后发现计算量简直爆炸……然后就gg……
后来的发现递推这种方法。首先假如用x个字母去填n，那所有的情况就是x^n，然后分别减去只用了1、2…到x-1种字母的情况，就是x个字母全用上的组合方式。
处理出来这些之后，答案就是第一个单词用1~m-1种，第二个字母用1~m-i种方案的组合累加。

注意取模的地方很多……

#include <bits/stdc++.h>
#define ll long long
using namespace std;

const int maxn= 2e3+10;
const int mod=1e9+7;
ll C[maxn][maxn],dp[maxn];
int n,m,t;
void get_c()
{
    C[1][0] = C[1][1] = 1;
    for (int i = 2; i < maxn; i++){
        C[i][0] = 1;
        for (int j = 1; j < maxn; j++)
            C[i][j] = (C[i - 1][j] + C[i - 1][j - 1])%mod;
    }
}

ll quick_power(ll a,ll b)
{
    ll ans=1;
    ll tmp=a%mod;
    while(b)
    {
        if(b&1)
        {
            ans=(ans*tmp)%mod;
        }
        tmp=(tmp*tmp)%mod;
        b>>=1;
    }
    return ans;
}

int main()
{
    cin>>t;
    get_c();
    while(t--)
    {
        cin>>n>>m;
        dp[1]=1;
        for(int i=2;i<=m;i++)
        {
            dp[i]=quick_power(i,n);
            for(int j=1;j<i;j++)
                dp[i]=(dp[i]-dp[j]*C[i][j]%mod+mod)%mod;
        }
        ll ans=0;
        for(int i=1;i<m;i++)
        {
            for(int j=1,k=m-i;j<=k;j++)
            ans=(ans+C[m][i]*dp[i]%mod*C[k][j]%mod*dp[j]%mod)%mod;
        }
        cout<<ans<<endl;

    }
    return 0;
}