POJ 1185 状压DP

legal[] 保存所有在当前行可显示的状态，由dfs得到，len[]保存legal[]对应下标状态中的 1 的个数 ， 也就是放置炮台的个数

state[i] 表示第 i 行这块区域的土地情况，H表示 1 ，P表示 0

那么每次加入一个legal状态  都要符合 !(legal[i] & state[k])

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int dp[105][80][80],vis[12],state[105],legal[80],len[80];
int cnt,n,m;

void dfs(int u , int _state , int l)
{
    if(u > m){
        legal[++cnt] = _state;
        len[cnt] = l;
        return;
    }
    int flag = 1;
    if(u>1){
        if(vis[u-1]) flag=0;
    }
    if(u>2){
       if(vis[u-2]) flag=0;
    }
    if(flag){
        vis[u] = 1;
        dfs(u+1 , _state<<1|1 , l+1);
        vis[u] = 0;
    }
    dfs(u+1 , _state<<1 , l);
}

void DP()
{
    for(int k=1;k<=n;k++){
        for(int i=1;i<=cnt;i++){
            for(int j=1;j<=cnt;j++){
                if(k == 1 && !(state[k] & legal[i])){
                    dp[k][j][i] = len[i];
                }
                for(int t=1;t<=cnt;t
                ++){
                    if(!(state[k] & legal[i]) && !(legal[i] & legal[j]) && !(legal[i] & legal[t]) && !(legal[t] & legal[j])){
                        dp[k][j][i] = max(dp[k][j][i] , dp[k-1][t][j] + len[i]);
                    }
                }
            }
        }
    }
}

int main()
{
    //freopen("POJ1185.test","rb",stdin);
    char c;
    while(~scanf("%d%d",&n,&m)){
        for(int i=1;i<=n;i++){
            state[i] = 0;
            for(int j=1 ; j<=m ; j++){
                cin>>c;
                if(c == 'H') state[i] = state[i]<<1|1;
                else state[i] = state[i]<<1;
            }
        }
        cnt=0;
        memset(vis,0,sizeof(vis));
        dfs(1,0,0);

        memset(dp,0,sizeof(dp));
        DP();

        int ans=0;
        for(int i=1;i<=cnt;i++){
            for(int j=1;j<=cnt;j++){
                ans = max(dp[n][i][j] , ans);
            }
        }
        cout<<ans<<endl;
    }

    return 0;
}