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+状态压缩+位运算+滚动数组
状态压缩:将每一行进行状态压缩为一个(0~2^m-1)的数字,这样可摆放状态(即左两格和右两格都不
摆放炮兵阵地),经估算大约有60个左右,再将山地的01编码转换成十进制数,用a数组记
下每一行山地的状态。
滚动数组:由于dp数组需要开三维,如果开成dp[n][1<<m][1<<m]占用内存太大,有可能爆内存,又dp
数组中dp[i][j][k] 只与 dp[i-1][k][g] 有关所以可以将数组简化成dp[2][1<<m][1<<m],
状态转移时利用迭代变量roll进行滚动,当前状态roll 与 (roll+1)%2 有关。
这题写好了之后,改正细节错误后,还是wa,后来看了下数据范围,发现会有几组特殊数据,测了下,发现
有问题,即当n<2时,dp的初始化需要进行改动,把这个改好了之后,果断AC了!!!
Codeing:
# include<stdio.h>
# include<string.h>
int max(int x,int y)
{
return (x>y?x:y);
}
int dp[2][1<<10][1<<10];
int main()
{
char ch;
int n,m,i,j,sc,k,g,roll;
int a[105],s[65],c[65];
while (~scanf("%d%d",&n,&m))
{
memset(a,0,sizeof(a));
for (i=0;i<n;i++)
{
int k,ss;
getchar();
for (j=0;j<m;j++)
{
scanf("%c",&ch);
if (ch=='P') k=0;
else k=1;
a[i]=a[i]*2+k; //记录每行山地的状态
}
}
sc=0;
for (i=0;i<(1<<m);i++)
{
int k=i;
if ((k<<1)&i||(k<<2)&i) continue;//判断该状态时候符合互不攻击的要求
c[sc]=k%2;
while (k=(k>>1)) c[sc]+=k%2; //记录此状态下的炮兵阵地数量
s[sc++]=i; //记录符合题意该状态
}
memset(dp,-1,sizeof(dp));
if (n>1)
for (i=0;i<sc;i++)
{
if (s[i]&a[1]) continue;
for (j=0;j<sc;j++)
{
if ((s[i]&s[j])||(s[j]&a[0])) continue;
dp[1][s[i]][s[j]]=c[i]+c[j];
}
}
else if (n>0)
{
for (j=0;j<sc;j++)
{
if (s[j]&a[0]) continue;
dp[1][s[j]][s[0]]=c[j];
}
} //dp的初始化,特别注意n<2时的情况
roll=1; //迭代变量
for (i=2;i<n;i++)
{
roll=(roll+1)%2;
for (j=0;j<sc;j++)
if ((s[j]&a[i])==0)
for (k=0;k<sc;k++)
if ((s[j]&s[k])==0&&(s[k]&a[i-1])==0)
for (g=0;g<sc;g++)
{
if ((s[j]&s[g])||(s[k]&s[g])||(s[g]&a[i-2])) continue;
if (dp[roll][s[j]][s[k]]==-1) dp[roll][s[j]][s[k]]=dp[(roll+1)%2][s[k]][s[g]]+c[j];
else dp[roll][s[j]][s[k]]=max(dp[roll][s[j]][s[k]],dp[(roll+1)%2][s[k]][s[g]]+c[j]); //状态转移核心方程
}
}
int sum=0;
for (i=0;i<sc;i++)
for (j=0;j<sc;j++)
if (dp[roll][s[i]][s[j]]>sum) sum=dp[roll][s[i]][s[j]];
printf("%d\n",sum);
}
return 0;
}
预祝AC!!!