HYSBZ - 4806 炮 [dp]

题意：MxN的棋盘中最多放多少个炮，使得互相不能被打。
dp解决，M，N的范围为100， 不考虑状压。
状态：dp[i][j][k]表示当前行为i，前i-1行有j列没有炮，有k列有1个炮的方案数。
决策：
if(j-1>=0) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-1][k]*(m-(j-1)-k))%p; ///在没炮的地方放一个
if(j+1 <= m && k >= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+1][k-1]*(j+1))%p; ///在一个炮的地方放一个
if(j >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-2][k]*C(m-k-(j-2)))%p; ///在没炮的地方放两个
if(k >= 1 && j>= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j][k-1](m-j-(k-1))(j))%p; ///在没炮的地方和一个炮的地方分别放一个
if(j+2 <= m && k >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+2][k-2]*C(j+2))%p; ///在两个一个炮的地方分别放一个炮

边界条件：dp[0][0][0] = 1, else = 0

#include<iostream>
#include<string>
#include<cstdio>
#include<cstring>
#include<bitset>
#include<algorithm>
#include<map>
#include<set>
#include<queue>
#include<vector>
#include<cstdlib>
#include<list>
#include<stack>
#include<cmath>
using namespace std;
typedef long long LL;
const int p = 999983;
int n, m;
LL dp[105][105][105], ans;
LL C(LL x)
{
    return x*(x-1)/2;
}
int main()
{
    ios::sync_with_stdio(false);
    #ifndef ONLINE_JUDGE
    freopen("input.txt", "r", stdin);
    #endif // ONLINE_JUDGE
    while(cin >> n >> m)
    {
        memset(dp, 0, sizeof(dp));
        dp[0][0][0] = 1;
        for(int i = 1; i <= n; i++)
            for(int j = 0; j <= m; j++)
                for(int k = 0; j+k <= m; k++)
                {
                    dp[i][j][k] += dp[i-1][j][k] %= p;///不放
                    if(j-1>=0) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-1][k]*(m-(j-1)-k))%p; ///在没炮的地方放一个
                    if(j+1 <= m && k >= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+1][k-1]*(j+1))%p; ///在一个炮的地方放一个
                    if(j >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-2][k]*C(m-k-(j-2)))%p;    ///在没炮的地方放两个
                    if(k >= 1 && j>= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j][k-1]*(m-j-(k-1))*(j))%p;   ///在没炮的地方和一个炮的地方分别放一个
                    if(j+2 <= m && k >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+2][k-2]*C(j+2))%p;    ///在两个一个炮的地方分别放一个炮
                }

        LL ans = 0;
        for(int i = 0; i <= m; i++)
            for(int j = 0; j+i <= m; j++)
                ans = (ans+dp[n][i][j])%p;
        cout << ans << endl;
    }
    return 0;
}