炮兵阵地(POJ 1185状压dp)

题意:n*m地图‘H‘能放‘p‘不能放,布兵的方格上下左右不能布兵,给你地图求最大布兵数

分析:关系到前两行,所以dp[i][j][k]第i行状态为j,i-1行状态为k时的最大布兵数, 先求出所有可行的状态,统计出其布兵数。

#include <map>
#include <set>
#include <list>
#include <cmath>
#include <queue>
#include <stack>
#include <cstdio>
#include <vector>
#include <string>
#include <cctype>
#include <complex>
#include <cassert>
#include <utility>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
typedef pair<int,int> PII;
typedef long long ll;
#define lson l,m,rt<<1
#define pi acos(-1.0)
#define rson m+1,r,rt<<11
#define All 1,N,1
#define read freopen("in.txt", "r", stdin)
const ll  INFll = 0x3f3f3f3f3f3f3f3fLL;
const int INF= 0x7ffffff;
const int mod =  1000000007;
int n,m,a[200],dp[110][200][200];
int num,sum[2000],cas[2000];
bool judge(int x){
    return ((x<<1)&x)||((x<<2)&x)||((x>>1)&x)||((x>>2)&x);
}
int countnum(int x){
int t=0;
    while(x){
        if(x&1)t++;
    x>>=1;
    }
    return t;
}
//求预处理所有可能的状态
void init(){
    num=0;
    for(int i=0;i<(1<<m);++i){
        if(!judge(i)){
            cas[num]=i;
            sum[num++]=countnum(i);
        }
    }
}
void solve(){
    init();
    memset(dp,0,sizeof(dp));
    for(int i=0;i<num;++i){
        if(a[0]&cas[i])continue;
        dp[0][i][0]=sum[i];
    }
    for(int i=1;i<n;++i){
        //枚举当前状态
        for(int j=0;j<num;++j){
            if(cas[j]&a[i])continue;
        //枚举符合条件的上一行状态
            for(int k=0;k<num;++k){
                if(cas[k]&a[i-1])continue;
                if(cas[k]&cas[j])continue;
            int maxv=-1;
            for(int l=0;l<num;++l){
                if((cas[l]&cas[k])||(cas[l]&cas[j]))continue;
                maxv=max(maxv,dp[i-1][k][l]);
            }
            dp[i][j][k]=maxv+sum[j];
            }
        }
    }
    int maxn=-1;
    for(int i=0;i<num;++i)
    for(int j=0;j<num;++j){
        maxn=max(maxn,dp[n-1][i][j]);
    }
    printf("%d\n",maxn);
}
int main()
{
    char ch;
    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){
            cin>>ch;
            if(ch==‘H‘){
                a[i]|=(1<<j);
            }
        }
        solve();
    }
return 0;
}

  

时间: 2024-08-07 19:36:06

炮兵阵地(POJ 1185状压dp)的相关文章

poj 1185 状压dp+优化

http://poj.org/problem?id=1185 炮兵阵地 Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 29176   Accepted: 11303 Description 司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队.一个N*M的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用"P"表示),如下图.在每一格平原地形上最多可以布置一支

POJ 1185 状压DP

legal[] 保存所有在当前行可显示的状态,由dfs得到,len[]保存legal[]对应下标状态中的 1 的个数 , 也就是放置炮台的个数 state[i] 表示第 i 行这块区域的土地情况,H表示 1 ,P表示 0 那么每次加入一个legal状态  都要符合 !(legal[i] & state[k]) 1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 using namesp

POJ 1038 状压DP

一个公司生产一种2*3规模的芯片,但是原材料上面有一些地方是不能用来当作芯片材料的,给出原料大小,及上面不能做原料的点,问你怎么分解,可以使生成芯片最大化. 对M进行三进制状压 last数组存储第i-1行和i-2行状态,cur数组存储i行和i-1行状态 cur[k]=2; // 本行k位置和上行k位置都不可用 cur[k]=1; // 本行k位置可用,上行k位置不可用 cur[k]=0; // 本行和上行位置k均可用 必须用滚动数组,否则爆内存 #include "stdio.h" #

hdu 1185 状压dp 好题 (当前状态与上两行有关系)

/* 状压dp 刚开始&写成&&看了好长时间T0T. 状态转移方程 dp[i][k][j]=Max(dp[i][k][j],dp[i-1][l][k]+num[i][j]);(第i行的第j个状态有上一行的第k个状态得到) num[i][j]有两个功能,第一:判断第i行第j个状态是否合法 第二:判断第i行第j个状态的数目 */ #include<stdio.h> #include<string.h> #define N 110 int dp[N][N][N];

POJ 3254 (状压DP) Corn Fields

基础的状压DP,因为是将状态压缩到一个整数中,所以会涉及到很多比较巧妙的位运算. 我们可以先把输入中每行的01压缩成一个整数. 判断一个状态是否有相邻1: 如果 x & (x << 1) 非0,说明有相邻的两个1 判断一个状态能否放在该行: 如果 (a[i] & state) != state,说明是不能放置的.因为a[i]中存在某个0和对应state中的1,与运算之后改变了state的值 判断相邻两行的状态是否有同一列相邻的1: 如果(state & _state)不

Best Sequence(poj 1699) 状压dp(TSP)

类似于前两天做的那个wordstack.状压的其实有时候爆搜+记忆化也差不多. 就是这个是要与之前的都重合,移位预处理要注意. 理解好第一个样例就行 /* *********************************************** Author :bingone Created Time :2014/12/9 22:48:56 File Name :a.cpp ************************************************ */ #inclu

POJ 3254 状压DP

题目大意: 一个农民有一片n行m列 的农场   n和m 范围[1,12]  对于每一块土地 ,1代表可以种地,0代表不能种. 因为农夫要种草喂牛,牛吃草不能挨着,所以农夫种菜的每一块都不能有公共边. 告诉你 n ,m 和那些地方能种菜哪些地方不能种菜,求农夫一共有多少种方案种菜 解法: 基本思想是状压 也就是用一个int 型的数代表一行的种菜策略,二进制的0代表该位不能种菜,1位代表能种菜,使用位运算使处理速度变快 对于单行行,最多有2^12 种情况,并且 2^12种情况里面还有很多不满足题意的

POJ 2411 状压DP经典

Mondriaan's Dream Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 16771   Accepted: 9683 Description Squares and rectangles fascinated the famous Dutch painter Piet Mondriaan. One night, after producing the drawings in his 'toilet series

POJ 2686(状压DP

第一次做状压感觉那一长串for显示了这是个多么暴力的算法呢...1A了倒是挺顺的 #include<iostream> #include<cstdio> #include<algorithm> #include<queue> #include<utility> #include<vector> #include<cstring> #include<cmath> #define INF 0x3fffffff #d