【NOI2001】炮兵阵地

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

一看到这种题，就有种熟悉的感觉——不会是状压dp吧。于是我们可以考虑往状压方面想，这道题最多有10行，于是就最多只有2^10种状态，但是题目有限制，所以绝对不可能有那么多种状态，最多不超过100种。于是，我们可以设F[i][j][k]表示第i行，前1行状态为j，当前状态为k，进行dp。
F[i][k][j]=max(F[i][k][j],F[i-1][l][k]+sum[j]);
其中可以用sum存下每种状态下他的炮塔数量。
l表示前2行的状态。

具体实现还要有一点技巧，包括处理是否与山峰发生冲突，与判断出合法的情况，求出sum数组。

具体代码如下：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#define N 200
#define fo(i,a,b) for (int i=a;i<=b;i++)
using namespace std;
int n,m,G[N],b[N],F[N][N][N],sum[N],tot=0;
char ch;
bool pd(int x)
{
    if (x&(x<<1)) return false;
    if (x&(x<<2)) return false;
    return true;
}
int GetSum(int x)//求这个状态中有多少个1，即求炮塔数。
{
    int S=0;
    while (x)
    {
        if (x&1) S++;
        x>>=1;
    }
    return S;
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,n) 
        fo(j,1,m)
        {
            cin>>ch;
            if (ch=='H') G[i]+=1<<(j-1);
        }   
    int p=(1<<m)-1;
    tot=0;
    fo(i,0,p)
        if (pd(i))
        {
            b[++tot]=i;
            sum[tot]=GetSum(i);
        }

    memset(F,255,sizeof(F));
    fo(i,1,tot)
    if (!(G[1]&b[i])) F[1][1][i]=sum[i];

    fo(i,1,n)
        fo(j,1,tot)
        if (!(G[i]&b[j]))
            fo(k,1,tot)
            if (!(b[j]&b[k]))
            {
                fo(l,1,tot) 
                if (!(b[j]&b[l]))
                {
                    int q=F[i-1][l][k]+sum[j];
                    if (F[i-1][l][k]!=-1) F[i][k][j]=max(F[i][k][j],q);
                }
            }
    int ans=0;
    fo(i,1,tot) 
        fo(j,1,tot) ans=max(ans,F[n][i][j]);
    printf("%d",ans);
}