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。
接下来的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 }