POJ - 1185 跑兵阵地 状态压缩DP

    状态压缩DP...很简单了...一行最多10个...用1代表放炮兵..0代表不放...每行的状态最多1024种..并且还要合法..这样一来..一行的状态最多60种了...

    能影响到当前行的只有上两行..所以用三维的DP...dp[t][x1][x2]....t代表哪一层了...x1代表上一层的状态..x2代表当前层的状态...

    dp[t][x1][x2]= max ( dp[t-1][y][x1] + cnt[x2]) ...y是t-2层的状态...cnt是当前层状态的炮兵数,也就是1的数量...这里看出约束条件..就是状态x1,x2,y不能冲突...

Program:

#include<iostream>#include<cmath>#include<stack>#include<queue>#include<set>#include<algorithm>#include<stdio.h>#include<string.h>#define ll long long#define oo 1000000007using namespace std;   int n,m,a[105],anum;char arc[105][12];int cnt[65],dp[105][65][65];bool ok1(int x)   // 判断当前状态是否合法{       int p,i;       p=-10;        while (x)       {             i++;             if (x%2)             {                   if (i-p<=2) return false;                   p=i;             }             x/=2;       }              return true;}bool ok2(int t,int x)  //判断当前状态能否放在当前的图上{       int i;        for (i=m;i>=1 && x;i--)       {             if (x%2 && arc[t][i]=='H') return false;             x/=2;       }       return true;}bool ok3(int y,int x) //判断状态y在上,状态x在下,能否不冲突{       while (y && x)       {             if (y%2 && x%2) return false;             y/=2;  x/=2;       }       return true;}int main(){              int i,t,j,x,ans;       while (~scanf("%d%d",&n,&m))       {                anum=0;               a[0]=0;               for (i=0;i<(1<<m);i++)                  if (ok1(i))                   {                          a[++anum]=i;                          cnt[anum]=0;                          j=i;                          while (j)                          {                                 j=j&(j-1);                                 cnt[anum]++;                          }                  }               for (i=0;i<=n;i++) gets(arc[i]+1);               memset(dp,0,sizeof(dp));               for (i=1;i<=anum;i++)                 if (ok2(1,a[i]))  dp[1][0][i]=cnt[i];                 for (t=2;t<=n;t++)                  for (i=0;i<=anum;i++)                    if (ok2(t,a[i]))                      for (j=0;j<=anum;j++)                        if (ok2(t-1,a[j]) && ok3(a[j],a[i]))                          for (x=0;x<=anum;x++)                            if (dp[t][j][i]<dp[t-1][x][j]+cnt[i] && ok3(a[x],a[i]))                                dp[t][j][i]=dp[t-1][x][j]+cnt[i];                                        ans=0;               for (i=0;i<=anum;i++)                 for (j=0;j<=anum;j++)                   if (ans<dp[n][i][j])                        ans=dp[n][i][j];               printf("%d\n",ans);       }        return 0;}