poj 1038 Bugs Integrated, Inc. __dp状态压缩

题目简述：

        一个生产芯片的工厂要生产一种规格为2*3的芯片，方法是先生产一块规格为n*m的矩形硅片，由n*m个正方形方块组成，但是硅片上存在一些损坏的方块，显示为黑色，它们不能用来制作芯片。现在给出硅片上每个损坏方块的位置，求用该硅片最多能切割出多少块芯片。

思路：

       为什么要用3进制来状态压缩呢，其实可以看成两个二进制结合而来的，如果题目改成只能横着放那就可以用二进制表示第i行第j个可用或不可用，如果题目改成只能竖着放那就可以用二进制表示第i行第j个和第i-1行第j个都可用或都不可用；把这两者结合起来就是一个三进制的数了。

    接下来看代码：

#include <iostream>#include <cstring>#include <cstdio>#include <cmath>#include <algorithm>const int maxn = 152;const int tri[11] = {1 ,3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049};using namespace std;int dp[2][59049], pre[12], now[12], mp[maxn][12];int v = 0, n, m;int getten(int* a) {    int result = 0, i;    for (i = 0; i < m; i++)        result += a[i] * tri[i];    return result;}void getthree(int k, int* a) {    int i;    for(i = 0; i < m; i++) {        a[i] = k % 3;        k /= 3;    }}void dfs(int j, int k, int ss) {    dp[v^1][k] = max(dp[v^1][k], ss);    if (j >= m - 1) return;    int kk;    if (!pre[j] && !pre[j+1] && !now[j] && !now[j+1]) {        now[j] = 2;        now[j+1] = 2;        kk = getten(now);        dfs(j+2, kk, ss + 1);        now[j] = 0;        now[j+1] = 0;    }    if (j < m - 2 && !now[j] && !now[j+1] && !now[j+2]) {        now[j] = 2;        now[j+1] = 2;        now[j+2] = 2;        kk = getten(now);        dfs(j+3, kk, ss + 1);        now[j] = 0;        now[j+1] = 0;        now[j+2] = 0;    }    dfs(j+1, k, ss); // 这里很重要，意思是跳过这一块；}int main() {    int k, D, i, j, t;    scanf("%d", &D);    while (D--) {        memset(mp, 0, sizeof(mp));        memset(dp[v], -1, sizeof(dp[v]));        scanf("%d%d%d", &n, &m, &k);        int x, y;        while (k--) {            scanf("%d%d", &x, &y);            mp[x-1][y-1] = 1;        }        for(i = 0; i < m; i++)            pre[i] = mp[0][i] + 1;        int tmp = getten(pre);        dp[v][tmp] = 0;//初始状态是可行的；        for(i = 1; i < n ; i++) {            memset(dp[v^1], -1, sizeof(dp[v^1]));            for(j = 0; j < tri[m]; j++) {                if (dp[v][j] == -1) continue;//这个状态没有；                getthree(j, pre);                for(t = 0; t < m; t++)                    if (mp[i][t]) now[t] = 2;                    else now[t] = max(0, pre[t] - 1);                dfs(0, getten(now), dp[v][j]);            }            v ^= 1;        }        int result = 0;        for(i = 0; i < tri[m]; i++)            if (result < dp[v][i])                result = dp[v][i];        cout<<result<<endl;    }    return 0;}