Killer Names
题目:(题目链接)
思路数学题
考验的是排列组合
思路:
因为一个字符不可以同时出现在名和性中,那么就先选择字符(字符可以不用完)
选完字符后乘以其所有的情况,然后把其加起来,过程中要mod
情况是提前先算出来了,下面代码中a数组存的就是此过程
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <map>
#include <set>
#include <algorithm>
#include <vector>
#include <bitset>
#include <stack>
#include <queue>
#include <functional>
using namespace std;
#define LL long long
const int INF=0x3f3f3f3f;
LL a[2005];
const LL mod=1e9+7;
LL myypow(LL x,LL y)//快速幂
{
LL ans=1;
while(y>=1)
{
if(y&1)
ans=(ans*x)%mod;
x=(x*x)%mod;
y>>=1;
}
return ans;
}
LL c[2022][2022];
void chss()//计算排列组合···Cij
{
c[0][0] = 1;
for(int i = 1; i < 2022; i++)
{
c[i][0] = c[i][i] = 1;
for(int j = 1; j < i; j++)
c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
}
}
int main()
{
int T;
LL m,n;
chss();//把所有的排列组合存起来
for(scanf("%d",&T); T--;)
{
scanf("%I64d%I64d",&n,&m);
for(int i=1; i<m; i++)
{
a[i]=myypow(i,n);
for(int j=1; j<i; j++)//a[i]表示用i个字符组成名(或者性)有多少种情况
a[i]=((a[i] - c[i][j]*a[j] % mod) % mod + mod) % mod;
}
LL ans = 0;
for(int i=1; i<m; i++)
{
for(int j = 1; j <= m - i; j++)
ans = (ans + c[m][i] * a[i] % mod * c[m - i][j] % mod * a[j] % mod) % mod;
}
printf("%I64d\n", (ans % mod + mod) % mod);//加上mod是防止出现负数
}
return 0;
}