码迷,mamicode.com
首页 > 其他好文 > 详细

poj1185 状压dp

时间:2015-05-15 01:11:52      阅读:141      评论:0      收藏:0      [点我收藏+]

标签:

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;
}
View Code

 

poj1185 状压dp

标签:

原文地址:http://www.cnblogs.com/--560/p/4504897.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!