[AHOI2009]中国象棋

题目描述

这次小可可想解决的难题和中国象棋有关，在一个N行M列的棋盘上，让你放若干个炮（可以是0个），使得没有一个炮可以攻击到另一个炮，请问有多少种放置方法。大家肯定很清楚，在中国象棋中炮的行走方式是：一个炮攻击到另一个炮，当且仅当它们在同一行或同一列中，且它们之间恰好 有一个棋子。你也来和小可可一起锻炼一下思维吧！
输入格式

一行包含两个整数N，M，之间由一个空格隔开。
输出格式

总共的方案数，由于该值可能很大，只需给出方案数模9999973的结果。

一开始就想着用状压dp，可这题数据量有点大，压不了，看了题解后才搞懂，很新奇的一个dp状态，dp[i][j][k]:
i表示前i行，j表示有j列棋子个数为1，k表示有k列棋子个数为2。
这样可以保证每一列的合法性，然后每次向i+1行递推时，只放0，1，2个棋子，从而保证每一行的合法性。
递推式的话就分别考虑每种情况，例如：放一个棋子，放在有0个棋子的行，放在有1个棋子的行。把所有情况考虑即可。
过程中会爆int，用long long

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define Max 0x3f3f3f3f
typedef long long LL;
using namespace std;
const int maxn=105;

LL dp[maxn][maxn][maxn];//前i行，有j *列放了1个棋子，k *列放了2个棋子

int main()
{
    int n,m;
    LL mod=9999973;
    scanf("%d%d",&n,&m);

    /*
    dp[1][0][0] = 1;
    dp[1][1][0] = m;
    dp[1][2][0] = m * (m-1) / 2;
    */
    dp[0][0][0] = 1;
    for(int i=0; i<=n-1; i++){
        for(int j=0; j<=m; j++){
            for(int k=0; k+j<=m; k++){
                dp[i+1][j][k] = (dp[i+1][j][k] + dp[i][j][k])%mod;
                if(j+1+k<=m)
                    dp[i+1][j+1][k] = (dp[i+1][j+1][k] + dp[i][j][k] * (m-j-k))%mod;
                if(j+k+1<=m && j>=1)
                    dp[i+1][j][k+1] = (dp[i+1][j][k+1] + dp[i][j][k] * (m-j-k) * j)%mod;
                if(j+2+k<=m)
                    dp[i+1][j+2][k] = (dp[i+1][j+2][k] + dp[i][j][k] * ((m-j-k)*(m-j-k-1)/2))%mod;
                if(j>=1)
                    dp[i+1][j-1][k+1] = (dp[i+1][j-1][k+1] + dp[i][j][k] * j)%mod;
                if(j>=2)
                    dp[i+1][j-2][k+2] = (dp[i+1][j-2][k+2] + dp[i][j][k] * (j*(j-1)/2))%mod;
            }
        }
    }
    LL ans=0;
    for(int j=0; j<=m; j++){
        for(int k=0; k+j<=m; k++){
            ans = (ans + dp[n][j][k])%mod;
        }
    }
    printf("%lld\n",ans);
}