洛谷 P2051 [AHOI2009]中国象棋【dp】

---

题目：

洛谷传送门

---

题意：

有一张$n\ast m$的图，在这其中放置若干象棋炮兵，在它们不会互相攻击的情况下，求有多少合法方案

---

分析：

因为题目中已经明确告诉我们是象棋炮兵，所以每一行和每一列都至多放置$2$个炮兵
针对炮兵每行每列放置的数量极其有限，我们想出了方程$f_{i,j,k}$，表示已经放好了前$i$行，其中有$j$列放了一个炮兵，有$k$列放了两个炮兵
随着行数的递增，我们发现会有$5$中可能的情况出现：
$1.$新的一行什么也不做：
$$f_{i,j,k}+=f_{i-1,j,k}$$
$2.$新的一行只放一个炮兵
$(1)$我们在没有炮兵的一列放上，也就是多了一列放了一个炮兵的列：
$$f_{i,j+1,k}+=f_{i-1,j,k}\ast (m-j-k)(m-j-k>0)$$
$(2)$我们在只有一个炮兵的列放上，也就是少了一列只放一个的，多了一个放两个的列：
$$f_{i,j-1,k+1}+=f_{i-1,j,k}\ast j(j>0)$$
$3.$新的一行放两个炮兵
$(1)$两个都放在空列上，即多了两列只放一个的：
$$f_{i,j+2,k}+=f_{i-1,j,k}\ast C_{m-j-k}^2(m-j-k>1)$$
$(2)$两个都放在放了一个的上面，即少了两列只放一个的，多了两列放两个的：
$$f_{i,j-2,k+2}+=f_{i-1,j,k}\ast C_j^2(j>1)$$
$(3)$一个放在空列上，一个放在只有一个的上面：
$$f_{i,j,k+1}+=f_{i-1,j,k}\ast j\ast (m-j-k)(j>0andm-j-k>0)$$

---

代码：

// luogu-judger-enable-o2
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define LL long long 
#define LZX 9999973
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
LL f[105][105][105];
LL C(LL x) {return x*(x-1)/2;}
int main()
{
    LL n=read(),m=read();
    f[0][0][0]=1;
    for(LL i=0;i<n;i++)
      for(LL j=0;j<=m;j++)
        for(LL k=0;k+j<=m;k++)
        {
        	(f[i+1][j][k]+=f[i][j][k])%=LZX;
            if(m-j-k) (f[i+1][j+1][k]+=f[i][j][k]*(m-j-k))%=LZX;
            if(j) (f[i+1][j-1][k+1]+=f[i][j][k]*j)%=LZX;
            if(m-j-k&&j) (f[i+1][j][k+1]+=f[i][j][k]*j*(m-j-k))%=LZX;
            if(m-j-k>1) (f[i+1][j+2][k]+=f[i][j][k]*C(m-j-k))%=LZX;
            if(j>1) (f[i+1][j-2][k+2]+=f[i][j][k]*C(j))%=LZX;
        }
    LL ans=0;
    for(LL i=0;i<=m;i++)
      for(LL j=0;j+i<=m;j++)
        (ans+=f[n][i][j])%=LZX;
    cout<<ans%LZX;
    return 0;
}