题目分析
首先我们可以利用坐标来观察,我们发现x+y坐标之和为偶数的颜色可以相同,奇数的也可以相同,所以我们就可以将图分成2部分,然后进行染色。当然这里我们用dp方程打表推出染色所需要用到的各种方法数。dp[i][j]表示用j种颜色染i个格子的方案数,注意这里面j中颜色必须全部用完,因此状态转移方程为dp[i][j] = dp[i-1][j]*j+dp[i-1][j-1]*j;预处理一下即可计算结果。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define LL long long
const int mod = 1e9+7;
const int maxn = 205;
LL C[maxn][maxn], dp[maxn][maxn];
//dp[i][j]表示用j种颜色染i个格子且每种颜色都出现过的方案数,则dp[i][j] = dp[i-1][j]*j + dp[i-1][j-1]*j;
void init(){
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;
}
for(int i = 1; i < maxn; i++){
for(int j = 1; j < maxn; j++){
if(j == 1) dp[i][j] = 1;
else dp[i][j] = dp[i-1][j]*j + dp[i-1][j-1]*j;
dp[i][j] %= mod;
}
}
}
int main(){
init();
int T, m, n, k;
scanf("%d", &T);
for(int kase = 1; kase <= T; kase++){
scanf("%d%d%d", &m, &n, &k);
if(!n && !m){
printf("Case %d: %d\n", kase, k);
continue;
}
int x = ((n+1)*(m+1)+1)/2;
int y = (n+1)*(m+1)-x;
LL ans = 0;
for(int i = 1; i < k; i++)
for(int j = 1; i+j <= k; j++)
ans = (ans + (C[k][i]*C[k-i][j])%mod *((dp[x][i]*dp[y][j])%mod))%mod;
printf("Case %d: %lld\n", kase, ans);
}
return 0;
}