poj 1185

 

炮兵阵地

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions: 9240		Accepted: 3270

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

Source

Noi 01

分析：状态压缩DP。由于第i行的状态只和前两行有关，所以用f[i][s1][s2]表示第i行状态为s1，i-1行状态为s2时的最优解，容易得到状态转移方程f[i][s1][s2]=max(f[i][s1][s2],f[i-1][s2][s3]+w[s1]),(s1&(s2|s3)==0),如果纯枚举状态s，肯定超时，所以先求出每一行可能的状态，并保存起来，还可以加一个优化，就是存储上下两行哪些状态是可以连接的，不过貌似效果不明显。。。‍

很挫的代码，谁帮我优化啊：

#include<cstdio> #include<cstring> using namespace std; inline int max(int a,int b){return a>b?a:b;} int p[11]={1,2,4,8,16,32,64,128,256,512,1024}; char w[1025]={0},map[11]; int f[101][61][61],g[101][61]={0},q[101][61][61]; int i,j,k,l,n,m; void dfs(int l,int s) { g[i][++g[i][0]]=s; if(l>=m)return; for(int k=l;k<m;++k) if(map[k]=='P')dfs(k+3,s|p[k]); } int main() { int i1,i2,j1,j2,ans; for(i=0;i<1025;++i) for(j=0;j<11;++j) if(i&p[j])++w[i]; while(scanf("%d%d",&n,&m)!=EOF) { memset(f,0,sizeof(f)); scanf("%s",map); g[0][0]=1,g[i=1][0]=0,dfs(0,0); for(j=1;j<=g[i][0];++j)f[i][j][1]=w[g[i][j]],q[i][j][0]=q[i][j][1]=1; for(i=2;i<=n;++i) { scanf("%s",map); g[i][0]=0,dfs(0,0); i1=i-1; i2=i-2; for(j=1;j<=g[i][0];++j) for(q[i][j][0]=0,k=1;k<=g[i1][0];++k) if((g[i][j]&g[i1][k])==0)q[i][j][++q[i][j][0]]=k; for(j=1;j<=g[i][0];++j) for(k=1;k<=q[i][j][0];++k) for(j1=q[i][j][k],l=1;l<=q[i1][j1][0];++l) if((g[i][j]&g[i2][j2=q[i1][j1][l]])==0) f[i][j][j1]=max(f[i][j][j1],f[i1][j1][j2]+w[g[i][j]]); } for(ans=0,i=1;i<=g[n][0];++i) for(j=1;j<=q[n][i][0];++j) ans=max(ans,f[n][i][q[n][i][j]]); printf("%d/n",ans); } return 0; }