poj1185 炮兵阵地 状压dp

题目链接:

https://vjudge.net/contest/159644#problem/B

题意:

在一个N*M的矩阵上布置炮兵部队，只有平原可以布置，然后每个炮兵部队都有一个攻击范围，它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。
问最多可以安排几个炮兵

题解:

http://www.cnblogs.com/buptLizer/archive/2012/08/23/2651881.html
状压dp
dp[i][j][k]:=考虑到第i行， 这一行是第i个合法状态，上一行是第j个合法状态，最多可以安排的炮兵
转移：dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][p]+num[j]); num[j]表示j这个合法状态可以安排的炮兵数

代码:

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
typedef long long ll;
#define MS(a) memset(a,0,sizeof(a))
#define MP make_pair
#define PB push_back
const int INF = 0x3f3f3f3f;
const ll INFLL = 0x3f3f3f3f3f3f3f3fLL;
inline ll read(){
    ll 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;
}
//////////////////////////////////////////////////////////////////////////
const int maxn = 80+10;

int dp[105][maxn][maxn];
int state[maxn],num[maxn],cnt;
int cur[105];
char mp[105][20];
int n,m;

int getnum(int x){
    int res = 0;
    while(x){
        res += (x&1);
        x = x>>1;
    }
    return res;
}

void init(){
    cnt = 0;
    for(int i=0; i<(1<<m); i++){
        if(!(i&(i<<1)) && !(i&(i<<2))){
            state[cnt] = i;
            num[cnt] = getnum(i);
            cnt++;
        }
    }
}

int main(){
    cin >> n >> m;
    for(int i=0; i<n; i++){
        cin >> mp[i];
        for(int j=0; j<m; j++){
            if(mp[i][j] == 'P')
                cur[i] += (1<<(m-j-1));
        }
    }

    memset(dp,-1,sizeof(dp));
    init();
    for(int i=0; i<cnt; i++)
        if(!(state[i]&(~cur[0]))) // 如果与能放炮兵的反有交集 就不可行
            dp[0][i][0] = num[i];
    for(int i=1; i<n; i++){
        for(int j=0; j<cnt; j++){
            if((~cur[i])&state[j]) continue;
            for(int k=0; k<cnt; k++){
                if((~cur[i-1])&state[k]) continue;
                if(state[j]&state[k]) continue;
                for(int p=0; p<cnt; p++){
                    if(i>=2 && ((~cur[i-2])&state[p])) continue;
                    if(state[j]&state[p]) continue;
                    if(dp[i-1][k][p] == -1) continue;
                    dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][p]+num[j]);
                }
            }
        }
    }

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

    return 0;
}

// http://www.cnblogs.com/buptLizer/archive/2012/08/23/2651881.html