BZOJ1801: [Ahoi2009]chess 中国象棋

http://www.lydsy.com/JudgeOnline/problem.php?id=1801

Description

在N行M列的棋盘上，放若干个炮可以是0个，使得没有任何一个炮可以攻击另一个炮。 请问有多少种放置方法，中国像棋中炮的行走方式大家应该很清楚吧.

/*递推，有点像扫雷Mine
每一行每一列最多放两个炮 
考虑状压记录每一列已经放了几个炮了
N M<=100 TLE MLE
发现状压记录的信息太多了，有些可以简化 
其实我们不用管到底哪些列有一个有两个炮，
只用看到现在有几列是一个炮，有几列没有炮、两个炮
转移的时候只用乘一个组合数即可 */
#include <cstdio>
#include <algorithm>
typedef long long LL;
const int MAX_N=105;
const int MOD=9999973;
int f[MAX_N][MAX_N][MAX_N];
int n,m;
int main(){
	scanf("%d%d",&n,&m);
	if (n<m) n^=m,m^=n,n^=m;
	f[0][m][0]=1;
	for (int i=0;i<n;++i){
		for (int j=0;j<=m;++j){
			for (int k=0;k<=m-j;++k){
				//不放
				(f[i+1][j][k]+=f[i][j][k])%=MOD; 
				//放 1:
				if (j-1>=0) (f[i+1][j-1][k+1]+=f[i][j][k]*j)%=MOD;//0->1
				if (k-1>=0) (f[i+1][j][k-1]+=f[i][j][k]*k)%=MOD;//1->2
				//放 2:
				if (j-2>=0) (f[i+1][j-2][k+2]+=((LL)f[i][j][k]*((LL)j*(j-1)/2))%MOD)%=MOD;//0,0->1,1
				if (j-1>=0&&k>0) (f[i+1][j-1][k]+=((LL)f[i][j][k]*j*k)%MOD)%=MOD;//0,1->1,2
				if (k-2>=0) (f[i+1][j][k-2]+=(f[i][j][k]*((LL)k*(k-1)/2))%MOD)%=MOD;//1,1->2,2
				}
			}
		}
	int ans=0;
	for (int j=0;j<=m;++j)
	 for (int k=0;k<=m-j;++k)
	 	(ans+=f[n][j][k])%=MOD;
	printf("%d\n",ans);
}