poj1185 状压dp

poj1185   状压dp

炮兵阵地

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions: 20940		Accepted: 8104

Description

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示： 

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。 
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。 

Input

第一行包含两个由空格分割开的正整数，分别表示N和M； 
接下来的N行，每一行含有连续的M个字符('P'或者'H')，中间没有空格。按顺序表示地图中每一行的数据。N <= 100；M <= 10。

Output

仅一行，包含一个整数K，表示最多能摆放的炮兵部队的数量。

Sample Input

5 4
PHPP
PPHH
PPPP
PHPP
PHHP

Sample Output

6

思路：状压dp。按行进行递推，每一行的状态只与上一行和上上一行有关，而每一行的结果只是由上一行推出，满足无后效性，可dp，将每一行的状态用二进制整数存，筛掉明显不符合条件的，进行dp。
　　　　dp(i,j,k)表示第i行的决策做出后，第i行的状态为j，第i-1行的状态为k的炮兵数。
　　　　　　dp(i,j,k)=max(dp(i,j,k),dp(i-1,k,t)+sum[j]) (i,j,k状态存在且不互相冲突)
　　　　注意初始化筛掉不存在或不符合条件的状态

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>

using namespace std;

const int maxn=2200;
const int INF=(1<<29);

int N,M;
char s[maxn];
int Map[maxn];
int sum[maxn*300];
int dp[2][maxn][maxn];
int state[maxn*300],cnt;

void InitSum()
{
    memset(sum,0,sizeof(sum));
    for(int i=0;i<cnt;i++){
        for(int j=0;j<M;j++){
            if(state[i]&(1<<j)) sum[i]++;
        }
    }
}

bool jud(int n)
{
    if((n&(n<<1))||(n&(n<<2))) return 0;
    return 1;
}

void Init()
{
    cnt=0;
    for(int i=0;i<(1<<M);i++){
        if(jud(i)) state[cnt++]=i;
    }
}

int main()
{
    cin>>N>>M;
    memset(Map,0,sizeof(Map));
    for(int i=1;i<=N;i++){
        scanf("%s",s);
        for(int j=0;j<strlen(s);j++){
            if(s[j]=='H') Map[i]|=(1<<j);
        }
    }
    Init();
    InitSum();
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=N;i++){
        for(int j=0;j<cnt;j++){
            if(state[j]&Map[i]) continue;
            for(int k=0;k<cnt;k++){
                if(i==1) dp[i%2][j][k]=sum[j];
                else if((state[k]&Map[i-1])||(state[j]&state[k])) continue;
                else if(i==2) dp[i%2][j][k]=max(dp[i%2][j][k],dp[(i+1)%2][k][0]+sum[j]);
                else{
                    for(int t=0;t<cnt;t++){
                        if(!(state[t]&Map[i-2])&&!(state[j]&state[t])&&!(state[k]&state[t])) dp[i%2][j][k]=max(dp[i%2][j][k],dp[(i+1)%2][k][t]+sum[j]);
                    }
                }
            }
        }
    }
    int ans=0;
    for(int j=0;j<cnt;j++){
        for(int k=0;k<cnt;k++){
            if(dp[N%2][j][k]>ans) ans=dp[N%2][j][k];
        }
    }
    cout<<ans<<endl;
    return 0;
}