HDU 6143 Killer Names

T题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6143
题意：
m 种颜色需要为两段长度为 n 的格子染色，且这两段之间不能出现相同的颜色，问总共有多少种情况。
思路：

截图于：http://blog.csdn.net/qq_28954601/article/details/77337256
大佬的思路非常清晰。 我也不用多说了。
就是调试代码的时候，不懂打的表是什么意思，看了看其实就是利用组合数性质：C（M-1，N-1）+C（M-1，N）=C（M，N），快速批量求出组合数。ps：高中的知识都还给老师了。T_T

#include <cstdio>
#include <algorithm>
#include <iostream>
using namespace std;

const int mod = 1e9 + 7;
const int maxn = 2000 + 10;
//快速幂
int power(int x, int times)
{
    int rt = 1;
    int base = x;
    while (times)
    {
        if (times & 1) rt = 1LL * rt * base % mod;
        base = 1LL * base * base % mod;
        times >>= 1;
    }
    return rt;
}
int c[maxn][maxn];
int f[maxn];
int main()
{
    //求的是组合数
    for (int i = 0; i < maxn; 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;
            /*if(i<=100)
            {
                printf("%d %d %d\n",i,j,c[i][j]);
            }*/
        }
    }

    int cases;
    int n, m;
    scanf("%d", &cases);
    while (cases--)
    {
        scanf("%d%d", &n, &m);
        int ans = 0;
        for (int i = 1; i < m; i++)
        {
            if (i > n) break;
            f[i] = power(i, n);
            for (int j = 1; j < i; j++)
            {
                f[i] = (f[i] - 1LL * c[i][j] * f[j]) % mod;
            }
            ans = (ans + 1LL * f[i] * c[m][i] % mod * power(m - i, n)) % mod;
        }
        ans = (ans + mod) % mod;
        printf("%d\n", ans);
    }
    return 0;
}

这题还可以用第二类Stirling数来做。
第二类Stirling数 S(p,k) ，S(p,k)的一个组合学解释是：将p个物体划分成k个非空的不可辨别的（可以理解为盒子没有编号）集合的方法数。
k!S(p,k)是把p个人分进k间有差别(如：被标有房号）的房间(无空房）的方法数。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=2e3+5;
const int mod=1e9+7;
ll str[maxn][maxn];
ll a[maxn][maxn],f[maxn];
///快速幂
ll q(ll x,ll n)
{
    ll res=1;
    while(n)
    {
        if(n&1) res=res*x%mod;
        n>>=1;
        x=x*x%mod;
    }
    return res;
}
void init()
{
    for(int i=1;i<=2001;i++)
    {
        str[i][0]=0;str[i][1]=1;
        for(int j=2;j<=i;j++)//第二类Stirling数
        {
            str[i][j]=(j*str[i-1][j]%mod+str[i-1][j-1])%mod;
        }
    }
    a[0][0]=1;
    for(int i=1; i<=maxn; i++)
    {
        a[i][0]=a[i][i]=1;
        for(int j=1; j<i; j++)//组合数
        {
            a[i][j]=(a[i-1][j-1]+a[i-1][j])%mod;
        }
    }
    f[1]=1;
    for(int i=2; i<maxn; i++)
    {
        f[i]=(f[i-1]*i)%mod;
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    init();
    while(t--)
    {
        ll n,m;
        scanf("%lld %lld",&n,&m);
        ll ans=0;
        if(m<=n)
            for(int i=1; i<m; i++)
            {
                ans+=(((a[m][i]*f[i])%mod*str[n][i])%mod*q(m-i,n))%mod;
                ans%=mod;
            }
        else
        {
            for(int i=1; i<=n; i++)
            {
                ans+=(((a[m][i]*f[i])%mod*str[n][i])%mod*q(m-i,n))%mod;
                ans%=mod;
            }
        }
        printf("%lld\n",ans);
    }
}