bzoj 1801（递推）

传送门
题解：状压dp只能得部分分，所以（有大佬）考虑dp[i][j][k]表示前i行有j列放了一个炮，k列放了两个炮。（这样刚好可以表示所有状态orzorz）然后分6种情况递推，最后把前n行所有合法情况加起来即可。（0表示列上没有炮，1表示有1个炮）
P.S.一开始死活弄不懂为什么dp[i][j][k]可以从dp[i-1][j+1][k-1]转移（为什么要+1？？？），后来发现”+1”就是此时放置一个新的炮得那一列，放完后，一个炮的列数就会-1，从而j+1变成当前枚举的j！
又一次智商被碾压orz( ＿ ＿)ノ｜扶墙

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll MOD=9999973;
int n,m;
ll dp[102][102][102],ans=0;
inline ll C2(int x) {
    return x*(x-1)>>1;
}
int main() {
    scanf("%d%d",&n,&m);
    dp[0][0][0]=1;
    for (int i=1;i<=n;++i)
        for (int j=0;j<=m;++j)
            for (int k=0;k<=m-j;++k) {
                dp[i][j][k]=dp[i-1][j][k];//不放
                if (j)          dp[i][j][k]=(dp[i][j][k]+dp[i-1][j-1][k]*(m-j-k+1))%MOD;//放在0上
                if (k&&j<m)     dp[i][j][k]=(dp[i][j][k]+dp[i-1][j+1][k-1]*(j+1))%MOD;//放在1上
                if (j>1)        dp[i][j][k]=(dp[i][j][k]+dp[i-1][j-2][k]*C2(m-j-k+2))%MOD;//放在0 0上
                if (k)          dp[i][j][k]=(dp[i][j][k]+dp[i-1][j][k-1]*(m-j-k+1)*j)%MOD;//放在0 1上
                if (k>1&&j<m-1) dp[i][j][k]=(dp[i][j][k]+dp[i-1][j+2][k-2]*C2(j+2))%MOD;//放在1 1上
            }
    for (int j=0;j<=m;++j)
        for (int k=0;k<=m-j;++k)
            ans+=dp[n][j][k],ans%=MOD;
    printf("%lld\n",ans);
    return 0;
}