炮兵阵地

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+状态压缩+位运算+滚动数组
状态压缩：将每一行进行状态压缩为一个(0~2^m-1)的数字，这样可摆放状态（即左两格和右两格都不
          摆放炮兵阵地），经估算大约有60个左右，再将山地的01编码转换成十进制数，用a数组记
          下每一行山地的状态。
滚动数组：由于dp数组需要开三维，如果开成dp[n][1<<m][1<<m]占用内存太大，有可能爆内存，又dp
           数组中dp[i][j][k] 只与 dp[i-1][k][g] 有关所以可以将数组简化成dp[2][1<<m][1<<m],
           状态转移时利用迭代变量roll进行滚动，当前状态roll 与 (roll+1)%2 有关。
这题写好了之后，改正细节错误后，还是wa，后来看了下数据范围，发现会有几组特殊数据，测了下，发现
有问题，即当n<2时，dp的初始化需要进行改动，把这个改好了之后，果断AC了！！！

Codeing:

# include<stdio.h>
# include<string.h>
int max(int x,int y)
{
    return (x>y?x:y);
}
int dp[2][1<<10][1<<10];
int main()
{
    char ch;
    int n,m,i,j,sc,k,g,roll;
    int a[105],s[65],c[65];
    while (~scanf("%d%d",&n,&m))
    {
        memset(a,0,sizeof(a));
        for (i=0;i<n;i++)
        {
            int k,ss;
            getchar();
            for (j=0;j<m;j++)
            {
                scanf("%c",&ch);
                if (ch=='P') k=0;
                else k=1;
                a[i]=a[i]*2+k;  //记录每行山地的状态
            }
        }
        sc=0;
        for (i=0;i<(1<<m);i++)
        {
            int k=i;
            if ((k<<1)&i||(k<<2)&i) continue;//判断该状态时候符合互不攻击的要求
            c[sc]=k%2;
            while (k=(k>>1)) c[sc]+=k%2;  //记录此状态下的炮兵阵地数量
            s[sc++]=i;                                 //记录符合题意该状态
        }
        memset(dp,-1,sizeof(dp));          
        if (n>1)
        for (i=0;i<sc;i++)
        {
            if (s[i]&a[1]) continue;
            for (j=0;j<sc;j++)
            {
                if ((s[i]&s[j])||(s[j]&a[0])) continue;
                dp[1][s[i]][s[j]]=c[i]+c[j];
            }
        }
        else if (n>0)
        {
            for (j=0;j<sc;j++)
            {
                if (s[j]&a[0]) continue;
                dp[1][s[j]][s[0]]=c[j];
            }
        }                                    //dp的初始化，特别注意n<2时的情况
        roll=1;                           //迭代变量
        for (i=2;i<n;i++)
        {
            roll=(roll+1)%2;
            for (j=0;j<sc;j++)
            if ((s[j]&a[i])==0)
            for (k=0;k<sc;k++)
            if ((s[j]&s[k])==0&&(s[k]&a[i-1])==0)
            for (g=0;g<sc;g++)
            {
                if ((s[j]&s[g])||(s[k]&s[g])||(s[g]&a[i-2])) continue;
                if (dp[roll][s[j]][s[k]]==-1) dp[roll][s[j]][s[k]]=dp[(roll+1)%2][s[k]][s[g]]+c[j];
                else dp[roll][s[j]][s[k]]=max(dp[roll][s[j]][s[k]],dp[(roll+1)%2][s[k]][s[g]]+c[j]);  //状态转移核心方程
            }
        }
        int sum=0;
        for (i=0;i<sc;i++)
        for (j=0;j<sc;j++)
        if (dp[roll][s[i]][s[j]]>sum) sum=dp[roll][s[i]][s[j]];
        printf("%d\n",sum);
    }
    return 0;
}

                   预祝AC！！！