B - 炮兵阵地 POJ - 1185

 刚开始接触状压dp， 感觉就形式上来说还是比较固定的，和数位dp差不多，但是各种位运算还是不够熟练，用起来有点不是得心应手

还是得多练练，这部分代码还不是很好写。

 

题意: 中文题， 不多说了。

思路： dp[n][prestate][preprestate] 表示第n行的时候的状态是由前两行的状态推倒出的 

转移方程: dp[i][j][k] = max(dp[i][j][k], dp[i-1][k][t] + num[j]);

接下来就是看代码了，里面注释都有

#include <string.h>
#include <algorithm>
#include <stdio.h>
#include <string>
#include <iostream>
using namespace std;

int dp[101][1025][1025];
int a[1025];
int n, m;
int vnum;
int num[1025], vst[1025];

int getsum(int i) {    // __builtin_popcount(i)
    int res = 0;
    while (i) {
        i &= i - 1;
        ++res;
    }
    return res;
}

void initialState() {//vst数组记录的是对于每一行符合题目要求的安防方法 num记录的是对应安防方法的安防炮兵个数
    vnum = 0;
    for (int i = 0; i < (1 << m); ++i) {
        if (!(i&(i << 1)) && !(i&(i << 2))) {            
            vst[vnum] = i;
            num[vnum++] = getsum(i);
            cout << vst[vnum - 1] << ends << num[vnum - 1] << endl;
        }
    }
}

int main() {
    ios::sync_with_stdio(false);
    while (cin >> n >> m) {
        for (int i = 0; i < n; ++i) {
            a[i] = 0;
            for (int j = 1; j <= m; ++j) {
                char c; cin >> c;
                a[i] = a[i] * 2 + (c == 'P');
            }
        }
        for (int i = 0; i <= n; ++i)
            for (int j = 0; j <= 61; ++j)
                for (int k = 0; k <= 61; ++k)
                    dp[i][j][k] = -1;
        initialState();
        for (int i = 0; i < vnum; ++i)
            if (!(vst[i] & (~a[0])))
                dp[0][i][0] = num[i];
        for (int i = 1; i < n; ++i) {
            for (int j = 0; j < vnum; ++j) {
                if ((~a[i]) & vst[j]) continue;            //判断当前行vit[j]中1的集合是不是属于a[i]的集合
                for (int k = 0; k < vnum; ++k) {
                    if ((~a[i - 1])&vst[k]) continue;    //对于上一行同样的判断
                    if (vst[j] & vst[k]) continue;        //上下相邻两行得符合题目要求
                    for (int t = 0; t < vnum; ++t) {
                        if (vst[j] & vst[t]) continue;
                        if (-1 == dp[i - 1][k][t]) continue;
                        dp[i][j][k] = max(dp[i][j][k], dp[i-1][k][t] + num[j]); //状态转移方程
                    }
                }
            }
        }
        int res = 0;
        for (int i = 0; i < vnum; ++i) {
            for (int j = 0; j < vnum; ++j) {
                if (dp[n - 1][i][j] > res)
                    res = dp[n - 1][i][j];
            }
        }
        cout << res << endl;
    }
    return 0;
}