【bzoj4806】炮

Description

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

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

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

棋子都是相同的。

Input

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

N<=100,M<=100

Output

一行，输出方案数mod 999983。

Sample Input

1 3

Sample Output

7

HINT

Source

By FancyCoder

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int mod=999983;
LL f[110][110][110];//f[i][j][k]表示前i行j列有一个棋子k列有两个棋子 
int main()
{
	int n,m;
	scanf("%d%d",&n,&m);
	f[0][0][0]=1;
	//每行每列只有可能放2个棋子
	for(int i=1;i<=n;i++)//放第i行
	{
		for(int j=0;j<=m;j++)//j列有一个棋子
		{
			for(int k=0;k<=m-j;k++)//k列有两个棋子
			{
				//第i行不放
				f[i][j][k]=( f[i][j][k]+f[i-1][j][k] )%mod;
				//放1个在第i行,放在已有1个棋子的列上
				if(k>=1) f[i][j][k]=( f[i][j][k]+f[i-1][j+1][k-1]*(j+1) )%mod;
				//放1个在第i行,放在没有棋子的列上
				if(j>=1) f[i][j][k]=( f[i][j][k]+f[i-1][j-1][k]*(m-j-k+1) )%mod;
				//放2个在第i行,都各放在空的列上
				if(j>=2) f[i][j][k]=( f[i][j][k]+f[i-1][j-2][k]*( (m-j-k+2)*(m-j-k+1)/2 ) )%mod;
				//放2个在第i行,有棋子和没棋子的列各放一个
				if(j>=1 && k>=0) f[i][j][k]=( f[i][j][k]+f[i-1][j][k-1]*j*(m-j-k+1) )%mod;
				//放2个在第i行,都放在有一个棋子的列上
				if(k>=2) f[i][j][k]=( f[i][j][k]+f[i-1][j+2][k-2]*( (j+1)*(j+2)/2 ) )%mod;
			}
		}
	}
	LL ans=0;
	for(int j=0;j<=m;j++)//j列有一个棋子
	{
		for(int k=0;k<=m-j;k++)//k列有两个棋子
		{
			ans=(ans+f[n][j][k])%mod;
		}
	}
	printf("%lld\n",ans);
	return 0;
}