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,花了很长时间终于做出来了(笑),一开始看错题目,还以为求总的方案数= =。

我的做法是先把图每一行都压缩成一个十进制数,'P'变成0,'Q'变成'1',这么写是为了之后判断能不能放炮兵提供方便(因为这样表示如果用&运算符算出来是1的只有放炮兵(1)和'Q'(1)这种情况,其他情况下运算后都是0,可以动手画一下),然后枚举所有每一行所有可能的摆放情况(j从0~(1<<m)-1),先判断是不是两格以内都是1,用if( (j&(j<<1)) || (j&(j<<2)))来判断,然后再判断当前的j和temp是不是有某位同时为1的情况,用if(j&temp)判断。用num1[i][j]来表示第i行第j种状态所压缩后的十进制数,用num1[i]来表示第i行所有状态的总个数,用soldier[i][j]来表示第i行第j种状态的该行炮兵总数,用dp[i][j][k]来表示第i行第j中状态下第i-1行第k中状态下的最大炮兵数。

先对第一行和第二行进行初始化,然后从第三行开始依次进行动规,转移方程是dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][h]+soldier[i][j]);

最后求值的时候要注意,如果n==1或者n==2要单独讨论,如果n>3的时候要注意最大值不一定是dp[n][j][k],因为可能最大值不是在最后一行取到的,所以要把ans的取值放在中间。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<vector>
#include<map>
#include<queue>
#include<stack>
#include<string>
#include<algorithm>
using namespace std;
int num1[120],num[120][70],dp[120][70][70],n,m,soldier[120][70];void hang(int i,int temp)
{int t=0,j,k,sum;for(j=0;j<(1<<m);j++){if( (j&(j<<1)) || (j&(j<<2)))continue; // || ((j<<1)&(j<<2))||(j&(j>>1)) || (j&(j>>2)) || ((j>>1)&(j>>2))  if(j&temp)continue;t++;num[i][t]=j;k=j;sum=0;while(k>0){if(k%2==1)sum++;k/=2;}soldier[i][t]=sum;}num1[i]=t;
}int main()
{int i,j,temp,k,h,sum,sum1,ans;char s[20];while(scanf("%d%d",&n,&m)!=EOF){memset(num1,0,sizeof(num1));memset(num,0,sizeof(num));memset(dp,0,sizeof(dp));memset(soldier,0,sizeof(soldier));for(i=1;i<=n;i++){scanf("%s",s);temp=0;for(j=0;j<m;j++){if(s[j]=='P'){temp=temp*2;}else if(s[j]=='H'){temp=temp*2+1;}}hang(i,temp);}ans=0;for(j=1;j<=num1[1];j++){dp[1][j][0]=soldier[1][j];}if(n==1){for(j=1;j<=num1[1];j++){ans=max(ans,dp[1][j][0]);}printf("%d\n",ans);continue;}ans=0;for(j=1;j<=num1[2];j++){for(k=1;k<=num1[1];k++){if(num[2][j]&num[1][k])continue;dp[2][j][k]=soldier[2][j]+soldier[1][k];ans=max(ans,dp[2][j][k]);}}if(n==2){printf("%d\n",ans);continue;}ans=0;for(i=3;i<=n;i++){for(j=1;j<=num1[i];j++){for(k=1;k<=num1[i-1];k++){if(num[i][j]&num[i-1][k])continue;for(h=1;h<=num1[i-2];h++){if( (num[i][j]&num[i-2][h]) || (num[i-1][k]&num[i-2][h])     )continue;dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][h]+soldier[i][j]);ans=max(ans,dp[i][j][k]);}}}}printf("%d\n",ans);}return 0;
}

poj1185 炮兵阵地相关推荐

  1. POJ1185 炮兵阵地 状压DP

    炮兵阵地 Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 31819   Accepted: 12295 Descriptio ...

  2. poj1185炮兵阵地【状压dp】

    Language: Default 炮兵阵地 Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 22142   Accepted ...

  3. POJ1185炮兵阵地(状压DP)

    炮兵阵地 Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 33176   Accepted: 12787 Descriptio ...

  4. poj1185炮兵阵地

    #include <iostream> #include <cstdio> #include <cmath> #include <algorithm> ...

  5. POJ-1185 炮兵阵地 动态规划+状态压缩

    由于递推的时候依赖于三个连续层的关系.一开始想着直接三重for循环,但是这里有个问题就是上一层的0位置上包括着上上层是0和1两种可能,而后者又对当前行有约束,因此该方法不行.当然有一个办法就是增加状态 ...

  6. POJ1185炮兵阵地(状态压缩 + dp)

    题目链接 题意:给出一张n * m的地图,其中 有的地方能放大炮,有的地方不能,大炮与上下左右两个单位范围内会相互攻击,问最多能放几个大炮 能放大炮为1不能放大炮为0,把每一行看做一个状态,要除去同一 ...

  7. hdu4539 郑厂长系列故事——排兵布阵 + POJ1158 炮兵阵地

    题意: 郑厂长系列故事--排兵布阵 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Other ...

  8. POJ 1185 炮兵阵地 (状压DP)

    炮兵阵地 Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 14869   Accepted: 5575 Description ...

  9. poj 1185 NYOJ 85 炮兵阵地(状态压缩dp)

    炮兵阵地 时间限制:2000 ms  |  内存限制:65535 KB 难度:6 描述 司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队.一个N*M的地图由N行M列组成,地图的每一格可能是山地 ...

最新文章

  1. 【linux】Valgrind工具集详解(十三):Helgrind(线程错误检测器)
  2. 基于SSM实现考研信息管理平台系统
  3. python数据分析numpy_利用python进行数据分析-NumPy高级应用
  4. Timer TimeTask Handler
  5. 用java和汇编开发一个Hello World系统内核
  6. OpenCV学习笔记之Mat
  7. ES9的新特性:异步遍历Async iteration
  8. Java打印三角形(双层for循环)
  9. matlab 类 继承,继承MATLAB中的密封类
  10. 基于Spring MVC的ECharts动态数据实时展示
  11. python调用nmap扫描全端口_python nmap实现端口扫描器教程
  12. C#实现微信扫码登录
  13. 使用前端方式挑战 Chrome 小恐龙游戏高分
  14. 谷歌学术搜索技巧,命令搜索
  15. 哪里可以在线观看《人再囧途之泰囧》
  16. 【移动安全高级篇】————1、Android沙盘原理与实现
  17. Windows下将文件夹下所有的文件名导出成生成TXT文件
  18. HKEY_LOCAL_MACHINE\Software\WOW6432Node
  19. 无法启动此程序因为计算机丢失gdiplus,处理系统提示无法启动此程序,因为计算机中丢失gdiplus.dll的方法...
  20. mongodb的NUMA问题

热门文章

  1. 在windows上设置代理服务器
  2. mysql根据表结构自动生成表设计文档
  3. php中警告提示A session had already been started – ignoring session_start() 解决方法
  4. 《The swirlds hashgraph consensus algorithm: Fair, fast, byzantine fault tolerance》Hashgraph论文的学习
  5. ConcurrentHashMap扩容原理
  6. 【Excel】分列操作,对一个单元格内的数据进行特定分列,如按照逗号,冒号,空格等等进行文本划分
  7. 计算机审计概念有统一标准,计算机审计的研究回顾.doc
  8. 文献笔记02 网络入侵检测技术综述(信息安全学报)
  9. 虚幻4漏光问题解决方法
  10. Nginx配置、原理学习