bzoj1801[AHOI2009]CHESS中国象棋

题意：在棋盘上放一些炮使得它们不互相攻击。其实就是一行/一列最多放两个。

50分的数据中n，m至少有一个不超过8，比较直接的想法是对n/m中较小的一维做状态压缩,状态f[i][S1][S2]表示在前i行/列中,S1集合中的列/行放了1炮，S2集合中的列/行放了2炮。转移的时候，需要枚举第i行/列怎么放炮。这个时间复杂度好像炸天了….

仔细观察转移的过程,会发现S1和S2集合中的元素到底是谁其实无关紧要，我们只需要知道S1和S2的元素数目即可完成转移。所以可以把状态改成f[i][j][k]表示前i行有j列放了1炮,k列放了2炮的方案数，依然讨论第i行的方法，不过不需要逐个枚举方案。

因为闲得蛋疼把数组开成了两维，我的DP顺序比较奇怪…某个状态一定是从总炮数小于等于它的状态转移过来(等于的情况只有不放炮一种，写的时候不需要考虑，直接从之前继承下来),所以只要按照总炮数从大到小DP就行了。不过这题并没有卡内存,所以开成三维数组的话总是从i-1那一层转移过来,循环顺序怎样都没关系.

转移的时候,f[i][j][k]可以从f[i-1][j][k-1]转移过来,看起来是放一个炮的列数不变，放两个炮的列数+1，但不能理解成在一个空列上放两炮使得它变成有两炮的列,实际的过程是在一个空列放了一炮，在另一个之前有1炮的列上再放了1炮，因为在第i行每一列只能放一炮

我加了滚动数组之后864kb,404ms,bzoj上有人176kb就过了，而且跑得飞起,30ms.甚至还有0msAC的…不知道别人怎么做的，只能ym。

#include<cstdio>
const int maxn=105,mod=9999973;
int f[maxn][maxn];
int main(){
    int n,m;scanf("%d%d",&n,&m);
    f[0][0]=1;
    int lim=2*m;
    int ans=0;
    for(int i=1;i<=n;++i){
        for(int tot=lim;tot>=0;--tot){//倒着枚举，类似于01背包的方法
            for(int k=0;k*2<=tot;++k){
                int j=tot-k*2;
                if(j>m)continue;

　　　　　　　　　 //这一行放一炮
                if(j!=0){
                    f[j][k]+=(m-(k)-(j-1))*1LL*f[j-1][k]%mod;
                }

                if(k!=0){

                    f[j][k]+=(j+1)*1LL*f[j+1][k-1]%mod;

                }

                //这一行放两炮
                if(j>=2){
                    f[j][k]+=(m-k-(j-2))*1LL*(m-k-(j-1))/2LL*f[j-2][k]%mod;
                }
                if(k>=2){
                    f[j][k]+=(j+2)*(j+1)/2LL*f[j+2][k-2]%mod;
                }
                if(k>=1){
                    f[j][k]+=(m-j-k+1)*1LL*(j)*f[j][k-1]%mod;
                }
                f[j][k]%=mod;
                if(i==n)ans=(ans+f[j][k])%mod;
            }
        }
    }
    printf("%d\n",ans);
    return 0;
}