poj1185 炮兵阵地 状压dp

题目链接:

https://vjudge.net/contest/159644#problem/B

题意:

在一个N*M的矩阵上布置炮兵部队，只有平原可以布置，然后每个炮兵部队都有一个攻击范围，它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。
问最多可以安排几个炮兵

题解:

http://www.cnblogs.com/buptLizer/archive/2012/08/23/2651881.html
状压dp
dp[i][j][k]:=考虑到第i行， 这一行是第i个合法状态，上一行是第j个合法状态，最多可以安排的炮兵
转移：dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][p]+num[j]); num[j]表示j这个合法状态可以安排的炮兵数

代码:

 1 #include <cstdio>
 2 #include <iostream>
 3 #include <cstring>
 4 using namespace std;
 5 typedef long long ll;
 6 #define MS(a) memset(a,0,sizeof(a))
 7 #define MP make_pair
 8 #define PB push_back
 9 const int INF = 0x3f3f3f3f;
10 const ll INFLL = 0x3f3f3f3f3f3f3f3fLL;
11 inline ll read(){
12     ll x=0,f=1;char ch=getchar();
13     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
14     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
15     return x*f;
16 }
17 //////////////////////////////////////////////////////////////////////////
18 const int maxn = 80+10;
19 
20 int dp[105][maxn][maxn];
21 int state[maxn],num[maxn],cnt;
22 int cur[105];
23 char mp[105][20];
24 int n,m;
25 
26 int getnum(int x){
27     int res = 0;
28     while(x){
29         res += (x&1);
30         x = x>>1;
31     }
32     return res;
33 }
34 
35 void init(){
36     cnt = 0;
37     for(int i=0; i<(1<<m); i++){
38         if(!(i&(i<<1)) && !(i&(i<<2))){
39             state[cnt] = i;
40             num[cnt] = getnum(i);
41             cnt++;
42         }
43     }
44 }
45 
46 int main(){
47     cin >> n >> m;
48     for(int i=0; i<n; i++){
49         cin >> mp[i];
50         for(int j=0; j<m; j++){
51             if(mp[i][j] == 'P')
52                 cur[i] += (1<<(m-j-1));
53         }
54     }
55 
56     memset(dp,-1,sizeof(dp));
57     init();
58     for(int i=0; i<cnt; i++)
59         if(!(state[i]&(~cur[0]))) // 如果与能放炮兵的反有交集 就不可行
60             dp[0][i][0] = num[i];
61     for(int i=1; i<n; i++){
62         for(int j=0; j<cnt; j++){
63             if((~cur[i])&state[j]) continue;
64             for(int k=0; k<cnt; k++){
65                 if((~cur[i-1])&state[k]) continue;
66                 if(state[j]&state[k]) continue;
67                 for(int p=0; p<cnt; p++){
68                     if(i>=2 && ((~cur[i-2])&state[p])) continue;
69                     if(state[j]&state[p]) continue;
70                     if(dp[i-1][k][p] == -1) continue;
71                     dp[i][j][k] = max(dp[i][j][k],dp[i-1][k][p]+num[j]);
72                 }
73             }
74         }
75     }
76 
77     int ans = 0;
78     for(int i=0; i<cnt; i++)
79         for(int j=0; j<cnt; j++)
80             ans = max(ans,dp[n-1][i][j]);
81     cout << ans << endl;
82 
83     return 0;
84 }
85 
86 // http://www.cnblogs.com/buptLizer/archive/2012/08/23/2651881.html

 

posted @ 2017-04-21 20:43  _yxg123  阅读(203)  评论(0)  编辑  收藏  举报