POJ 1185 炮兵阵地 (状压DP)

题目地址：POJ 1185

这题的一个判断条件写错了。。导致错了好多遍。。sad。。

这题跟3254差不多。

可以发现，对于每一行放大炮的状态，只与它上面一行和上上一行的状态有关，每一行用状态压缩的表示方法，0表示不放大炮，1表示放大炮，同样的，先要满足硬件条件，即有的地方不能放大炮，然后就是每一行中不能有两个1的距离小于2（保证横着不互相攻击），这些要预先处理一下。然后就是状态表示和转移的问题了，因为是和前两行的状态有关，所以要开个三维的数组来表示状态，当前行的状态可由前两行的状态转移而来。即如果当前行的状态符合前两行的约束条件（不和前两行的大炮互相攻击），则当前行的最大值就是上一个状态的值加上当前状态中1的个数（当前行放大炮的个数） 

代码如下：

#include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL __int64
#define pi acos(-1.0)
const int mod=100000000;
const int INF=0x3f3f3f3f;
const double eqs=0.01;
char s[110][12];
int dp[110][100][100], num[1100];
vector<int>vec[110];
int Cal(int x, int m)
{
        int i, ans=0;
        for(i=0; i<m; i++) {
                if(x&(1<<i))
                        ans++;
        }
        return ans;
}
void init(int n,int m, int tot)
{
        int i, j, flag, k;
        for(i=0; i<n; i++) {
                for(j=0; j<tot; j++) {
                        flag=0;
                        for(k=0; k<m; k++) {
                                if(j&(1<<k)) {
                                        if(s[i][k]=='H') {
                                                flag=1;
                                                break;
                                        }
                                }
                        }
                        if(!flag&&!(j&(j<<1))&&!(j&(j<<2))) {
                                vec[i].push_back(j);
                        }
                }
        }
        for(i=0; i<tot; i++) {
                num[i]=Cal(i,m);
        }
}
int main()
{
        int n, m, i, j, tot, k, h, max1;
        //("1.txt","r",stdin);
        //freopen("2.txt","w",stdout);
        while(scanf("%d%d",&n,&m)!=EOF) {
                for(i=0;i<n;i++){
                    vec[i].clear();
                }
                for(i=0; i<n; i++) {
                        scanf("%s",s[i]);
                }
                tot=1<<m;
                memset(dp,0,sizeof(dp));
                init(n,m,tot);
                max1=0;
                for(i=0; i<vec[0].size(); i++) {
                        max1=max(max1,num[vec[0][i]]);
                }
                for(i=1; i<n; i++) {
                        for(j=0; j<vec[i].size(); j++) {
                                for(k=0; k<vec[i-1].size(); k++) {
                                        if(vec[i][j]&vec[i-1][k]) continue ;
                                        if(i==1) {
                                                dp[1][j][k]=num[vec[0][k]]+num[vec[1][j]];
                                                max1=max(max1,dp[i][j][k]);
                                                continue ;
                                        }
                                        for(h=0; h<vec[i-2].size(); h++) {
                                                if((vec[i][j]&vec[i-2][h])||(vec[i-1][k]&vec[i-2][h])) continue ;
                                                dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][h]+num[vec[i][j]]);
                                                max1=max(max1,dp[i][j][k]);
                                        }
                                }
                        }
                }
                printf("%d\n",max1);
        }
        return 0;
}