poj 1185 炮兵阵地 （状态压缩DP）

炮兵阵地

Time Limit: 2000MS Memory Limit: 65536KTotal Submissions: 19690 Accepted: 7602

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 4PHPPPPHHPPPPPHPPPHHP

Sample Output

6

不容易啊，debug了一下午。。。

</pre><pre name="code" class="cpp">#include<iostream>#include<stdio.h>#include<math.h>#include<string.h>#include<algorithm>#include<iostream>using namespace std;#define N 105const int inf=0x3fffffff;char g[N][15];int num[65],cnt[65],cur[N];int top,dp[N][65][65];int getsum(int x){    int t=0;    while(x)    //得到X转换成为二进制后1的个数    {        x&=(x-1);        t++;    }    return t;}void inti(int m){    int i;    top=0;    memset(num,-1,sizeof(num));    for(i=0;i<(1<<m);i++)    {        if(i&(i<<1)) continue;        if(i&(i<<2)) continue;        num[top]=i;        cnt[top++]=getsum(i);    }}int main(){    int i,j,k,r,n,m;    inti(10);    while(scanf("%d%d",&n,&m)!=-1)    {        for(i=0;i<n;i++)        {            scanf("%s",g[i]);            cur[i]=0;            for(j=0;j<m;j++)            {                if(g[i][j]=='H')                    cur[i]|=(1<<j);  //得到每一行的不可放置大炮的信息，              }        }        for(top=0;num[top]!=-1&&num[top]<(1<<m);top++)            ;                     //因为当m=10时会访问没有赋值的数组元素        memset(dp,-1,sizeof(dp));        for(i=0;i<top;i++)        {            if(num[i]&cur[0]) continue;            dp[0][i][0]=cnt[i];        }        for(r=1;r<n;r++) //枚举剩下的每一行          {            for(i=0;i<top;i++) //每一行可取的每种组合              {                if(cur[r]&num[i]) continue;                for(j=0;j<top;j++) //上一行（i-1）的每种组合                 {                    if(num[i]&num[j]) continue;                    for(k=0;k<top;k++)  //枚举i-2行的每种组合                      {                        if(num[i]&num[k]) continue;                        if(num[j]&num[k]) continue;   //改成dp[r-1][j][k]==-1速度更快                        dp[r][i][j]=max(dp[r][i][j],dp[r-1][j][k]+cnt[i]);                    }     //当前行和第j、k行不冲突                  }            }        }        int ans=0;        for(i=0;i<top;i++)        {            for(j=0;j<top;j++)                ans=max(ans,dp[n-1][i][j]);        }        printf("%d\n",ans);    }    return 0;}