bzoj1801 [Ahoi2009]chess 中国象棋

【bzoj1801】[Ahoi2009]chess 中国象棋

2014年6月7日1,2581

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

 

一行最多放两个，一列最多放两个。

f[i][j][k]表示前i行有j行有一个，k行有两个

然后转移显然。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;

const int mod=9999973,NN=107;

int n,m;
ll ans,f[NN][NN][NN];

inline ll C(ll x)
{
    return x*(x-1)/2;
}
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<=m-j;k++)
            {
                f[i][j][k]=f[i-1][j][k];                                    //什么都不放、 
                if(j)f[i][j][k]+=f[i-1][j-1][k]*(m-j-k+1);                    //在0个中再放一个。 
                if(j>1)f[i][j][k]+=f[i-1][j-2][k]*C(m-j-k+2);                //在0中放两个。 
                if(k&&j<m)f[i][j][k]+=f[i-1][j+1][k-1]*(j+1);                //在1个中放1个。 
                if(j&&k)f[i][j][k]+=f[i-1][j][k-1]*(m-j-k+1)*j;                //在0个和1个割放一个。 
                if(k>1&&j<m-1)f[i][j][k]+=f[i-1][j+2][k-2]*C(j+2);            //在1个中放两个。 
                f[i][j][k]%=mod;                                            //取模数。 
            }
    for(int i=0;i<=m;i++)
        for(int j=0;j<=m-i;j++)
            ans+=f[n][i][j];
    ans%=mod;
    
    printf("%lld",ans);
}