POJ_1185_状态压缩dp

http://poj.org/problem?id=1185

一次考虑两行，比一行略为复杂。sta保存每种状态炮兵位置，sum保存每种状态当行炮兵总数，a保存地形，dp[i][j][k]表示到第i行当前行j状态上一行k状态的最大炮兵数。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
//int main()
//{
//    int num = 0;
//    for(int i = 0;i <= (1<<10);i++)
//    {
//        if(i & (i<<2) || i & (i<<1))  continue;
//        num++;
//    }
//    printf("%d",num);
//}

int dp[][][],a[] = {},n,m,sta[],sum[] = {};

int main()
{
    memset(dp,-,sizeof(dp));
    scanf("%d%d",&n,&m);
    getchar();
    if(n ==  && m == )
    {
        printf("0\n");
        return ;
    }
    for(int i = ;i < n;i++)
    {
        char temp[];
        gets(temp);
        for(int j = ;j < m;j++)
        {
            if(temp[j] == 'H')  a[i] |= (<<j);
        }
    }
//    for(int i = 0;i < n;i++)    printf("%d\n",a[i]);
    int cnt = ;
    for(int i = ;i < (<<m);i++)
    {
        if(i & (i<<))  continue;
        if(i & (i<<))  continue;
        sta[cnt] = i;
        int temp = i;
        while(temp)
        {
            sum[cnt] += temp%;
            temp /= ;
        }
        cnt++;
    }
//    for(int i = 0;i < cnt;i++)    printf("%d %d\n",sta[i],sum[i]);
    for(int i = ;i < cnt;i++)
    {
        if(a[] & sta[i])   continue;
        dp[][i][] = sum[i];
    }
//    for(int i = 0;i < cnt;i++)    printf("%d\n",dp[0][i][0]);
    for(int i = ;i < cnt;i++)
    {
        if(a[] & sta[i])   continue;
        for(int j = ;j < cnt;j++)
        {
            if(sta[i] & sta[j]) continue;
            dp[][i][j] = max(dp[][i][j],dp[][j][]+sum[i]);
        }
    }
    for(int i = ;i < n;i++)
    {
        for(int j = ;j < cnt;j++)
        {
            if(a[i] & sta[j])   continue;
            for(int k = ;k < cnt;k++)
            {
                if(sta[j] & sta[k])   continue;
                for(int l = ;l < cnt;l++)
                {
                    if(sta[j] & sta[l]) continue;
                    if(dp[i-][k][l] == -) continue;
                    dp[i][j][k] = max(dp[i][j][k],dp[i-][k][l]+sum[j]);
                }
            }
        }
    }
//    for(int i = 0;i < cnt;i++)
//    {
//        for(int j = 0;j < cnt;j++)  printf("%d ",dp[n-1][i][j]);
//        printf("\n");
//    }
    int ans = -;
    for(int i = ;i < cnt;i++)
    {
        for(int j = ;j < cnt;j++)  ans = max(ans,dp[n-][i][j]);
    }
    printf("%d\n",ans);
    return ;
}