poj1185 炮兵布阵 状压dp

题目链接：http://poj.org/problem?id=1185

题意： $n\ast m$ 的矩阵， 字符$P$表示空地（能放炮台）， 字符$H$表示山地（不能放炮台）,每个炮台的攻击范围是它上下左右两格之内的所有格子（见题图），问最多能摆几个炮台。

题解：状压dp。
首先分析，炮台的约束条件：上下左右两格内不能有炮台。
$H$表示 1，为不能放的点，$P$表示0，能放的点。将每一行的字符转换成二进制数，再压缩成一个十进制状态。
预处理出每一行的状态。$con[i]$数组表示第$i$行的状态。
然后再预处理对于某一行所有合法的状态中零的个数，$st[i]$表示状态为$i$，$w[i]$表示状态$i$的二进制表示中0的个数。
那么如何判断合法？ （见代码）
最后再预处理第一行的方案数。

$dp[i][j][k]$表示对于第$i$行，状态为$j$，前$i-1$行状态为$k$时的方案数
状态转移方程就是：
$dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][h]+w[j])$

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<stdio.h>
#include<string.h>
#include<queue>
#include<cmath>
#include<map>
#include<set>
#include<vector> 
using namespace std;
#define inf 0x3f3f3f3f
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define mem(a,b) memset(a,b,sizeof(a));
#define lowbit(x)  x&-x;  
#define debugint(name,x) printf("%s: %d\n",name,x);
#define debugstring(name,x) printf("%s: %s\n",name,x);
typedef long long ll;
typedef unsigned long long ull;
const double eps = 1e-6;
const int maxn = 1e5+5;
const int mod = 1e9+7;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m;
char s[105][105];
int a[105][105];
int w[205],st[205],con[205];
ll dp[105][205][205]; //第i行 第j种状态的方案数

int cal(int x){
    int res = 0;
    while(x){
        if(x&1) res++;
        x >>= 1;
    }
    return res;
}
int main() {
    scanf("%d%d",&n,&m);
    for(int i = 1; i <= n; i++){
        scanf("%s",s[i]+1);
        for(int j = 1; j <= m; j++){
            if(s[i][j] == 'H') a[i][j] = 1;
            else a[i][j] = 0;
        }
    }
    for(int i = 1; i <= n; i++){
        int tmp = 0;
        for(int j = 1; j <= m; j++){
            tmp = (tmp<<1) + a[i][j];
        }
        con[i] = tmp;
    }
    int sz = 1<<m;
    int tot = 0;
    for(int i = 0; i < sz; i++){
        if(i&(i<<2)) continue; //左移两位如果相与为1，表示冲突，不合法
        if(i&(i<<1)) continue; //同上
        st[++tot] = i;
        w[tot] = cal(i);
    }
    for(int i = 1; i <= tot; i++){
        if(con[1]&st[i]) continue;  //第 1 行的状态和 状态 i 有冲突，所以第1行的状态不能选 i 
        dp[1][i][1] = w[i];
    }
    for(int i = 2; i <= n; i++){  //当前行
        for(int j = 1; j <= tot; j++){  //状态为j
            if(st[j]&con[i]) continue;
            for(int k = 1; k <= tot; k++){ //前一行状态为k
                if(st[k]&con[i-1]) continue;
                if(st[j]&st[k]) continue;
                for(int h = 1; h <= tot; h++){ //前两行状态为h
                    if(st[h]&con[i-2]) continue;
                    if(st[h]&st[k]) continue;
                    if(st[h]&st[j]) continue;
                    dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][h]+w[j]);
                }
            }
        }
    }
    ll ans = 0;
    for(int i = 1; i <= tot; i++)
        for(int j = 1; j <= tot; j++)
            ans = max(dp[n][i][j],ans);
    printf("%lld\n",ans);
}