今天下午多校赛的痛苦经历……
下午的时候在选出一些数字放入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;
}