方格有m*n个格子,一共有2^(m+n)种排列,很显然不能使用暴力法,因而选用动态规划求解.

求解DP问题一般有3步,即定义出一个状态 求出状态转移方程 再用算法实现.多数DP题难youguan点在于第2步,而在状态压缩DP中,定义状态也是很关键的一个步骤.有关位运算的基础知识,按位与,按位或,异或等可自行查阅资料,这里仅作简单说明.

<<n == 2的n次方
>>n == /2的n次方
(n>>k) & // 取出整数n在二进制下的第k位
n & ((<<k)-) // 取出整数n在二进制下的后k位
(i>>j) & k // i右移j位后和k与运算

很容易想到用二进制数来表示方格,1表示放炮兵,0表示不放.在同一行中,只要没有出现两个炮兵紧邻或者两个炮兵只间隔1个位置的情况,均是合法的状态.故在二进制表示的行01串中删除字串含有"11","101"的原串即可,预处理出合法的01串并存于legal中.

vector<int> legal;
void init() { // 找到合法的摆放总数
for(int i = ; i < (<<m); i++) { // 1<<m == pow(2,m),遍历所有情况
int c1 = , c2 = ; // 3 -> (11) , 5 -> (101)
bool sub = ;
for(int j = ; j < m - ; j++) {
if(((i >> j) & c1) == c1) {
sub = ;
}
}
for(int j = ; j < m - ; j++) {
if(((i >> j) & c2) == c2) {
sub = ;
}
}
if(sub) legal.push_back(i);
}
}

并用count函数计算每行中的1(炮兵数目):

int count(int x) {
int cnt = ;
for(int i = ; i < m; i++) {
if(((x>>i)&) == ) {
cnt++;
}
}
return cnt;
}

接下来设计DP状态和状态转移

每一行的状态受该行前两行摆放状态的影响,因此选择dp[i][j][k]表示可行方案数.dp[i][j][k]表示第i行压缩后状态为j,第i-1行压缩后状态为k的情况下前i行最多放多少个炮兵.同时由于dp[102][1050][1050]会MLE,只有80分不过也知足了, 需要将dp改为滚动数组,将第一维每处均%3即可,只记录前两次的状态

读入地图时把每行压缩成一个二进制数,为了便于后续查找可行状态,读入时H对应1,P对应0. 后续遍历了legal中的合法状态时,若与该行的01表示与运算值不为0,即与地图存在冲突.

状态转移方程为dp[本行][前一行][再前一行] = max{dp[前一行][再前一行][再前一行的前一行] + count(本行) , dp[本行][前一行][再前一行]}

完整的代码如下:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + ; int n, m;
string s;
int mp[];
vector<int> legal; // 储存所有的合法01串
int dp[][][]; void init() { // 找到合法的摆放总数
for(int i = ; i < (<<m); i++) { // 1<<m == pow(2,m)
int c1 = , c2 = ; // 3 -> (0011) , 5 -> (0101)
bool sub = ;
for(int j = ; j < m - ; j++) {
if(((i >> j) & c1) == c1) {
sub = ;
}
}
for(int j = ; j < m - ; j++) {
if(((i >> j) & c2) == c2) {
sub = ;
}
}
if(sub) legal.push_back(i);
}
} int count(int x) {
int cnt = ;
for(int i = ; i < m; i++) {
if(((x>>i)&) == ) {
cnt++;
}
}
return cnt;
} int main() {
cin>>n>>m;
init();
// cout<<legal.size()<<endl;
for(int i = ; i <= n + ; i++) { // i初值2 避免越界(需考虑到前两行)
cin>>s;
for(int j = ; j < m; j++)
if(s[j] == 'H') mp[i] |= (<<j);
}
for(int i = ; i <= n + ; i++) {
for(auto step : legal) {
if((step & mp[i]) != ) continue;
for(auto bst : legal) {
if((step & bst) != ) continue;
if((bst & mp[i-]) != ) continue;
for(auto bbst : legal) {
if((step & bbst) != ) continue;
if((bbst & mp[i-]) != ) continue;
if((bbst & bst) != ) continue;
dp[i%][step][bst] = max(dp[(i-)%][bst][bbst] + count(step), dp[i%][step][bst]);
}
}
}
}
int res = ;
for(auto step : legal)
for(auto bst : legal)
res = max(res, dp[(n+)%][step][bst]);
cout<<res<<endl;
}

