状态压缩dp_POJ1185炮兵阵地

    位压缩dp，单行有效状态最多只有60种。

    0.压缩：用j的二进制表示单行状态（0<=j<=2^m-1)。j的二进制表示中，0的数位表示该位置不放炮兵，1的数位表示该位置放置了一个炮兵。

    1.状态表示：dp[i][j][k]表示第i行放置方式为j，第i-1行放置方式为k时所能得到的最大炮兵数量

    2.状态转移：dp[i][j][k] = max{ dp[i-1][k][s]} + num[j] }

       {

        j，k，s均为有效状态且分别属于map[i],map[i-1],map[i-2]的子集且j,k,s互不相交。num[j]是j的二进制表示中1的个数。

       }

      j为有效状态 iff (j&(j>>1)) == 0 && (j&(j>>2)) == 0 ;（炮兵间距离不能小于2）

      j为map[i]的子集 iff (map[i] & j) == j ; (map[i]中平原‘P’才可以布置炮兵）

      j与k不相交 iff （j & k）== 0

    3.初始化：dp = 0;

      map[0] = (1<<m)-1;

      dp[1][j][0] = 1 (j是有效状态且是map[1]的子集)。

    4.Answer：∑dp[n][j][k]。

//M<=10单行最多只有60种状态，但是M==11时，状态暴增至1075种！。

//先预处理出单行所有可能状态，再dp。

//细节：state[]最后添加一个极大值，防止state[i]<=Max时越界。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <vector>
#include <cstdlib>
#define  clr(A) memset(A,0,sizeof(A))
#define LL long long

using namespace std;

int const maxn = 3000;
int const D = 70;
int state[D],sum;
int num[D];
int dp[102][D][D];
int map[102];
int hash[D][D];
int n,m;

void Getstate()
{
    sum = 0;
    for(int i = 0;i<=(1<<10)-1;i++)
    {
        if ((i&(i>>1)) != 0 || (i&(i>>2)) != 0) continue;
        state[sum] = i;
        num[sum] = 0;
        int x = i;
        while(x)
        {
            num[sum] += (x&1);
            x >>= 1;
        }
        sum++;
    }
    state[sum++] = (1<<20);
}

void Dp()
{
    int Max = (1<<m)-1;

    Getstate();

    clr(hash);
    for(int i = 0;state[i]<=Max;i++)
    for(int j = 0;state[j]<=Max;j++)
    if((state[i]&state[j]) == 0) hash[i][j] = hash[j][i] = 1;

    clr(dp);

    for(int i = 0;state[i]<=Max;i++)
    if((map[1]&state[i]) == state[i])
    dp[1][i][0] = num[i];
    map[0] = Max;

    for(int i = 2;i<=n;i++)
    for(int j = 0;state[j]<=Max;j++)
    if((map[i]&state[j]) == state[j])
    for(int k = 0;state[k]<=Max;k++)
    if((map[i-1]&state[k]) == state[k] && hash[j][k])
    for(int s = 0;state[s]<=Max;s++)
    if((map[i-2]&state[s]) == state[s] && hash[j][s]&&hash[k][s])
    {
        dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][s]+num[j]);
    }
}

int main()
{

    while(~scanf("%d%d\n",&n,&m))
    {
        clr(map);
        char st[20];
        for(int i = 1; i<=n;i++)
        {
            scanf("%s",st);
            for(int j = 0; st[j]!='\0';j++)
            map[i] = map[i]*2+(st[j]=='P');
        }

        Dp();
        int ans = 0;
        for(int i = 0;i<sum;i++)
        for(int j = 0;j<sum;j++)
        ans = max(ans,dp[n][i][j]);
        printf("%d\n",ans);
    }
    return 0;
}