POJ 1185 (状态DP)

                                                                                                                            炮兵阵地

Time Limit: 2000MS Memory Limit: 65536KTotal Submissions: 11753 Accepted: 4239

Description

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示：

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。

Input

第一行包含两个由空格分割开的正整数，分别表示N和M；
接下来的N行，每一行含有连续的M个字符('P'或者'H')，中间没有空格。按顺序表示地图中每一行的数据。N <= 100；M <= 10。

Output

仅一行，包含一个整数K，表示最多能摆放的炮兵部队的数量。

第一次做状态压缩DP，A的很不容易呀。。。。。

F[I][J][K] 第I行状态为J第I-1行状态为K时可得到的最大值

#include<stdio.h>#include<string.h>#define clear(a,b) memset(a,b,sizeof(a))#define max(a,b) ((a)>(b)?(a):(b))#define maxn1 70#define maxn2 110int cas[maxn1],map[maxn2],sum[maxn1],f[maxn2][maxn1][maxn1],m,n,l = 0;int cal(int x){  int s = 0,i;  for(i = 0;i < m;i++)     s += (x>>i)&1;  return s;}int judge(int x){  int i;  for(i = 0;i < m;i++) {     if (x & (x>>1) & 1) return 0;     if (x & (x>>2) & 1) return 0;     x = x >> 1;     }  return 1;}int judge2(int x1,int x2,int x3){  if (x1 & x2) return 0;  if (x1 & x3) return 0;  if (x2 & x3) return 0;  return 1;}void pre(){   int i,j;   clear(f,0);   clear(cas,0);   clear(sum,0);   clear(map,0);   scanf("%d%d",&n,&m);   for(i = 1;i <= n;i++)      for(j = 1,getchar();j <= m;j++)         if (getchar() == 'H') map[i] = map[i] | (1<<(m-j));   int mi = (1 << m) - 1;   for (i = 0;i <= mi;i++)      if (judge(i)) cas[++l] = i;   for(i = 1;i <= l;i++)      sum[i] = cal(cas[i]);}int work(){  int i,j,k,t,ans = 0;   for(i = 1;i <= l;i++)      for(j = 1;j <= l;j++)         if ((map[1] & cas[i]) == 0)  {f[1][i][j] = sum[i];ans = max(f[1][i][j],ans);}    for(i = 2;i <= n;i++)       for(j = 1;j <=l;j++)          for(k = 1;k <= l;k++)            for(t = 1;t <= l;t++) {              if (!(cas[j] & map[i])  && judge2(cas[t],cas[j],cas[k])) f[i][j][k] = max(f[i][j][k],f[i-1][k][t] + sum[j]);              ans = max(ans,f[i][j][k]);              }    printf("%d\n",ans);}int debug(){  int i,j;  for(i = 1;i <= n;i++,printf("\n"))    for(j = 1;j <= m;j++) printf("%d",(map[i]>>(m-j))&1);  for(i = 1;i <= l;i++){    for(j = 1;j <= m;j++) printf("%d",(cas[i]>>(m-j))&1);    printf(" ");    printf(" %d\n",sum[i]);    }}int main(){  //freopen("xx.out","w",stdout);  pre();  work();  //debug();  return 0;}