BZOJ1801 Ahoi2009 chess 中国象棋

1801: [Ahoi2009]chess 中国象棋

Time Limit: 10 Sec   Memory Limit: 64 MB
Submit: 1911   Solved: 1105
[ Submit][ Status][ Discuss]

Description

在N行M列的棋盘上，放若干个炮可以是0个，使得没有任何一个炮可以攻击另一个炮。 请问有多少种放置方法，中国像棋中炮的行走方式大家应该很清楚吧.

Input

一行包含两个整数N，M，中间用空格分开.

Output

输出所有的方案数，由于值比较大，输出其mod 9999973

Sample Input

1 3

Sample Output

7

HINT

除了在3个格子中都放满炮的的情况外，其它的都可以.

100%的数据中N,M不超过100
50%的数据中，N,M至少有一个数不超过8
30%的数据中，N,M均不超过6

【题解】题目大意，n*m的棋盘中放棋子，每行每列不得超过2个

      dp即可，好像可以滚一下？没必要。。懒。。

#include <bits/stdc++.h>
using namespace std;
#define maxn 105
#define p 9999973
#define ll long long
ll dp[105][105][105];
ll C(int n)
{
	return (n*(n-1)/2);
}
int main()
{
	int n,m;
	scanf("%d%d",&n,&m);
	dp[0][0][0]=1;
	for (int i=0;i<n;i++)
		for (int j=0;j<=m;j++)
			for (int k=0;j+k<=m;k++)
			{
				 dp[i+1][j][k]=(dp[i+1][j][k]+dp[i][j][k])%p;
				 if (m-j-k>0)dp[i+1][j+1][k]=(dp[i+1][j+1][k]+dp[i][j][k]*(m-j-k))%p;
				 if (m-k-j>1)dp[i+1][j+2][k]=(dp[i+1][j+2][k]+dp[i][j][k]*C(m-j-k))%p;
				 if (j>0)dp[i+1][j-1][k+1]=(dp[i+1][j-1][k+1]+dp[i][j][k]*j)%p;
				 if (j>0 && m-j-k>0)dp[i+1][j][k+1]=(dp[i+1][j][k+1]+dp[i][j][k]*j*(m-j-k))%p;
				 if (j>1)dp[i+1][j-2][k+2]=(dp[i+1][j-2][k+2]+dp[i][j][k]*C(j))%p;
			}
	ll ans=0;
	for (int i=0;i<=m;i++)
		 for (int j=0;(i+j)<=m;j++)
			 (ans+=dp[n][i][j])%=p;
	cout<<ans<<endl;
	return 0;
}