【bzoj1801】【AHOI2009】【chess中国象棋】【组合数学】

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

题解：首先状压一定超时。

           考虑一下每行每列最多只能放2个。

           所以我们可以设f[i][j][k]为前i行放1个的有j列。放2个的有k列。

           转移的时候分类讨论即可。具体看代码吧。

代码：

#include<iostream>
#include<cstdio>
#define P 9999973
using namespace std;
int n,m,ans;
long long f[110][110][110];
int cal(int x){return (x*(x-1)/2)%P;} 
int main(){
   scanf("%d%d",&n,&m);
   f[0][0][0]=1;
   for (int i=1;i<=n;i++)
     for (int j=0;j<=m;j++)
       for (int k=0;k+j<=m;k++){
           f[i][j][k]+=f[i-1][j][k];//不放。 
           if (k>=1) (f[i][j][k]+=f[i-1][j+1][k-1]*(j+1))%=P;//在已经放一个的地方放。 
           if (j>=1) (f[i][j][k]+=f[i-1][j-1][k]*(m-k-j+1))%=P;//在没有放的地方放。 
           if (j>=2) (f[i][j][k]+=f[i-1][j-2][k]*cal(m-k-j+2))%=P;//在没有放的地方放两个
           if (k>=2) (f[i][j][k]+=f[i-1][j+2][k-2]*cal(j+2))%=P;//在放一个的地方放两个。
           if (k>=1)(f[i][j][k]+=f[i-1][j][k-1]*j*(m-j-k+1))%=P;//在放一个和没放的地方各放一个。 
         }
    for(int i=0;i<=m;i++)
      for (int j=0;j+i<=m;j++)
        (ans+=f[n][i][j])%=P;
    cout<<ans<<endl; 
}