Description
司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用"P"表示),如下图。在每一格平原地形上最多可以布置一支炮兵部队(山地上不能够部署炮兵部队);一支炮兵部队在地图上的攻击范围如图中黑色区域所示:
如果在地图中的灰色所标识的平原上部署一支炮兵部队,则图中的黑色的网格表示它能够攻击到的区域:沿横向左右各两格,沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在,将军们规划如何部署炮兵部队,在防止误伤的前提下(保证任何两支炮兵部队之间不能互相攻击,即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内),在整个地图区域内最多能够摆放多少我军的炮兵部队。
![](https://i-blog.csdnimg.cn/blog_migrate/ba1d4fa736b99550a7ee43784c287df9.jpeg)
如果在地图中的灰色所标识的平原上部署一支炮兵部队,则图中的黑色的网格表示它能够攻击到的区域:沿横向左右各两格,沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在,将军们规划如何部署炮兵部队,在防止误伤的前提下(保证任何两支炮兵部队之间不能互相攻击,即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内),在整个地图区域内最多能够摆放多少我军的炮兵部队。
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
题解
状态压缩DP。在状态压缩的过程中要考虑到以下几点:
1、对于一个位置,如果它部署了一支部队,那么会对它的前后左右的2格位置产生影响,所以对于第i行的状态,我们要从i-2和i-1行转移过来,所以我们考虑用三维数组作为dp数组。
2、对于每一行,它的状态不能和地图冲突(即不能在山地上部署)。所以,我们在考虑每一行的状态时,可用位运算判断该状态是否合法。
3、如果一行只有“空地”,则每一行最多有2^10种状态,显然100*1024*1024时间和空间难以承受。
但是,可以发现:如果一行中的一个点部署了部队,则其前后两格范围内均不能再部署部队。所以在2^10种状态中有很多是不合法的。于是,我们可以证明:一行里最多有60种合法的状态。我们只要暴力枚举出所有的合法状态state[x],并记录该状态下第i行所增加的部队数count[x]即可。
4、而且在动态规划中,对于第i行的状态,我们要且仅要从i-2和i-1行转移过来,所以我们只需要开2*60*60的数组即可。
综上我们可用f[i][j][k]表示第i行用j状态,第i-1行用k状态,转移方程为f[i][j][k]=max(f[i][j][k],f[(i+1)%2][k][l]+count[j]);
方程如下:
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;
int n,m,map[102]/*记录那些地方为山地*/;
int state[65]/*记录一行内所有可行的状态*/,num/*一行内可行的状态数*/,
S/*总状态数*/,count[65]/*每个状态包含多少个炮兵*/;
int f[2][102][102];
char ch[15];
void deal(int tot)//将一行内所有可能的状态记录下来
{
for(int t=0;t<=tot;t++)
{int k=t;
if((k<<1)&t||(k<<2)&t) continue;//防止一行内互相攻击
while(k)
{count[num]+=k%2;
k=k>>1;
}
state[num++]=t;
}
//for(int i=0;i<num;i++)
//printf("%d %d\n",state[i],count[i]);
}
void dp()
{
int h=0;
for(int i=0;i<n;i++)
{for(int j=0;j<num;j++)//第i行的状态是否跟地形冲突
{if(map[i]&state[j]) continue;
if(i==0) f[h][j][0]=count[j];
else if(i==1)
{for(int k=0;k<num;k++)
{if(map[i-1]&state[k]) continue;//i-1行是否冲突
if(state[j]&state[k]) continue;//i行和i-1行是否冲突
f[h][j][k]=max(f[h][j][k],f[(h+1)%2][k][0]+count[j]);
}
}
else
{for(int k=0;k<num;k++)
{if(state[k]&map[i-1]) continue;
if(state[j]&state[k]) continue;
for(int l=0;l<num;l++)
{if(state[l]&map[i-2]) continue;
if( (state[k]&state[l]) | (state[j]&state[l]) ) continue;
f[h][j][k]=max(f[h][j][k],f[(h+1)%2][k][l]+count[j]);
}
}
}
}
h=(h+1)%2;
}
h=(h+1)%2;//求出最大值
int ans=0;
for(int i=0;i<num;i++)
for(int j=0;j<num;j++)
ans=max(ans,f[h][i][j]);
printf("%d\n",ans);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=0;i<n;i++)
{scanf("%s",ch);
int l=strlen(ch);
for(int j=0;j<l;j++)
{if(ch[j]=='H')
map[i]=map[i]+(1<<j);
}
}
num=0; S=(1<<m)-1;
deal(S);
dp();
return 0;
}