K - 炮 bzoj- 4806 (DP)

众所周知，双炮叠叠将是中国象棋中很厉害的一招必杀技。炮吃子时必须隔一个棋子跳吃，即俗称"炮打隔子"。 

炮跟炮显然不能在一起打起来，于是rly一天借来了许多许多的炮在棋盘上摆了起来……他想知道，在N×M的矩形

方格中摆若干炮（可以不摆）使其互不吃到的情况下方案数有几种。

棋子都是相同的。

Input

一行，两个正整数N和M。

N<=100,M<=100

Output

一行，输出方案数mod 999983。

Sample Input

1 3

Sample Output

7

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <queue>
#define mem(p,k) memset(p,k,sizeof(p));
#define pb push_back
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define inf 0x3f3f3f3f
#define ll long long
using namespace std;
const int N=110;
const int mod=999983;
int n,m,T;
ll dp[N][N<<1][N<<1];
int main()
{
    while(~scanf("%d%d",&n,&m)){
        mem(dp,0);
        dp[0][m][0]=1;
        int ans=0;
        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];//no zhuan
                    if(k>=1)dp[i][j][k]+=dp[i-1][j+1][k-1]*(j+1);//one from 0
                    dp[i][j][k]+=dp[i-1][j][k+1]*(k+1);//one from 1
                    if(k>=2)dp[i][j][k]+=dp[i-1][j+2][k-2]*(j+1)*(j+2)/2;//two from 0
                    dp[i][j][k]+=dp[i-1][j][k+2]*(k+2)*(k+1)/2;//two from 1
                    dp[i][j][k]+=dp[i-1][j+1][k]*(j+1)*k;//two from 0 and 1
                    dp[i][j][k]%=mod;
                    if(i==n){
                        ans+=dp[i][j][k];
                        ans%=mod;
                    }
                    //cout<<j<<k<<" = "<<dp[i][j][k]<<endl;
                }
            }

        }
        cout<<ans<<endl;

    }


    return 0;
}
/*


*/