luogu P2704 炮兵阵地(经典状态压缩DP)的更多相关文章

  1. poj 1185 炮兵阵地 [经典状态压缩DP]

    题意:略. 思路:由于每个大炮射程为2,所以如果对每一行状态压缩的话,能对它造成影响的就是上面的两行. 这里用dp[row][state1][state2]表示第row行状态为state2,第row- ...

  2. 洛谷 P2704 [NOI2001]炮兵阵地 (状态压缩DP+优化)

    题目描述 司令部的将军们打算在NM的网格地图上部署他们的炮兵部队.一个NM的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用"P" ...

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

    POJ飞翔.数据弱 ZQOJ飞翔 数据强 Description 司令部的将军们打算在N×M的网格地图上部署他们的炮兵部队.一个N×M的地图由N行M列组成,地图的每一格可能是山地(用"H&q ...

  4. POJ 1185 炮兵阵地 (状态压缩DP)

    题目链接 Description 司令部的将军们打算在NM的网格地图上部署他们的炮兵部队.一个NM的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用& ...

  5. POJ 1185 炮兵阵地(状态压缩DP)

    题解:nState为状态数,state数组为可能的状态 代码: #include <map> #include <set> #include <list> #inc ...

  6. poj 1185 炮兵阵地(三维状态压缩dP)

    题目:http://poj.org/problem?id=1185 思路: d[i][j][k]表示第i行的状态为第k个状态,第i-1行的状态为第j个状态的时候 的炮的数量. 1表示放大炮, 地形状态 ...

  7. POJ 2411 Mondriaan's Dream [经典状态压缩dp]

    题意:略. 思路:这一题开始做的时候完全没有思路,便去看了别人的题解. 首先,对于这个题目解法想有一个初步的了解,请看这里:http://www.2cto.com/kf/201208/146894.h ...

  8. poj -1185 炮兵阵地 (经典状压dp)

    http://poj.org/problem?id=1185 参考博客:http://poj.org/problem?id=1185 大神博客已经讲的很清楚了,注意存状态的时候是从1开始的,所以初始化 ...

  9. 【Luogu】P2704炮兵阵地(状压DP)

    题目链接 话说还真没见过能影响两行的状压.想了半天想出来f数组再多一维就能表示,但是没想到怎么才能不爆空间…… 也是从这道题里学到的一个妙招. 可以把合法状态存到一个数组里,然后用数组下标来映射状态. ...

随机推荐

  1. 大小端(MSB & LSB)

    谈到字节序的问题,必然牵涉到两大CPU派系.那就是Motorola的PowerPC系列CPU和Intel的x86系列CPU.PowerPC系列采用big endian方式存储数据,而x86系列则采用l ...

  2. java判断两个时间相差得天数

    方法一:通过Calendar类得日期比较,在这需要考虑闰年和平年,也要考虑跨年份 /** * date2比date1多的天数 * @param date1 * @param date2 * @retu ...

  3. Azure IoT Hub 十分钟入门系列 (1)- 10分钟带你了解Azure IoT Hub 并创建IoT Hub

    建议您先对<Azure 上 IoT 整体解决方案概览 >进行了解. 本文主要分享一个案例: 10分钟-了解Azure IoT Hub并创建Azure IoT Hub 本文主要有如下内容: ...

  4. OracleDBConsoleorcl 服务无法启动:Agent process exited abnormally during initialization.

    OracleDBConsoleorcl 服务无法启动 在事件查看器里看到 Agent process exited abnormally during initialization.的记录.知道是因为 ...

  5. Django+Celery+redis kombu.exceptions.EncodeError:Object of type is not JSON serializable报错

    在本文中例子中遇到问题的各种开发版本如下: Python3.6.8 Django==2.2 celery==4.4.0 kombu==4.6.7 redis==3.3.0 大概的报错如下截图: 是在开 ...

  6. ehcache注解全面解析

    通过ehcache以编程方式使用缓存: 跟上面的方式相同,但是缓存通过ehcache去管理,当然比使用map有N多种好处,比如缓存太大了快达到上限之后,将哪一部分缓存清除出去.这种方式完全是通过代码的 ...

  7. TensorFlow使用RNN实现手写数字识别

    学习,笔记,有时间会加注释以及函数之间的逻辑关系. # https://www.cnblogs.com/felixwang2/p/9190664.html # https://www.cnblogs. ...

  8. 【PAT甲级】1099 Build A Binary Search Tree (30 分)

    题意: 输入一个正整数N(<=100),接着输入N行每行包括0~N-1结点的左右子结点,接着输入一行N个数表示数的结点值.输出这颗二叉排序树的层次遍历. AAAAAccepted code: # ...

  9. VMware15下载、安装、激活

    1.VMware15下载 链接:https://pan.baidu.com/s/1bI8LReRY-5k81O3rrNgg-A  提取码:6c03 2.VMware15安装 3.VMware15激活

  10. mybatis Plus 多表联合查询

    //实体类package com.sk.skkill.entity; import com.baomidou.mybatisplus.annotation.TableField;import com. ...