poj1185——炮兵阵地

题目大意：n*m的网格分布着山地H和平原P，只有平原可以放一只炮兵部队，每个部队能攻击自己上下左右四个方向相邻两个网格，任何一只部队不能放在其他部队的攻击区域内，问最多能放多少只炮兵部队

输入：n m

           网格HP分布

输出：炮兵部队的最大数量

分析：动态规划，经典的状态压缩（二进制压缩成十进制）。第i行能放多少炮兵部队，取决于第i-1、i-2行，而对于每一行，每个网格只有放与不放（01）两种选择，如果简单的暴力dp[i-1][][][][][]...dp[i-2][][][][][]...太大了，所以选择状态压缩，将01的二进制状态串压缩成十进制数。

           原地图中：0表示平原1表示山地，即平原0可以在放置部队后变为1

           先不考虑会不会放到山上，先二进制枚举放置部队的合法情况。int state[MaxM]：保存部队的合法放置状态，即保证部队之间不相互攻击。显然这些状态不能取到某些值，因为某些二进制表示部队相邻，不符合要求。

           枚举方法：for(i=0;i<(i<<col);i++) 

                                 if(!(i&(i<<1))&!(i&(i<<2)))

                                     i为合法状态（可以转化为十进制储存在state数组中）

           int soldier[MaxM]：在state[]状态下放了多少个部队（数出state十进制对应的二进制有几个1）

           int base[MaxR]：原地图中第i行的状态

           然后再判断某个state[i]状态能不能放到base[r]（第r行）上。判断方法：state[i]&base[r]=0表示可以放。

           同理判断i和i-1、i-2行冲突与否：state[i]&state[i-1]=0,state[i]&state[i-2]=0,state[i-1]&state[i-2]=0，三者同时满足才是不冲突。

           做题步骤：1.计算原地图每行的状态base[r]

                             2.枚举部队不相互攻击的状态

                             3.在以上摆放方式中选出部队不在山上的

           接下来解释动态规划：

           状态dp[r][i]：第r行状态为state[i]时的最大值

           状态转移方程：dp[r][i]=max{dp[r-1][j]+dp[r-2][k]}+soldier[i]（满足ijk三种状态不相互攻击）

           可以改成三维形式：dp[r][i][j]=max{dp[r-1][j][k]}+soldier[i]表示第r行状态为state[i]，第r-1行状态为state[j]时的最大值

代码：转载自http://www.cnblogs.com/scau20110726/archive/2013/02/27/2935256.html

#include <cstdio>
#include <cstring>
#define MAXR 110 //行数
#define MAXC 15  //列数
#define MAXM 70  //状态数
#define max(a,b) a>b?a:b //返回较大值
#define CL(a) memset(a,0,sizeof(a)) //初始化清空数组
#define legal(a,b) a&b //判断两个状态共存时是否合法，合法为0，不合法为非0

int row,col;  //行列
int nums;  //仅是两个炮兵不互相攻击的条件下，符合条件的状态个数
int base[MAXR];  //第i行的原地图压缩成的一个状态
int state[MAXM]; //仅是两个炮兵不互相攻击的条件下，符合条件的状态（一个十进制数）
int soldier[MAXM]; //对应着，在state[i]状态下能放多少个士兵
int dp[MAXR][MAXM][MAXM];
//dp[i][j][k] 表示第i行状态为state[j]，第i-1行状态为state[k]时的最优解
char g[MAXR][MAXC];

int main()
{
    CL(base); CL(state); CL(soldier); CL(dp);
    nums=0;

    scanf("%d%d",&row,&col);
    for(int i=0; i<row; i++)  //先计算原始地图的状态数
    {
        scanf("%s",g[i]);
        for(int j=0; j<col; j++)
        if(g[i][j]=='H') base[i]+=1<<j; //像0110000，这里计算为6
    }

    for(int i=0; i<(1<<col); i++)  //仅是两个炮兵不互相攻击的条件下计算所有状态
    {
        if( legal(i,i<<1) || legal(i,i<<2)) continue; //i这个状态出现了士兵两两攻击
        int k=i;
        while(k) //这个循环计算状态i的二进制形式里面有多少个1，也就是放了多少个士兵
        {
            soldier[nums]+=k&1; //等价于k%2，相当于判断k的二进制形式里面有多少个1
            k=k>>1;
        }
        state[nums++]=i; //保存这个合法的状态
    }

    /***************************************/
    //for(int i=0; i<nums; i++) printf("%d %d\n",state[i],soldier[i]);
    /***************************************/

    for(int i=0; i<nums; i++) //先初始化dp[0][i][0]，即初始化第1行的情况
    {
        if(legal(state[i],base[0])) continue;
        //在state[i]的基础上，还要满足士兵不能放在山上，这个判断就是处理这个问题的
        dp[0][i][0]=soldier[i];
    }
    for(int i=0; i<nums; i++) //接着初始化dp[1][i][j]，即第2行的情况
    {
        if(legal(state[i],base[1])) continue;
        for(int j=0; j<nums; j++) //枚举第1行的状态
        {
            if(legal(state[j],base[0])) continue;
            if(legal(state[i],state[j])) continue;
            dp[1][i][j]=max(dp[1][i][j] , dp[0][j][0]+soldier[i]);
            //状态转移方程
        }
    }

    for(int r=2; r<row; r++) //第3行开始DP直到最后
        for(int i=0; i<nums; i++) //枚举第r行的状态
        {
            if(legal(state[i],base[r])) continue;
            for(int j=0; j<nums; j++) //枚举第r-1行的状态
            {
                if(legal(state[j],base[r-1])) continue;
                if(legal(state[i],state[j])) continue;
                //第r行的士兵和第r-1行的士兵相互攻击
                for(int k=0; k<nums; k++) //枚举第r-2行的状态
                {
                    if(legal(state[k],base[r-2])) continue;

                    if(legal(state[j],state[k])) continue;
                    //第r-1行的士兵和第r-2行的士兵相互攻击
                    if(legal(state[i],state[k])) continue;
                    //第r行的士兵和第r-1
                    dp[r][i][j]=max(dp[r][i][j] , dp[r-1][j][k]+soldier[i]);
                }
            }
        }

    int ans=0;
    for(int i=0; i<nums; i++)
        for(int j=0; j<nums; j++) //枚举dp[row-1][i][j]
            ans=max(ans,dp[row-1][i][j]);
    printf("%d\n",ans);

    return 0;
}