poj 1185 状态压缩dp

题目链接：http://poj.org/problem?id=1185

第i行的状态由第I-1行和第i-2行决定

dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][p]+cou[j]); 

初始化第i=0行就够了，因为不会访问i-2

#include <cstdio>//改成c++就过原来用c的，wal

<span style="font-family: Arial, Helvetica, sans-serif;">#include <cstring></span>

int n,m;
int sta[60],cou[100]; //60是因为满足题目条件的，即相隔2个单位的，同时小于1<<10的，刚好只有60个，可以自己打打看，
int dp[110][60][60];
int map[110];
char str[12];

int max(int a,int b)
{
    if(a > b)
        return a;
    return b;
}

int main()
{
    int i,j,k,p;
    int sp,tmp;

    while(~scanf("%d %d",&n,&m))
    {
        memset(map,0,sizeof(map));
        memset(dp,-1,sizeof(dp));
        memset(cou,0,sizeof(cou));

        for(i = 0; i < n; i++)
        {
            scanf("%s",str);
            for(j = 0; j < m; j++)
            {
                if(str[j] == 'H')
                    map[i]+=(1<<j);
            }
        }

        sp = 0;
        tmp = (1<<m);
        for(i = 0; i < tmp; i++)
        {
            if(!(i&(i<<2)) && !(i&(i<<1)))
            {
                p = i;
                while(p)
                {
                    cou[sp]++;
                    p&=(p-1);
                }
                sta[sp++] = i;
            }
        }

        for(i = 0; i < sp; i++)
        {
            if(!(map[0]&sta[i]))
            {
                dp[0][i][0] = cou[i];      
            }
        }

        for(i = 1; i < n; i++)
        {
            for(j = 0; j < sp; j++)
            {
                if(!(map[i]&sta[j]))
                {
                    for(k = 0; k < sp; k++)
                    {
                        if(!(map[i-1]&sta[k]) && !(sta[j]&sta[k]))
                        {
                            for(p = 0; p < sp; p++)
                            {
                                if(!(sta[p]&sta[j]) && !(sta[p]&sta[k]) && dp[i-1][k][p]!=-1)
                                {
                                    dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][p]+cou[j]);
                                }
                            }
                        }
                    }
                }
            }
        }
        tmp = 0;
        for(i = 0; i < n; i++)
        {
            for(j = 0; j < sp; j++)
            {
                for(k = 0; k < sp; k++)
                {
                    tmp = max(tmp,dp[i][j][k]);
                }
            }
        }
        printf("%d\n",tmp);
    }
    return 0;
}