Poj1038 Bugs Integrated, Inc.【状态压缩DP】

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

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

要知道能否在第k-2行至第k行这三行中切割出芯片及切割出多少芯片，需要知道这3行上每个方块的具体状态：是否损坏，是否已经被其他芯片使用，这两种情况都称为不可用。那么可以用一个三进制数来表征某一列第k-1行至第k行的状态：0表示两个方块都可用，1表示仅有第k行的方块可用，2表示第k行的方块不可用。用dp[k][t]表示第k-1行至第k行所有列的三进制特征数t的每一位的状态最多能够切割出的芯片数。那么在状态转移时，凭借上一层的状态和本层的状态就能知道这3行的全部情况。

该题中父状态和子状态是多对多的关系，应该对每个可行的子状态（有些状态是不可能形成的），求解父状态并更新它的最大值。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <cmath>
#include <queue>
using namespace std;
template <class T> void checkmin(T &t,T x) {if(x < t) t = x;}
template <class T> void checkmax(T &t,T x) {if(x > t) t = x;}
template <class T> void _checkmin(T &t,T x) {if(t==-1) t = x; if(x < t) t = x;}
template <class T> void _checkmax(T &t,T x) {if(t==-1) t = x; if(x > t) t = x;}
typedef pair <int,int> PII;
typedef pair <double,double> PDD;
typedef long long ll;
#define foreach(it,v) for(__typeof((v).begin()) it = (v).begin(); it != (v).end ; it ++)
int mp[155][15];
int tri[12] = {0,1,3,9,27,81,243,729,2187,6561,19683,59049};
int dp[2][59049];  //用滚动数组记录每一行的状态
int pre[15],now[15]; //上一行与本行的三进制状态
int n , m;
int getten(int a[]) {
    //把三进制状态转化为十进制数
    int sum = 0;
    for(int i=1;i<=m;i++) sum += tri[i] * a[i];
    return sum;
}
void gettri(int k,int a[]) {
    //把十进制数转化为三进制数组
    for(int i=1;i<=m;i++) {
        a[i] = k % 3;
        k /= 3;
    }
    return;
}
void dfs(int i,int j,int last,int key) {
    //四个参数：i当前行，j当前列，当前状态值，当前状态
    int k;
    checkmax(dp[i%2][key],last);
    if(j >= m) return;
    if(!pre[j]&&!pre[j+1]&&!now[j]&&!now[j+1]) { //竖着切割
        now[j] = now[j+1] = 2;
        k = getten(now);
        dfs(i,j+2,last+1,k);
        now[j] = now[j+1] = 0;
    }
    if(j<m-1&&!now[j]&&!now[j+1]&&!now[j+2]) { //横着切割
        now[j] = now[j+1] = now[j+2] = 2;
        k = getten(now);
        dfs(i,j+3,last+1,k);
        now[j] = now[j+1] = now[j+2] = 0;
    }
    dfs(i,j+1,last,key); //不作改动，搜索下一列
    return;
}
int main() {
    int T , k;
    scanf("%d",&T);
    while(T--) {
        scanf("%d%d%d",&n,&m,&k);
        for(int i=0;i<tri[m+1];i++) dp[1][i] = -1; //先初始化第1行的所有状态都是不可能形成的
        for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) mp[i][j] = 0;
        for(int i=0;i<k;i++) {
            int a , b;
            scanf("%d%d",&a,&b);
            mp[a][b] = 1;
        }
        for(int i=1;i<=m;i++) pre[i]=mp[1][i]+1; //计算第1行的状态，第0行的所有方块视为不可用
        int tmp = getten(pre);
        dp[1][tmp] = 0; //设置第1行的原本状态是可以形成的，状态值为0
        for(int i=2;i<=n;i++) {          //枚举子状态
            for(int j=0;j<tri[m+1];j++)   //先设置父状态值为-1，等待更新
                dp[i%2][j] = -1;
            for(int j=0;j<tri[m+1];j++) {
                if(dp[(i+1)%2][j]==-1) continue; //跳过不可能的子状态
                gettri(j,pre);
                for(int t=1;t<=m;t++) {
                    if(mp[i][t]) now[t] = 2;
                    else now[t] = max(pre[t]-1,0);
                } //根据上一行的状态得到本行的状态
                tmp = getten(now);
                dfs(i,1,dp[(i+1)%2][j],tmp);
                //神搜这一层所有可能的状态
            }
        }
        int ans = 0;
        for(int i=0;i<tri[m+1];i++)
            checkmax(ans , dp[n%2][i]);
        printf("%d\n" , ans);
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/aiiYuu/archive/2013/04/05/3000463.html