【状压】 炮兵阵地

炮兵阵地

Time Limit: 2000MS

Memory Limit: 65536K

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

连续做了几道状压的题，终于有点感觉了。。。。。。。

这一题和ＳＣＯＩ２００５那年的互不侵犯有点相似

每一行的状态会受到它前面两行的影响，那么我的状态就定义为 f [ i ] [ s ] [ t ] ，表示当前在第 i 行，当前行的状态为 s ，上一行的状态为 t

不是说上两行吗，怎么状态里只有上一行？我们在看看转移

f [ i ] [ s ] [ t ] = max ( f [ i ] [ s ] [ t ] , f [ i - 1 ] [ t ] [ k ] + count ( s ) ) 

这样就涉及到两行了

判断可行还是比较好判断，判断 s , t , k 不互相重复即可

最直观的想法是直接用 s 和 t 来存那个01序列，但是空间就会有 100*1024*1024=100MB，但是题目要求是64MB！

由于一行内每两个1之间至少有两个空格，所以我们可以用dfs生成出所有一行可能的情况，我们把它们编号，就会发现极限情况才60种，我们用a[]数组将它存下，所以我们s和t开成60+就ok了，然后它们对应的01序列分别为 a [ s ] 和 a [ t ] ，这样就可以解决空间问题了

其他的没什么，看看代码，注释什么的都没有删

/*http://blog.csdn.net/jiangzh7 
By Jiangzh*/  
#include<cstdio>  
#include<cstring>  
#include<iostream>  
#include<algorithm>  
using namespace std;  
typedef long long LL;  
const int N=100+10;  
const int M=10+3;  
  
int n,m;  
int map[N];//map[i]存储一个01序列，1表示不能走  
bool hash[1<<M];//dfs用的标号  
int a[1<<M],Cnt;//存储所有一行内可行的01方案  
LL f[N][60+10][60+10];//用dfs验证最多60种情况，所以开60即可  
LL num[1<<M];//种01序列中1的个数  
int all;//储存全为1的序列  
  
void PT(int s)  
{  
    int a[M];  
    for(int i=m;i>=1;i--)  
    {  
        a[i]=s&1;s>>=1;  
    }  
    for(int i=1;i<=m;i++) printf("%d",a[i]);  
    puts("");  
}  
  
void read()  
{  
    scanf("%d%d",&n,&m);  
    char ss[M];  
    for(int i=1;i<=n;i++)  
    {  
        scanf("\n%s",ss);  
        for(int j=0;j<m;j++)  
        if(ss[j]=='H') map[i]|=(1<<j);  
    }  
    all=(1<<m)-1;  
}  
  
LL count_1(int s)//统计s中1的个数  
{  
    LL res=0;  
    while(s)  
    {  
        res+=s&1;  
        s>>=1;  
    }  
    return res;  
}  
  
void dfs(int S,int x)  
{  
    if(x>=m)  
    {  
        if(hash[S]) return;  
        //printf("%d :",S);PT(S);  
        hash[S]=true;  
        a[++Cnt]=S;  
        num[S]=count_1(S);  
        return;  
    }  
    dfs((S<<1)&((1<<m)-1),x+1);  
    if(!(S&1)) dfs(((S<<2)+1)&((1<<m)-1),x+1);  
}  
  
bool could(int s,int t,int k)  
{  
    if((s&t)!=0) return false;  
    if((s&k)!=0) return false;  
    if((t&k)!=0) return false;  
    return true;  
}  
  
bool could2(int s,int t){return (s&t)==0;}  
  
void work()  
{  
    dfs(0,0);//dfs出一行内所有可能情况  
    for(int s=1;s<=Cnt;s++)//单独处理第一行  
    {  
        //printf("s : "); PT(a[s]);  
        if(a[s]&map[1]) continue;  
        f[1][s][1]=num[a[s]];  
    }  
    //puts("-----------------");  
    for(int s=1;s<=Cnt;s++)//单独处理第二行  
    {  
        //printf("s : "); PT(a[s]);  
        if(a[s]&map[2]) continue;  
        for(int t=1;t<=Cnt;t++)  
        if(could2(a[s],a[t]))  
        {  
            f[2][s][t]=num[a[s]]+num[a[t]];  
        }  
    }  
    //puts("-----------------");  
    for(int i=3;i<=n;i++)  
    {  
        //printf("%d : \n",i);  
        for(int s=1;s<=Cnt;s++)  
        {  
            if(a[s]&map[i]) continue;//不能在山上驻扎 下同  
            //printf("s=%d : ",a[s]);PT(a[s]);  
            for(int t=1;t<=Cnt;t++)  
            {  
                if(a[t]&map[i-1]) continue;  
                //printf("t=%d : ",a[t]);PT(a[t]);  
                for(int k=1;k<=Cnt;k++)  
                {  
                    if(a[k]&map[i-2]) continue;  
                    //printf("k=%d : ",a[k]);PT(a[k]);  
                    if(could(a[s],a[t],a[k]))  
                    {  
                        f[i][s][t]=max(f[i][s][t],f[i-1][t][k]+num[a[s]]);  
                        //printf("%d\n",f[i][s][t]);  
                    }  
                }  
            }  
        }  
    }  
    LL res=0;  
    for(int s=1;s<=Cnt;s++)  
    {  
        for(int t=1;t<=Cnt;t++)  
        {  
            res=max(res,f[n][s][t]);  
            //printf("%d ",f[n][s][t]);  
        }  
        //puts("");  
    }  
    cout<<res<<endl;  
}  
  
int main()  
{  
    freopen("poj1185.in","r",stdin);  
    freopen("poj1185.out","w",stdout);  
    read();  
    work();  
    return 0;  
}