码迷,mamicode.com
首页 > 其他好文 > 详细

poj1185炮兵阵地状压dp

时间:2014-08-13 14:23:46      阅读:140      评论:0      收藏:0      [点我收藏+]

标签:style   blog   color   os   io   for   ar   art   

  压前两行的状态很容易想到,但是 直接搞  (1<<10) * (1<<10)  空间时间都明显受不了, 但是经过高人指点,你会发现:枚举每一行可行的状态,其实并不多,预先打表处理,不用 1->(1<<10)枚举每一种状态。。

然后记忆化搜就ok了。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <climits>
#include <string>
#include <iostream>
#include <map>
#include <cstdlib>
#include <list>
#include <set>
#include <queue>
#include <stack>
using namespace std;

int deal(int a,int b)
{
    return (a&(1<<b));
}
int m,n;
int chart[1000];
int vis[111];
int len;
int Map[111];
int dp[111][222][222];
int Max(int a,int b)
{
    return a>b? a:b;
}

int Min(int a,int b)
{
    return a>b?b:a;
}
int judge(int x)
{
    for(int i=0;i<m;i++)if(deal(x,i)){
        int t= Max(0,i-2);int t1=Min(m-1,i+2);
        for(int j=t;j<=t1;j++){
        if(i==j) continue;
        if(deal(x,j)) return 0;
        }
    }
    return 1;
}
void Init()
{
    int gg=1<<m;
    for(int i=0;i<gg;i++){
        if(judge(i))chart[len++]=i;
    }
}

int panduan(int x,int yici,int erci,int mat)
{
    if(yici&x) return 0;
    if(erci&x) return 0;
    if(x&(~mat)) return 0;
   // for(int i=0;i<m;i++) if(deal(x,i)) if(!deal(mat,i)) return 0;
    return 1;
}
int gao(int x,int yici,int erci)
{
    if(~dp[x][yici][erci]) return dp[x][yici][erci];
    if(x==n) return dp[x][yici][erci]= 0;
    int ret=0;
    for(int i=0;i<len;i++){
        if(panduan(chart[i],chart[yici],chart[erci],Map[x])) ret=Max(ret,gao(x+1,i,yici)+vis[i]);
    }
    return dp[x][yici][erci]= ret;
}
void Init1()
{
    memset(vis,0,sizeof(vis));
    for(int i=0;i<len;i++){
        int tem=chart[i];int ans=0;
        for(int j=0;j<m;j++){
            if(tem&1) ans++;
            tem>>=1;
        }
        vis[i]=ans;
    }
}
int main()
{
    char str[1000];
    while(cin>>n>>m){
        len=0;
        Init();Init1();
        memset(dp,-1,sizeof(dp));
        memset(Map,0,sizeof(Map));
        for(int i=0;i<n;i++){
            scanf("%s",str);
            for(int j=0;j<m;j++){
                if(str[j]==P)Map[i]|=(1<<j);
            }
        };
        cout<<gao(0,0,0)<<endl;
    }
    return 0;
}

 

poj1185炮兵阵地状压dp,布布扣,bubuko.com

poj1185炮兵阵地状压dp

标签:style   blog   color   os   io   for   ar   art   

原文地址:http://www.cnblogs.com/yigexigua/p/3909779.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!