容易发现每行和每列炮数不超过2的状态是合法的
所以把行和列当成阶段都可以
f[i][j][k]
f
[
i
]
[
j
]
[
k
]
表示前i行有j列只有一个棋子,k列有两个棋子,这样就能算出没有棋子的列数
然后大力分类讨论
把行作为阶段的话,答案是
∑f[n][i][j]
∑
f
[
n
]
[
i
]
[
j
]
dp的时候忘记取模wa了几发,样例过了真的没有什么用
/**************************************************************
Problem: 1801
User: YJMSTR
Language: C++
Result: Accepted
Time:308 ms
Memory:8872 kb
****************************************************************/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define Max(_A,_B) (_A>_B?_A:_B)
#define Min(_A,_B) (_A<_B?_A:_B)
#define ll long long
int read(){
int s = 0;bool f = 1;char c = getchar();
while(c > '9' || c < '0'){if(c == '-') f = 0; c = getchar();}
while(c >= '0' && c <= '9') {
s = s * 10 + c - '0';
c = getchar();
}
return f == 1 ? s : -s;
}
const int md = 9999973, maxn = 101;
int n, m;
ll f[101][101][101];
/*
某一行和某一列中的炮的数量最多为2个
f[i][j][k] 表示前i行有j列只有一个棋子,k列有两个棋子
大力分类讨论:如果i + 1 行不放棋子 f[i+1][j][k] += f[i][j][k]
如果在i+1行只放了一个棋子的某一列放一个棋子 记该列编号为rua
f[i + 1][j - 1][k + 1] += f[i][j][k] * j
如果在i+1行的没放棋子的某一列放一个棋子
f[i + 1][j + 1][k] += f[i][j][k] *(m - j - k)
如果在i + 1行两个只放一个棋子的列放棋子 要求j >= 2
f[i + 1][j - 2][k + 2] += f[i][j][k] * C(j, 2)
如果在i + 1 行 一个只放了一个棋子的列 与一个未放棋子的列
f[i + 1][j][k + 1] += f[i][j][k] * C(j + m - k - j ,2), j >= 1 && j + k < m
如果在i + 1 行的两个空列放棋
f[i + 1][j + 2][k] += f[i][j][k] * C(m - k - j, 2);
边界: f[i][0][0] = 1;
*/
ll C(int x){//x里选2
return (ll)x*(x - 1) / 2;
}
int main() {
n = read(), m = read();
f[0][0][0] = 1;
for(int i = 0; i < n; i++) {
for(int j = 0; j <= m; j++) {
for(int k = 0; k + j <= m; k++) {//放了两个棋的列数
if(f[i][j][k]) {//转移不放的状态
f[i + 1][j][k] = (f[i + 1][j][k] + f[i][j][k]) % md;
}
if(j >= 1) f[i + 1][j - 1][k + 1] = (f[i + 1][j - 1][k + 1] + f[i][j][k] * j) % md;
if(m - j - k >= 1) f[i + 1][j + 1][k] = (f[i + 1][j + 1][k] + f[i][j][k] * (m - j - k)) % md;
if(j >= 2) f[i + 1][j - 2][k + 2] = (f[i + 1][j - 2][k + 2] + f[i][j][k] * C(j)) % md;
if(m - j - k >= 1 && j >= 1) f[i + 1][j][k + 1] = (f[i + 1][j][k + 1] + f[i][j][k] * (m - j - k) * j) % md;
if(m - j - k >= 2) f[i + 1][j + 2][k] = (f[i + 1][j + 2][k] + f[i][j][k] * C(m - j - k)) % md;
}
}
}
ll res = 0;
for(int i = 0; i <= m; i++)
for(int j = 0; j + i <= m; j++)
res = (res + f[n][i][j]) % md;
printf("%lld\n", res);
return 0;
}