hdu 6143 Killer Names(容斥原理)

枚举要为这两段分配的颜色数目分别为  i,j i,j ，则在第一段总共有  (mi) (mi) 种选取方案，在第二段总共有  (m−ij) (m−ij) 种选取方案。

而在每段内部，我们设  F(n,x) F(n,x) 为长度为  n n 的格子使用  x x 种颜色（等于  x x ）染色的方案数。

则根据容斥原理F(n,x)=xn−(m1)(x−1)n+(m2)(x−2)n−(m3)(x−3)n+...F(n,x)=xn−(m1)(x−1)n+(m2)(x−2)n−(m3)(x−3)n+...

//注意细节，优化很重要，很重要，在1s多点一直在卡，卡到boom

#include<bits/stdc++.h>
typedef long long LL;
using namespace std;
const int maxn =2005;
const int mod = 1e9+7;

LL res[maxn];
LL C[maxn][maxn];
LL mul[maxn][maxn];

LL mult(LL a,LL b)
{
    if(mul[a][b])return mul[a][b];
    LL ans=1,aa=a,bb=b;
    a%=mod;
    while(b)
    {
        if(b&1)ans=(ans*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    mul[aa][bb]=ans;
    return ans;
}

void init()
{
    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-1]+C[i-1][j])%mod;
}

LL getget(LL n,LL m)
{
    LL cnt=0;
    for(int k=0; k<=m; k++)
    {
        if(k&1)cnt=(cnt-C[m][k]*mult(m-k,n)%mod)%mod;
        else cnt=(cnt+C[m][k]*mult(m-k,n)%mod+mod)%mod;
    }
    return cnt;
}

template <class T>
inline void scan_d(T &ret)
{
    char c;
    ret = 0;
    while ((c = getchar()) < '0' || c > '9');
    while (c >= '0' && c <= '9')
        ret = ret * 10 + (c - '0'), c = getchar();
}

void Out(LL a)
{
    if (a < 0)
    {
        putchar('-');
        a = -a;
    }
    if (a >= 10)
        Out(a / 10);
    putchar(a % 10 + '0');
}

int main()
{
    init();
    int T;
    scan_d(T);
    while(T--)
    {
        memset(res,0,sizeof(res));
        int n,m;
        LL ans=0;
        scan_d(n);
        scan_d(m);
        for(int i=1; i<m; i++)
            res[i]=getget(n,i);
        for(int i=1; i<m; i++)
            for(int j=i; j<=m-i; j++)
            {
                LL aa=C[m][i]*res[i]%mod;
                LL bb=C[m-i][j]*res[j]%mod;
                LL cnt = (aa*bb)%mod;
                if(i!=j) cnt = (cnt+cnt)%mod;
                ans=(cnt+ans)%mod;
            }
        Out(ans);
        putchar('\n');
    }
    return 0;
}

F(n,x)=