poj 1185 火炮 （减少国家DP）

Description

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。

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

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

Input

Output

Sample Input

5 4
PHPP
PPHH
PPPP
PHPP
PHHP

Sample Output

6

不easy啊，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 105
const 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;
}