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写。
状态压缩dp是指用一个数的二进制表示状态,例如“PHPP”表示为1011,即11,因为每行不超过10个,这样子最多用2的10次方即1024可以表示所有的状态。
(1)对于每一行的状态j,需要保证周围两格内均为0,即j的二进制数中任何一个1左右至少两个数为0,即满足:
j & (j<<1) == 0 && j & (j<<2) == 0
(2)因为要考虑上下两行的状态,所以我们用dp[i][k][j]表示第i行处于j状态,且上一行状态为k的时候最大能放置的炮兵数目。
(3)对于任意三行的状态j、k、t,要保证不能有同一列不能有重复的1,即:
(j & k == 0) && (k & t == 0) && (j &t == 0)
(4)因为要开三维数组,所以使用1024会超过内存限制,要先处理每一行可以满足条件j & (j<<1) == 0 && j & (j<<2) == 0的数保存到一个数组中,以数组的下标作为对应的状态,可以计算满足这样子的条件的数字为60,则不会超过内存限制。
(5)dp[i][k][j] = max(cot(num[j]) + dp[i-1][t][k]);
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
int a[105];
int num[70];
int dp[105][70][70];
int top;
void cal(int x)
{
top = 0;
for(int i = 0; i < (1<<x); i++)
{
if(!(i&(i<<1))&&!(i&(i<<2)))
num[top++] = i;
}
}
int cot(int x)
{
int cnt = 0;
while(x)
{
cnt++;
x &= (x-1);
}
return cnt;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int n,m;
char t;
while(cin >> n >> m)
{
memset(a,0,sizeof(a));
for(int i = 1; i <= n; i++)
{
for(int j = 1; j <= m; j++)
{
cin >> t;
if(t == 'P')
a[i] += 1<<(m-j);
}
}
cal(m);
memset(dp,0,sizeof(dp));
for(int i = 1; i <= n; i++)
{
for(int j = 0; j < top; j++)
{
if((num[j]|a[i])<=a[i])
{
if(i==1)
dp[i][0][j] = cot(num[j]);
else if(i == 2)
{
for(int k = 0; k < top; k++)
{
if(!(num[j]&num[k])&&(num[k]|a[i-1])<=a[i-1])
{
dp[i][k][j] = cot(num[k]) + cot(num[j]);
}
}
}
else
{
for(int k = 0; k < top; k++)
{
if(!(num[j]&num[k])&&(num[k]|a[i-1])<=a[i-1])
{
for(int t = 0; t < top; t++)
{
if(!(num[j]&num[t]) && !(num[t]&num[k]) &&(num[t]|a[i-2])<=a[i-2])
{
int temp = cot(num[j]) + dp[i-1][t][k];
if(temp > dp[i][k][j])
dp[i][k][j] = temp;
}
}
}
}
}
}
}
}
int ans = 0;
for(int i = 0; i < top; i++)
{
for(int j = 0; j < top; j++)
{
if(dp[n][j][i] > ans)
ans = dp[n][j][i];
}
}
cout << ans << endl;
}
return 0;
}