poj1185-炮兵阵地(状态压缩dp)

炮兵阵地
Time Limit: 2000MS   Memory Limit: 65536K
Total Submissions: 25647   Accepted: 9892

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
  1 #include <map>
  2 
  3 #include <cstdio>
  4 
  5 #include <iostream>
  6 
  7 #include <string.h>
  8 
  9 #include <algorithm>
 10 
 11 using namespace std;
 12 
 13 const int maxn=100;
 14 
 15 int n,m,cnt,ans;
 16 
 17 char ch;
 18 
 19 int a[maxn+1];
 20 
 21 int vald[maxn+1][1<<10+1];
 22 
 23 int valdcount[maxn+1];
 24 
 25 int dp[maxn+1][maxn+10][maxn+10];
 26 
 27 //map<int,int> ma[maxn+10];
 28 
 29 bool judge(int x){
 30     return (!(x&(x<<1))&&!(x&(x<<2)));
 31 }
 32 
 33 int BitCount(unsigned int n)
 34 {
 35     unsigned int c =0 ;
 36     for (c =0; n; ++c)
 37     {
 38         n &= (n -1) ;
 39     }
 40     return c ;
 41 }
 42 
 43 int main()
 44 {
 45         scanf("%d%d",&n,&m);
 46         for(int i=1;i<=n;i++){
 47             for(int j=1;j<=m;j++){
 48                 scanf(" %c",&ch);
 49                 if(ch=='H'){
 50                     a[i]|=(1<<(j-1));
 51                 }
 52             }
 53         }
 54         for(int i=1;i<=n;i++){
 55             cnt=0;
 56             for(int j=0;j<=((1<<m)-1);j++){
 57                 if((!(j&a[i]))&&judge(j)){
 58                     vald[i][cnt]=j;
 59                     //ma[i][j]=cnt;
 60                     cnt++;
 61                 }
 62             }
 63             valdcount[i]=cnt-1;
 64         }
 65         for(int i=0;i<=valdcount[1];i++){
 66             int x=vald[1][i];
 67             for(int j=0;j<=(1<<(m)-1);j++) dp[1][i][j]=BitCount(x);
 68         }
 69         for(int i=0;i<=valdcount[2];i++){
 70             int x=vald[2][i];
 71             int q=BitCount(x);
 72             for(int j=0;j<=valdcount[1];j++) {
 73                 int y=vald[1][j];
 74                 if(!(x&y)) dp[2][i][j]=max(dp[2][i][j],q+BitCount(y));
 75             }
 76         }
 77         for(int i=3;i<=n;i++){
 78             for(int j=0;j<=valdcount[i];j++){
 79                    int x=vald[i][j];
 80                    int q=BitCount(x);
 81                   for(int k=0;k<=valdcount[i-1];k++){
 82                       int y=vald[i-1][k];
 83                     for(int h=0;h<=valdcount[i-2];h++){
 84                             int z=vald[i-2][h];
 85                            if(!(x&y)&&!(y&z)&&!(x&z)){
 86                                  dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][h]+q);
 87                            }
 88                     }
 89                   }
 90             }
 91         }
 92         for(int i=0;i<=valdcount[n];i++){
 93             int x=vald[n][i];
 94             for(int j=0;j<=valdcount[n-1];j++) {
 95                     ans=max(ans,dp[n][i][j]);
 96             }
 97         }
 98         printf("%d\n",ans);
 99     return 0;
100 }

 

posted @ 2016-08-07 14:47  GeniusYang  阅读(149)  评论(0编辑  收藏  举报