poj 1185 炮兵 状态压缩dp

司令部的将军们打算在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

转：http://blog.sina.com.cn/s/blog_6635898a0100qe3w.html

思路：状态压缩DP。很经典的状态压缩DP。用int型来表示每行的状态（如果int型的二进制的第i位为1，则表示这一行的第i列有安装大炮）。这样的话由于最多有10列，故由计算可得最多有60种状态。DP部分：dp[r][i][k]表示第r行的状态为k，第r-1行的状态为i时候，前r行最多能够安装的大炮数量。

解读了一下大佬的代码。。 强。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int dp[1010][70][70];
int mp[1010];
int cnt=0;
int sum[100];
int stk[1100];
int ok(int x)
{
    if(x&(x<<1)) return false;  //如果两个炮台相差距离小于2 那么不行。
    if(x&(x<<2)) return false;
    return true;
}



int getSum(int x)
{
    int num=0;
    while(x>0)
    {
        if(x&1) num++;
        x>>=1;
    }
    return num; //得到当前状态可以放多少个炮台
}


void findStk(int n)
{

    for(int i=0;i<(1<<n);i++)
    {
        if(ok(i))
        {
            stk[cnt]=i; //把可行的状态存起来
            sum[cnt++]=getSum(i);//存放可行的状态的炮台数。
        }
    }
}
int main()
{
    int row,col,r,c;
    cin>>row>>col;
    memset(dp,-1,sizeof(dp));
    for(int r=0;r<row;r++)
    {
        for(c=0;c<col;c++)
        {
            char cmp;
            cin>>cmp;
            if(cmp=='H') mp[r]|=(1<<c); //逆向读表 答案和正向是一样的
        }
    }
    findStk(col); //每一行的所有状态是一样，然后再处理每一个可行的状态，先预处理
    for(int i=0;i<cnt;i++)
    {
        if(!(stk[i]&mp[0]))
            dp[0][0][i]=sum[i];
    }//第一行特殊 需要预处理。
    //当前状态的炮台数可以加上上面不冲突的状态的炮台数。
    //而当前状态的炮台冲突只与上两行有关系
    //所以用三维数组 dp[r][j][i] 代表 第r行状态是i，r-1行是j状态时的炮台数。
    //为什么要用三维而不是四维呢，三维已经可以代表状态的唯一性。然后转移即可
    for(int r=1;r<row;r++)
    {
        for(int i=0;i<cnt;i++)
        {
            if(stk[i]&mp[r]) continue; //当前状态的放炮台数 是否与不能放炮台的地方冲突。
            for(int j=0;j<cnt;j++)
            {
                if(stk[i]&stk[j]) continue;// 这两行放炮台的状态有冲突或者重叠。
                for(int k=0;k<cnt;k++)
                {
                    if(stk[i]&stk[k]) continue;
                    if(dp[r-1][k][j]==-1) continue;//上面没放炮台那么
                    dp[r][j][i]=max(dp[r][j][i],dp[r-1][k][j]+sum[i]);
                    //这样是确定每一行的状态最优。
                }
            }
        }
    }
    int ans=0;
    for(int i=0;i<cnt;i++)
        for(int j=0;j<cnt;j++)
            ans=max(ans,dp[row-1][i][j]);
        printf("%d\n",ans );
}