hdu6143-多校8&三种方法-组合数|递推|容斥-Killer Names

http://acm.hdu.edu.cn/showproblem.php?pid=6143
这道题这么多写法qwq
真的好强悍
1 递推。当时也想来着，不过我想的是打表那种，但其实m和n不一样，结果也是不一样的。在有限的颜色下。
i为 当前段，j为使用过的颜色
满足递推方程dp[i][j]=dp[i-1][j]j+dp[i-1][j-1](n-j+1)…
后面计算另一个的时候，就反着求。
开始没开dp数组没开longlong。。虽然1e9+7，但是乘起来可能就爆longlong了，出现负数，皆大悲伤。

#include <bits/stdc++.h>
using namespace std;
/* 方法有三种，
 1 容斥原理
 2 组合数递推。（我们没推出来这个，GG）
  3 dp，老实说。这个我感觉最可以写。
  这个最扯淡。。 直接dp出 i行j列的情况，
  然后取一个相反的。我tm。。
*/
typedef long long ll;
const int maxn=2002;
const  int  mod=1e9+7;
ll quickpow(ll m,ll n){//快速求幂
   ll ans=1;
   while(n>0){
       if(n&1) ans=(ans*m)%mod;
        m=m*m%mod;
        n>>=1;
   }
   return ans;
}
ll  dp[2002][2002];
int main()
{    int t,m,n;
     scanf("%d",&t);
  int a=1e9+7;
     while(t--){
           scanf("%d%d",&m,&n);
           dp[1][1]=n;
           for(int i=2;i<=m;i++){
               for(int j=1;j<=n;j++)
                   dp[i][j]=((dp[i-1][j]*j)%mod+(dp[i-1][j-1]*(n-j+1))%mod)%mod;
        }
           ll ank=0;
           for(int i=1;i<n;i++){
             ank=(ank+1ll*dp[m][i]*quickpow(1ll*(n-i),1ll*m))%mod;
           }
          printf("%lld\n",ank);
     }
    return 0;
}

2 组合数，不推荐这个，不好写qwq。
3 容斥（离散里的。）
容斥真的好高大上qwq
1 先求 F【n,i】 n个板子，用最多i种颜色（递推也可以好么qwq,但是感觉递推还是没有容斥猛）
2 然后在 c(n,i)F[m,i]+c(n-i,j)F[m,j].
注意这个枚举的时候没有策略。
但是最好这样写，否则会超时。。

#include<bits/stdc++.h>
using namespace std;
/* 根据容斥原理。。
（我觉得递推更加好一点qwq）。
n长的板子  用最多k种颜色染（）
 把 快速幂 和组合数给暴力打表一下。
？？我们比赛的时候在想什么？？
*/
typedef long long ll;
const int maxn=2002;
const int mod=1e9+7;
ll c[maxn][maxn];
ll mi[maxn][maxn];
ll quickpow(ll m,ll n){
    ll ans=1;
    while(n){
        if(n&1) ans=ans*m%mod;
        m=m*m%mod;
        n>>=1;
    }
   return ans;
}
int zhs(){
     c[0][0]=1;
     for(int i=1;i<maxn;i++){
         for(int j=0;j<=i;j++){
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
          //求解组合数。qwq
         }
     }
}
int m,n;
ll ans[maxn];
int main()
{   int t;
    for(int i=1;i<maxn;i++)
       mi[i][1]=quickpow(i,1);
    for(int i=1;i<maxn;i++){
        for(int j=2;j<maxn;j++)
            mi[i][j]=(mi[i][j-1]*1ll*i)%mod;
    }
    zhs();
    scanf("%d",&t);
    while(t--){
          scanf("%d%d",&m,&n);
          ll all=0;
          memset(ans,0,sizeof(ans));
          for(int i=1;i<n;i++){
              for(int j=0;j<=i;j++)
                  ans[i]=(ans[i]+c[i][j]*mi[i-j][m]*(int)pow(-1,j)+mod)%mod;
          }
          for(int i=1;i<n;i++){
              for(int j=i;j<=n-i;j++){
                  ll aa=(ans[i]*c[n][i]+mod)%mod;
                   ll bb=(ans[j]*c[n-i][j]+mod)%mod;

                   ll kk=(aa*bb)%mod;
                  if(i!=j) kk=(kk+kk)%mod;
                 all=(all+kk)%mod;
              }
          }
          printf("%lld\n",all);

     }

    return 0;
}