POJ 1185 炮兵阵地 (状压DP)

炮兵阵地
Time Limit: 2000MS Memory Limit: 65536K
Total Submissions: 25954 Accepted: 10004
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开个二维的就ok。注意初始化，因为初始化一直T了很久。

#include "cstring"
#include "cstdio"
#include "iostream"
#include "string.h"
#include "algorithm"
#include "cmath"
using namespace std;
char mmap[200][200];
int n,m;
int now_status[2000],now_sum[2000],nowid,preid_1,preid_2,dp[2000][2000],pre_status1[2000],pre_status2[2000] ,pre_sum[2000][2000];
void dfs(int row,int column,int status,int sum)
{
    if(column>=m)
    {
        now_status[++nowid]=status;
        now_sum[nowid]=sum;
        return;
    }
    if(mmap[row][column]=='P')
        dfs(row,column+3,status|1<<column,sum+1);
    dfs(row,column+1,status,sum);
    return;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=0;i<n;i++)
                scanf("%s",mmap[i]);
        /*
        memset(pre_status1,0,sizeof(pre_status1));
        memset(pre_status2,0,sizeof(pre_status2));
        memset(now_status,0,sizeof(now_status));
        memset(now_sum,0,sizeof(now_sum));
        memset(pre_sum,0,sizeof(pre_sum));
        memset(dp,0,sizeof(dp));*/
        pre_sum[0][0]=pre_status1[0]=pre_status2[0]=0;
        preid_1=preid_2=0;
        for(int k=0;k<n;k++)
        {
            nowid=-1;
            dfs(k,0,0,0);
            memset(dp,0,sizeof(dp));
            /*
            for(int i=0;i<=nowid;i++)
                for(int j=0;j<=preid_1;j++)
                    dp[i][j]=0;*/
            for(int i=0;i<=nowid;i++)
            {
                for(int j=0;j<=preid_1;j++)
                {
                    for(int q=0;q<=preid_2;q++)
                    {
                        if(now_status[i]&pre_status1[j])
                            continue;
                        if(now_status[i]&pre_status2[q])
                            continue;
                        dp[i][j]=max(dp[i][j],pre_sum[j][q]+now_sum[i]);

                    }
                }
            }
            for(int i=0;i<=nowid;i++)
                for(int j=0;j<=preid_1;j++)
                    pre_sum[i][j]=dp[i][j];
            for(int i=0;i<=preid_1;i++)
                pre_status2[i]=pre_status1[i];
            for(int i=0;i<=nowid;i++)
                pre_status1[i]=now_status[i];
            preid_2=preid_1;
            preid_1=nowid;
        }
        int ans=0;
        for(int i=0;i<=preid_1;i++)
            for(int j=0;j<=preid_2;j++)
                ans=max(ans,dp[i][j]);
        printf("%d\n",ans);
    }
    return 0;
}