hdu4539之状态压缩dp

郑厂长系列故事——排兵布阵

Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 1524    Accepted Submission(s): 564

Problem Description

　　郑厂长不是正厂长
　　也不是副厂长
　　他根本就不是厂长
　　事实上
　　他是带兵打仗的团长

　　一天，郑厂长带着他的军队来到了一个n*m的平原准备布阵。
　　根据以往的战斗经验，每个士兵可以攻击到并且只能攻击到与之曼哈顿距离为2的位置以及士兵本身所在的位置。当然，一个士兵不能站在另外一个士兵所能攻击到的位置，同时因为地形的原因平原上也不是每一个位置都可以安排士兵。
　　现在，已知n,m 以及平原阵地的具体地形，请你帮助郑厂长计算该阵地,最多能安排多少个士兵。

 

Input

输入包含多组测试数据；
每组数据的第一行包含2个整数n和m (n <= 100, m <= 10 )，之间用空格隔开；
接下来的n行，每行m个数，表示n*m的矩形阵地，其中1表示该位置可以安排士兵，0表示该地形不允许安排士兵。

 

Output

请为每组数据计算并输出最多能安排的士兵数量，每组数据输出一行。

 

Sample Input

   
   
   
   

    
    
    
    6 6
0 0 0 0 0 0
0 0 0 0 0 0
0 0 1 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
   
   
   
   

 

Sample Output

   
   
   
   

    
    
    
    2
   
   
   
   
   
   
   
   


   
   
   
   

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <queue>
#include <algorithm>
#include <map>
#include <cmath>
#include <iomanip>
#define INF 99999999
typedef long long LL;
using namespace std;

const int MAX=169+10;
int n,m,lastsize,nowsize,lastlastsize;
int Map[110][15],last[MAX],lastlast[MAX],ans[MAX];
int dp[MAX][MAX],temp[MAX][MAX],now[MAX];

void dfs(int id,int k,int p,int sum){
    if(k>=m){now[++nowsize]=p;ans[nowsize]=sum;return;}
    if(k>=2 && Map[id][k] == 1 && (p & (1<<(k-2))) == 0)dfs(id,k+1,p|(1<<k),sum+1);
    else if(k<2 && Map[id][k] == 1)dfs(id,k+1,p|(1<<k),sum+1);
    dfs(id,k+1,p,sum);
}

void DP(){
    temp[1][1]=last[1]=lastlast[1]=0;
    lastsize=nowsize=lastlastsize=1;
    for(int k=0;k<n;++k){
        nowsize=0;
        dfs(k,0,0,0);
        for(int i=1;i<=nowsize;++i)for(int j=1;j<=lastsize;++j)dp[i][j]=0;
        for(int i=1;i<=nowsize;++i){
            for(int j=1;j<=lastsize;++j){
                for(int t=1;t<=lastlastsize;++t){
                    if(now[i] & lastlast[t])continue;
                    if(now[i] & (last[j]<<1))continue;
                    if(now[i] & (last[j]>>1))continue;
                    dp[i][j]=max(dp[i][j],temp[j][t]+ans[i]);
                }
            }
        }
        for(int i=1;i<=nowsize;++i)for(int j=1;j<=lastsize;++j)temp[i][j]=dp[i][j];
        for(int i=1;i<=lastsize;++i)lastlast[i]=last[i];
        for(int i=1;i<=nowsize;++i)last[i]=now[i];
        lastlastsize=lastsize,lastsize=nowsize;
    }
}

int main(){
	while(~scanf("%d%d",&n,&m)){
		for(int i=0;i<n;++i){
			for(int j=0;j<m;++j)scanf("%d",&Map[i][j]);
		}
		DP();
		int sum=0;
		for(int i=1;i<=lastsize;++i)for(int j=1;j<=lastlastsize;++j)sum=max(sum,temp[i][j]);
		cout<<sum<<endl;
	}
	return 0;
}