hdu-4539 郑厂长系列故事――排兵布阵(状态压缩)

来源:互联网 发布:微软keeper是什么软件 编辑:程序博客网 时间:2024/05/22 04:39

类似于poj-1185

//曼哈顿距离 d = |x1 - x2| + |y1 - y2|//解决曼哈顿距离冲突 s[i]&(s[j] << 1)||s[i]&(s[j] >> 1)#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;int n,m,num;int s[200],p[200],a[200];int dp[110][222][222];bool ok(int x){    if(x&(x<<2)) return false;    return true;}void init(){    num = 0;    for(int i = 0;i < (1<<m);i++){        if(ok(i)) s[num++] = i;    }}int cal(int x){    int count = 0;    while(x){        count ++;        x &= (x-1);    }    return count;}bool find(int x,int y){    if(x&y) return false;    return true;}bool find1(int x,int y){    if((x&(y<<1))||(x&(y>>1))) return false;    return true;}int main(){    int temp;    while(~scanf("%d%d",&n,&m)){        memset(a,0,sizeof(a));        for(int i = 0;i < n;i++){            for(int j = 0;j < m;j++){                scanf("%d",&temp);                if(!temp) a[i] += (1<<j);            }        }        init();        memset(dp,0,sizeof(dp));        for(int i = 0;i < num;i++){            p[i] = cal(s[i]);            if(find(s[i],a[0])) dp[0][i][0] = p[i];        }        for(int i = 1;i < n;i++){            for(int j = 0;j < num;j++){                if(s[j]&a[i]) continue;                for(int k = 0;k < num;k++){                    if(!find1(s[j],s[k])) continue;                    for(int l = 0;l < num;l++){                        if(s[l]&s[j]) continue;                        if(!find1(s[j],s[k])) continue;                        dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][l] + p[j]);                    }                }            }        }        int ans = 0;        for(int i = 0;i < num;i++){            for(int j = 0;j < num;j++){                ans = max(ans,dp[n-1][i][j]);            }        }        printf("%d\n",ans);    }}


0 0
原创粉丝点击