1801: [Ahoi2009]chess 中国象棋

#include<iostream>
#include<cstdio>
using namespace std;
inline long long c(long long x){return x*(x-1)/2;}
int n,m;
long long ans,f[101][101][101];
int main(){
	scanf("%d%d",&n,&m);
	f[1][0][0]=1;f[1][1][0]=m;f[1][2][0]=c(m);
	for(int i=2;i<=n;i++)
		for(int j=0;j<=m;j++)
			for(int k=0;k<=m-j;k++){
				f[i][j][k]=f[i-1][j][k];
				if(j)f[i][j][k]+=f[i-1][j-1][k]*(m-j-k+1);
				if(j>1)f[i][j][k]+=f[i-1][j-2][k]*c(m-j-k+2);
				if(k&&j<m)f[i][j][k]+=f[i-1][j+1][k-1]*(j+1);
				if(j&&k)f[i][j][k]+=f[i-1][j][k-1]*(m-j-k+1)*j;
				if(k>1&&j<m-1)f[i][j][k]+=f[i-1][j+2][k-2]*c(j+2);
				f[i][j][k]%=9999973;
			}
	for(int i=0;i<=m;i++)
		for(int j=0;j<=m-i;j++)	
			ans+=f[n][i][j];
	printf("%lld",ans%9999973);
	return 0;	
}

f[i][j][k]表示前i行有j列放了1个炮,k列放了2个炮的方案数。