BZOJ1801: [Ahoi2009]chess 中国象棋

容易发现每行和每列炮数不超过2的状态是合法的
所以把行和列当成阶段都可以
$f[i][j][k]$表示前i行有j列只有一个棋子，k列有两个棋子，这样就能算出没有棋子的列数
然后大力分类讨论
把行作为阶段的话，答案是$\sum{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;
}