[AHOI2009]中国象棋——[计数DP]

【题目描述】

这次小可可想解决的难题和中国象棋有关，在一个N行M列的棋盘上，让你放若干个炮（可以是0个），使得没有一个炮可以攻击到另一个炮，请问有多少种放置方法。大家肯定很清楚，在中国象棋中炮的行走方式是：一个炮攻击到另一个炮，当且仅当它们在同一行或同一列中，且它们之间恰好 有一个棋子。你也来和小可可一起锻炼一下思维吧！

【输入格式】

一行$n,m(n,m\le 100)$

【输出格式】

一行， 表示方案数（对9999973取模）

$SampleInput$

1 3

$SampleOutput$

7

【题意分析】
感谢大佬：顾z

非常像状压，但是100你是压不了的，考虑计数
根据规则转化题意：求每行每列棋子不超过两个的方案数

$dp[i][j][k]$表示前 $i$ 行，有 $j$ 列放了一个棋子，有 $k$ 列放了两个棋子的方案数
那么答案就是$\sum _{i=0}^m\sum _{j=0}^{m}dp[n][i][j]$

显而易见当前第 $i$ 行最多放 $2$ 个棋子，那么我们可以分类讨论了：

• 初始情况
  $$dp[0][0][0]=1$$
• 不放任何棋子
  继承上一行状态
  $$dp[i][j][k]+=dp[i-1][j][k]$$
• 放一个棋子
  放在没有棋子的一列，会使 $j$ 增加1，所以要寻找 $j-1$ 列的状态
  因为放在空列，所以有$m-(j-1)-k$种放法
  $$dp[i][j][k]+=dp[i-1][j-1][k]\ast (m-(j-1)-k)$$放在有一个棋子的一列，会使 $k$ 增加1， $j$减少1，所以要寻找$j+1,k-1$的状态
  那么有$j+1$种放法。
  $$dp[i][j][k]+=dp[i-1][j+1][k-1]\ast (j+1)$$
• 放两个棋子
  放在两个空列，会使$j$增加2，所以要寻找$j-2$的状态
  在空列中选两个，即$C_{(m-(j-2)-k)}^2$种放法
  $$dp[i][j][k]+=dp[i-1][j-2][k]\ast c(m-(j-2)-k,2)$$放在一个空列，一个有一个棋子的列，使$j$增加1又减少1，k增加1，所以要寻找$j,k-1$的状态
  根据乘法原理有$j\ast (m-j-(k-1))$种放法
  $$dp[i][j][k]+=dp[i-1][j][k-1]\ast j\ast (m-j-(k-1))$$放在两个有一个棋子的列，使$k$增加2，$j$减少2，所以要寻找$j+2,k-2$的状态
  在有一个棋子的列种选两个，即$C_{j+2}^2$种放法
  $$dp[i][j][k]+=dp[i-1][j+2][k-2]\ast c(j+2,2)$$

注意边界：对于需要减少的参数，要判断是否大于等于0
别忘了取%

Code:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <cstring>
#include <algorithm>
#define rep(x,y,z) for (register int x = y; x <= z; x++)
#define MAXN 200
#define qy 9999973
using namespace std;

long long dp[MAXN][MAXN][MAXN];

inline long long C (int x) {
	return (long long) (x * (x - 1) / 2) % qy;
} 

int main () {
	int n, m; scanf ("%d %d", &n, &m), dp[0][0][0] = 1;
	rep (i, 1, n) rep (j, 0, m) rep (k, 0, m - j) {
		dp[i][j][k] = dp[i-1][j][k];
		if (j >= 1) dp[i][j][k] += dp[i-1][j-1][k] * (m-(j-1)-k);
		if (k >= 1) dp[i][j][k] += dp[i-1][j+1][k-1] * (j+1);
		if (j >= 2) dp[i][j][k] += dp[i-1][j-2][k] * C (m-(j-2)-k);
		if (k >= 2) dp[i][j][k] += dp[i-1][j+2][k-2] * C (j+2);
		if (k >= 1) dp[i][j][k] += dp[i-1][j][k-1] * j * (m-j-(k-1));
		dp[i][j][k] %= qy;
	}
	long long ans = 0;
	rep (i, 0, m) rep (j, 0, m) ans += dp[n][i][j], ans %= qy;
	printf ("%lld", ans);
	return 0;
}