poj -1185 炮兵阵地 (经典状压dp)

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

参考博客:http://poj.org/problem?id=1185

大神博客已经讲的很清楚了,注意存状态的时候是从1开始的,所以初始化的时候也是dp[1][1][state],从0开始的话,状态就是dp[1][0][state]了.

dp[i][j][k]表示第i行状态为k第i-1行状态为j时的方案数.

dp[i][j][k]=max(dp[i][j][k],dp[i-1][t][j]+num[k]); (num[k]为k状态中1的个数)

边界条件：dp[1][1][i]=num[i],状态i可以满足第一行的条件。

还有就是为什么每一行最多只有60种状态,poj题目讨论里面有人给出了枚举的代码。

 #include <iostream>
 using namespace std;
 bool isok( int c ) {
     return !(c&(c<<)||c&(c<<));//同一行中不能有相邻的1距离小于3
 }
 int main() {
     int count=;
     for( int i=; i<; i++ )
         count += isok(i);
     cout<<count<<endl;
     return ;
 }

 #include <iostream>
 #include <cstdio>
 #include <cmath>
 #include <vector>
 #include <cstring>
 #include <string>
 #include <algorithm>
 #include <string>
 #include <set>
 #include <functional>
 #include <numeric>
 #include <sstream>
 #include <stack>
 //#include <map>
 #include <queue>
 #include <deque>
 //#pragma comment(linker, "/STACK:102400000,102400000")
 #define CL(arr, val)    memset(arr, val, sizeof(arr))

 #define ll long long
 #define INF 0x7f7f7f7f
 #define lc l,m,rt<<1
 #define rc m + 1,r,rt<<1|1
 #define pi acos(-1.0)

 #define L(x)    (x) << 1
 #define R(x)    (x) << 1 | 1
 #define MID(l, r)   (l + r) >> 1
 #define Min(x, y)   (x) < (y) ? (x) : (y)
 #define Max(x, y)   (x) < (y) ? (y) : (x)
 #define E(x)        (1 << (x))
 #define iabs(x)     (x) < 0 ? -(x) : (x)
 #define OUT(x)  printf("%I64d\n", x)
 #define lowbit(x)   (x)&(-x)
 #define Read()  freopen("a.txt", "r", stdin)
 #define Write() freopen("b.txt", "w", stdout);
 #define maxn 110
 #define maxv 5010
 #define mod 1000000000
 using namespace std;
 int n,m;
 char map[][],num[],top;
 int stk[],cur[];
 int dp[][][];

 inline bool ok(int x) //判断该状态是否合法,即同一行不存在相邻1之间的距离小于3的
 {
     if(x&(x<<)||x&(x<<)) return ;
     return ;
 }
 inline void jnite() //找到所有可能合法的状态
 {
     top=;
     int total=<<m;
     for(int i=;i<total;i++)
         if(ok(i)) stk[++top]=i;
 }

 inline bool fit(int x,int k) //判断状态x是否与第k行匹配
 {
     if(cur[k]&x) return ;
     return ;
 }
 inline int jcount(int x) //计算一个整型数x的二进制中1的个数(用于初始化)
 {
     int cnt=;
     while(x)
     {
         cnt++;
         x&=(x-); //很精炼,每次都会与掉一个1
     }
     return cnt;
 }
 int main()
 {
     //Read();
     while(~scanf("%d%d",&n,&m))
     {
         if(n==&&m==) break;
         jnite();
         for(int i=;i<=n;i++) scanf("%s",map[i]+);
         for(int i=;i<=n;i++)
         {
             cur[i]=;
             for(int j=;j<=m;j++)
             {
                 if(map[i][j]=='H') cur[i]+=(<<(j-));
             }
             //printf("%d\n",cur[i]);
         }
         memset(dp,,sizeof(dp));
         for(int i=;i<=top;i++) //初始化第一行
         {
             num[i]=jcount(stk[i]);
             //printf("%d\n",num[i]);
             if(fit(stk[i],)) dp[][][i]=num[i];
         }
         for(int i=;i<=n;i++) {
             for(int t=;t<=top;t++) {
                 if(!fit(stk[t],i)) continue;//第i行是否冲突
                 for(int j=;j<=top;j++) {
                     if(stk[t]&stk[j]) continue;//第i行和第i-2行是否冲突
                     for(int k=;k<=top;k++) {
                         if(stk[t]&stk[k]) continue;//第i行和第i-1行是否冲突
                         dp[i][k][t]=max(dp[i][k][t],dp[i-][j][k]+num[t]);
                       //  printf("%d\n",dp[i][k][t]);
                     }
                 }
             }
         }
         int ans=; //得到最大值
         for(int i=;i<=top;i++)
             for(int j=;j<=top;j++)
             ans=max(ans,dp[n][i][j]);
         printf("%d\n",ans);
     }
     return ;
 }