洛谷P2051 [AHOI2009]中国象棋（状压DP）

题目链接：https://www.luogu.com.cn/problem/P2051

洛谷P2051 [AHOI2009]中国象棋

状压DP
想明白了其实思路挺简单的一道DP。
转换一下题目大意，为在n*m的棋盘上放置放置若干个棋子，使得每行每列的棋子数目都不超过2，求所有放置的方案数。
题目中最关键的限制条件->每行每列的棋子数目都不超过2。
若按照行的顺序依次放置棋子，可以得知当前这一行棋子的放置会受到前面棋子放置方案的影响，并且后面放置的行不会影响到前面行的种类数，可以考虑将其转换成DP来求解，每行为一个状态，并且在放置时需要记录前面过程的放置方案即状态。
由于每一行最多只能放置2个棋子，考虑用数组f[i][j][k]表示前i行棋盘中，有j列放置了一个棋子，有k列放置了2个棋子的方案数，这样就完成了前i行放置方案状态的表示。
考虑将状态由第i-1行转移至第i行。
第i行有三种放置方案：
①不放置
②放置1个
③放置2个
若不放置
则f[i][j][k]+=f[i-1][j][k]
若不放置1个
则这一个可以放置在没有棋子的列中即f[i][j][k]+=f[i-1][j-1][k]*(m-j-k+1)
也可以放置在有一个棋子的列中即f[i][j][k]+=f[i-1][j+1][k-1]*(j+1)
若不放置2个
则这两个可以都放置在没有棋子的列中即f[i][j][k]+=f[i-1][j-2][k]*c(m-j-k+2,2)
也可以都放置在有一个棋子的列中即f[i][j][k]+=f[i-1][j+2][k-2]*c(j+2,2)
也可以一个放置在没有棋子的列中，一个放置在有棋子的列中即f[i][j][k]+=f[i-1][j][k-1]*j*(m-j-k+1)
计算过程中注意取模，最后结果即为${\sum }_{j=0}^m$${\sum }_{k=0}^{m-j}f[n][j][k]$。

#include<bits/stdc++.h>
#define next next_
#define y1 yy
#define hash hash_
#define size size_
#define complex complex_
using namespace std;

using ll=long long;
using ull=unsigned long long;
using pii=pair<int,int>;
using pll=pair<ll,ll>;

const ll mod=9999973;
const double pi=acos(-1.0);
int _;

ll n,m,f[110][110][110],ans;
ll fac[110],inv[110];

ll c2(ll n){
	return n*(n-1)/2;
}

void work(){
	scanf("%lld%lld",&n,&m);
	f[0][0][0]=1;
	for(ll i=1;i<=n;i++){
		for(ll j=0;j<=m;j++)
			for(ll k=0;k<=m-j;k++){
				f[i][j][k]+=f[i-1][j][k];
				f[i][j][k]%=mod;
				if(j>=1) f[i][j][k]+=f[i-1][j-1][k]*(m-j+1-k)%mod;
				f[i][j][k]%=mod;
				if(k>=1) f[i][j][k]+=f[i-1][j+1][k-1]*(j+1)%mod;
				f[i][j][k]%=mod;
				if(j>=2) f[i][j][k]+=f[i-1][j-2][k]*c2(m-j-k+2)%mod;
				f[i][j][k]%=mod;
				if(k>=2) f[i][j][k]+=f[i-1][j+2][k-2]*c2(j+2)%mod;
				f[i][j][k]%=mod;
				if(k>=1) f[i][j][k]+=f[i-1][j][k-1]*j*(m-j-k+1)%mod;
				f[i][j][k]%=mod;
			}
	}
	for(ll j=0;j<=m;j++)
		for(ll k=0;k<=m-j;k++) ans+=f[n][j][k],ans%=mod;
	printf("%lld",ans);
}

int main(){
//	scanf("%d",&_);
	_=1;
	while(_--){
		work();
	}
	return 0;
}