bzoj 1801 [Ahoi2009]chess 中国象棋

http://www.elijahqi.win/archives/3242
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

Source

Day2

分析题意 一行只可以有连个 一列只可以有两个

设 dp[i][j][k]表示当前处于第i行我 一列上有一个的有j个 一列上有两个的有k个

那么仔细分析就有六种转移直接写就okay

有一种转移 就是可以在空的地方放一个 可以在1的地方放一个就可以表现出凭空增加一个2的错觉

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int mod=9999973;
const int N=110;
inline void inc(int &x,ll v){x=(x+v)%mod;}
int n,m,dp[N][N][N];
int main(){
    freopen("bzoj1801.in","r",stdin);
    n=read();m=read();dp[0][0][0]=1;
    for (int i=0;i<n;++i){
        for (int j=0;j<=m;++j){
            for (int k=0;k<=m;++k){
                if (!dp[i][j][k]) continue;
                inc(dp[i+1][j][k],dp[i][j][k]);
                if (j>=1&&m-j-k>=1) inc(dp[i+1][j][k+1],dp[i][j][k]*(ll)j*(m-k-j)); 
                if (j+k+1<=m) inc(dp[i+1][j+1][k],dp[i][j][k]*(ll)(m-k-j));
                if (j+k+2<=m) inc(dp[i+1][j+2][k],dp[i][j][k]*(ll)(m-k-j)*(m-k-j-1)>>1);
                if(j>1) inc(dp[i+1][j-2][k+2],(ll)dp[i][j][k]*j*(j-1)>>1);
                if(j>=1&&k+1<=m) inc(dp[i+1][j-1][k+1],(ll)dp[i][j][k]*j);
            }
        }
    }
    int ans=0;
    for (int j=0;j<=m;++j)
        for (int k=0;k<=m-j;++k) inc(ans,dp[n][j][k]);
    printf("%d\n",ans);
    return 0;
}