POJ1185 - 炮兵阵地(状态压缩DP)
题目大意
中文的。。直接搬过来。。。
司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用"P"表示),如下图。在每一格平原地形上最多可以布置一支炮兵部队(山地上不能够部署炮兵部队);一支炮兵部队在地图上的攻击范围如图中黑色区域所示:
如果在地图中的灰色所标识的平原上部署一支炮兵部队,则图中的黑色的网格表示它能够攻击到的区域:沿横向左右各两格,沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在,将军们规划如何部署炮兵部队,在防止误伤的前提下(保证任何两支炮兵部队之间不能互相攻击,即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内),在整个地图区域内最多能够摆放多少我军的炮兵部队。
题解
经典的状态压缩DP啦~~~~总算把这题过了~~~~~做完POJ3254再搞这题还是挺容易的~~~
对于当前第i行的放置情况,会影响到i-1以及i-2行的放置情况,我们可以先预处理一行的所有合法状态(用数组s表示)并记录此状态放置的大炮个数(数组sum表示),那怎么判断合法呢?很简单,对于状态s,只需要把s向右移一位即可,即(s>>1),然后判断一下s&(s>>1)是否为0,如果是,则不存在相邻的两个大炮,同理,还需要判断一下s&(s>>2)是否为0,即不存在距离为2的两个大炮。接下来就是状态转移了,由于会受前面两行的影响,所以状态转移方程得是三维的:dp[i][a][b]=max(dp[i][a][b],dp[i-1][b][c]+sum[a])(s[a]&s[b]==0 && s[b]&s[c]==0 && s[a]&s[c]==0) 表示第i行状态为s[a],并且i-1行状态为s[b]能够放置最多的大炮数量,如果i-2行的状态为s[c],并且三个状态能够兼容(即不存在有在同列的大炮),那么就i-1行状态为s[b],i-2行为状态s[c]的这个决策(dp[i-1][b][c])就能够转移到决策dp[i][a][b]
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
using namespace std;
#define MAXN 11
int dp[MAXN*10][100][100];
int line[MAXN*10],sum[MAXN*10],s[MAXN*10];
string str[MAXN*10];
int t;
bool OK(int x)
{
if(x&(x>>1)) return false;
if(x&(x>>2)) return false;
return true;
}
int getsum(int x)
{
int ret=0;
while(x)
{
ret+=(x&1);
x>>=1;
}
return ret;
}
void pre_solve(int m)
{
for(int i=0; i<(1<<m); i++)
if(OK(i))
{
s[t]=i;
sum[t]=getsum(i);
t++;
}
}
int main()
{
int n,m;
while(cin>>n>>m)
{
memset(dp,0,sizeof(dp));
memset(line,0,sizeof(line));
for(int i=1; i<=n; i++) cin>>str[i];
for(int i=1; i<=n; i++)
for(int j=0; j<m; j++)
line[i]|=(str[i][j]=='P'?1:0)<<j;
t=0;
pre_solve(m);
for(int i=0; i<t; i++)
if((s[i]&line[1])==s[i])
dp[1][i][0]=sum[i];
for(int i=2; i<=n; i++)
for(int a=0; a<t; a++)
if((s[a]&line[i])==s[a])
{
for(int b=0; b<t; b++)
if((s[b]&line[i-1])==s[b])
for(int c=0; c<t; c++)
if((s[c]&line[i-2])==s[c])
{
if(!(s[a]&s[b])&&!(s[a]&s[c])&&!(s[b]&s[c]))
{
dp[i][a][b]=max(dp[i][a][b],dp[i-1][b][c]+sum[a]);
}
}
}
int ans=0;
for(int a=0; a<t; a++)
for(int b=0; b<t; b++)
ans=max(ans,dp[n][a][b]);
cout<<ans<<endl;
}
return 0;
}
POJ1185 - 炮兵阵地(状态压缩DP)的更多相关文章
- POJ1185炮兵阵地(状态压缩 + dp)
题目链接 题意:给出一张n * m的地图,其中 有的地方能放大炮,有的地方不能,大炮与上下左右两个单位范围内会相互攻击,问最多能放几个大炮 能放大炮为1不能放大炮为0,把每一行看做一个状态,要除去同一 ...
- POJ1185 炮兵阵地 状态压缩
因为不知道不同的博客怎么转,就把别人的复制过来了,这个题解写的非常好,原地址为: http://hi.baidu.com/wangxustf/item/9138f80ce2292b8903ce1bc7 ...
- luogu2704 炮兵阵地 状态压缩DP
题目大意:一个N*M的地图由N行M列组成,地图的每一格可能是山地(用“H” 表示),也可能是平原(用“P”表示),在每一格平原地形上最多可以布置一支炮兵部队,能攻击到的区域:沿横向左右各两格,沿纵向上 ...
- poj 1185 炮兵阵地 状态压缩dp
思路:定义一个三维数组dp[x][i][j]其中x为now和pre两种状态,now表示当前两行最优解,pre表示出了本行外,前两行的最优解.那么状态转移方程为 dp[now][j][k]=max(dp ...
- POJ 3254 炮兵阵地(状态压缩DP)
题意:由方格组成的矩阵,每个方格可以放大炮用P表示,不可以放大炮用H表示,求放最多的大炮,大炮与大炮间不会互相攻击.大炮的攻击范围为两个方格. 分析:这次当前行的状态不仅和上一行有关,还和上上行有关, ...
- POJ - 1185 炮兵阵地 (状态压缩)
题目大意:中文题目就不多说大意了 解题思路: 1.每行最多仅仅有十个位置,且不是山地就是平原,那么就能够用1表示山地,0表示平原,将每一行的状态进行压缩了 2.接着找出每行能放炮兵的状态.先不考虑其它 ...
- POJ1185 炮兵阵地 —— 状压DP
题目链接:http://poj.org/problem?id=1185 炮兵阵地 Time Limit: 2000MS Memory Limit: 65536K Total Submissions ...
- [P2704][NOI2001]炮兵阵地 (状态压缩)
最近抄状压的代码…… 然后盯上了这个题目 调试了一个晚上,终于A了 但是是对着宝典打的,我依然不懂状态压缩 那么下一步先把装压放一放,学一下树形DP吧 #include<cstdio> # ...
- poj1185 炮兵阵地 状压dp
司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队.一个N*M的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用"P"表示) ...
随机推荐
- 中国首个 SaaS 模式的云告警平台安卓版 APP 上线
今年一月底,国内首个 SaaS 模式的云告警平台 OneAlert 正式发布了 iOS 版 App 客户端,今天上午,安卓版 App 客户端也正式上线了!每个安卓用户,无需电脑,都可以通过手机全程跟踪 ...
- C++11新特性:自动类型推断和类型获取
声明:本文是在Alex Allain的文章http://www.cprogramming.com/c++11/c++11-auto-decltype-return-value-after-functi ...
- Servlet课程0426(十二)Servlet MV模式下用户登录及查看用户表中所有用户
Welcome.java //登录界面 package com.tsinghua; import javax.servlet.http.*; import java.io.*; import java ...
- 数据库存储安全之(MD5+盐)加密
一般系统数据库密码加密方式: MD5后存入数据库 SHA1 Hash后存入数据库 缺点:黑客可以通过密码暴力破解获取密码信息,具体做法是将常用密码进行Hash后做成一个字典, 破解的时候,只需要查字典 ...
- RTDX target application does not match emulation protocol!
2013-06-20 10:19:22 在CCS2.0 的emulator写dsp/bios 的程序,编译链接无错误,而点击LOAD Program下载xxx.out完成时弹出如下对话框: RTDX ...
- Scrapy在win7 32位的安装及依赖包
Scrapy,一个网络爬虫的框架,首先第一步肯定是安装. 参考网上的文章. 安装过程中需要用到pip工具,请自行安装. 1.安装python 这个是必须的,既然都用到scrapy了,肯定已经安装了py ...
- 解决ImportError: cannot import name HTTPConnection的方法
在写python程序的时候,使用from httplib import HTTPConnection,在run的时候提示ImportError: cannot import name HTTPConn ...
- Emmet快速开发
标签元素关系展开 div.wrap>div.content>(div.inner_l+div.inner_r)^div.sider ------缩写展开如下---------------- ...
- poj3592Instantaneous Transference(tarjan+spfa)
http://poj.org/problem?id=3592提交了30多次了 受不了了 两份的代码基本上一样了 一个AC一个WA 木办法 贴份别人的吧 改得跟我得一样 人家能A 我是WA.. 强连通 ...
- 结构体buf_block_t
/** Buffer block for which an uncompressed page exists */ typedef struct buf_block_struct buf_block_ ...