HDU 4539 郑厂长系列故事——排兵布阵 (状压DP)

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

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

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

Source
2013腾讯编程马拉松复赛第二场（3月30日）

Recommend
liuyiding

普通姿势的状压没啥好说的

#include "cstring"
#include "cstdio"
#include "iostream"
#include "string.h"
#include "algorithm"
#include "cmath"
using namespace std;
int mmap[200][200];
int n,m;
int now_status[2000],now_sum[2000],nowid,preid_1,preid_2,dp[2000][2000],pre_status1[2000],pre_status2[2000] ,pre_sum[2000][2000];
void dfs(int row,int column,int status,int sum)
{
    if(column>=m)
    {
        now_status[++nowid]=status;
        now_sum[nowid]=sum;
        return;
    }
    if(column<=1&&mmap[row][column]==1)
    {
        dfs(row,column+1,status|1<<column,sum+1);
    }
    if(column>1&&mmap[row][column]==1&&!(status&(1<<(column-2))))
        dfs(row,column+1,status|1<<column,sum+1);
    dfs(row,column+1,status,sum);
    return;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                scanf("%d",&mmap[i][j]);
        pre_sum[0][0]=pre_status1[0]=pre_status2[0]=0;
        preid_1=preid_2=0;
        for(int k=0;k<n;k++)
        {
            nowid=-1;
            dfs(k,0,0,0);
            for(int i=0;i<=nowid;i++)
             for(int j=0;j<=preid_1;j++)
                 dp[i][j]=0;
            for(int i=0;i<=nowid;i++)
            {
                for(int j=0;j<=preid_1;j++)
                {
                    for(int q=0;q<=preid_2;q++)
                    {
                        if(now_status[i]&(pre_status1[j]<<1))
                            continue;
                        if(now_status[i]&(pre_status1[j]>>1))
                            continue;
                        if(now_status[i]&pre_status2[q])
                            continue;
                        dp[i][j]=max(dp[i][j],pre_sum[j][q]+now_sum[i]);
                    }
                }
            }
            for(int i=0;i<=nowid;i++)
                for(int j=0;j<=preid_1;j++)
                    pre_sum[i][j]=dp[i][j];
            for(int i=0;i<=preid_1;i++)
                pre_status2[i]=pre_status1[i];
            for(int i=0;i<=nowid;i++)
                pre_status1[i]=now_status[i];
            preid_2=preid_1;
            preid_1=nowid;
        }
        int ans=0;
        for(int i=0;i<=preid_1;i++)
            for(int j=0;j<=preid_2;j++)
                ans=max(ans,dp[i][j]);
        printf("%d\n",ans);
    }
    return 0;
}