[NOI2001]炮兵阵地

题目大意：
　　一个n*m(n<=100,m<=10)的格子图，里面有些地方可以放炮兵，有些地方不行，炮兵能像上下左右攻击到两格远的位置。
　　问在炮兵不会互相攻击的情况下，最多能放多少炮兵？

思路：
　　状压DP。
　　f[i][j][k]表示DP到第i行，当前行状态为j，上一行状态为k。
　　然后枚举当前行，这一行状态，上一行状态和上两行状态，判断是否合法并转移即可。
　　复杂度很悬，然而洛谷和百练上随便A，但是POJ确实过不去。
　　考虑横向互相攻击的情况，这对每一行都是一样的，而且似乎无效状态会有很多。
　　那我们就预处理出同一行中所有的合法情况，极限情况下不会超过60个。
　　这样枚举的时候就能节约不少时间。

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=100,M=10;
int a[N];
int f[2][1<<M][1<<M];
int sit[61];
int main() {
    int n=getint(),m=getint();
    for(register int i=0;i<n;i++) {
        for(register int j=0;j<m;j++) {
            char c[2];
            scanf("%1s",c);
            if(*c=='H') a[i]|=1<<j;
        }
    }
    for(register int i=0;i<(1<<m);i++) {
        for(register int t=2;t<m;t++) {
            if(((i>>t)&(i>>(t-1)))||((i>>(t-1))&(i>>(t-2)))||((i>>t)&(i>>(t-2)))) goto Next;
        }
        sit[++sit[0]]=i;
        Next:;
    }
    for(register int j=1;j<=sit[0];j++) {
        if(sit[j]&a[0]) continue;
        f[0][sit[j]][0]=__builtin_popcount(sit[j]);
    }
    for(register int j=1;j<=sit[0];j++) {
        if(sit[j]&a[1]) continue;
        for(register int k=1;k<=sit[0];k++) {
            if(sit[k]&(sit[j]|a[0])) continue;
            f[1][sit[j]][sit[k]]=std::max(f[1][sit[j]][sit[k]],f[0][sit[k]][0]+__builtin_popcount(sit[j]));
        }
    }
    for(register int i=2;i<n;i++) {
        memset(f[i&1],0,sizeof f[i&1]);
        for(register int j=1;j<=sit[0];j++) {
            if(sit[j]&a[i]) continue;
            for(register int k=1;k<=sit[0];k++) {
                if(sit[k]&(sit[j]|a[i-1])) continue;
                for(register int l=1;l<=sit[0];l++) {
                    if(sit[l]&(sit[j]|sit[k]|a[i-2])) continue;
                    f[i&1][sit[j]][sit[k]]=std::max(f[i&1][sit[j]][sit[k]],f[!(i&1)][sit[k]][sit[l]]+__builtin_popcount(sit[j]));
                }
            }
        }
    }
    int ans=0;
    for(register int i=1;i<=sit[0];i++) {
        for(register int j=1;j<=sit[0];j++) {
            ans=std::max(ans,f[!(n&1)][sit[i]][sit[j]]);
        }
    }
    printf("%d\n",ans);
    return 0;
}