JZYZOJ1390noi2001炮兵阵地 状压DP

Posted 鲸头鹳

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了JZYZOJ1390noi2001炮兵阵地 状压DP相关的知识,希望对你有一定的参考价值。

http://172.20.6.3/Problem_Show.asp?id=1390

需要储存该行和上一行两个状态。通过观察规则可以发现条件允许的状态很少(相邻两个至少空两格),据此可以减少状态数量,从而极大压缩空间和时间。
代码
技术分享
 1 #include<iostream>  
 2 #include<cstdio>  
 3 #include<cstring>  
 4 #include<algorithm>  
 5 #include<cmath>
 6 using namespace std;
 7 const int maxn=110;
 8 int n,m,ma;
 9 char ch[12]={};
10 int ke[60]={},d[60]={},tot=0;
11 int a[maxn]={},f[maxn][60][60]={};
12 bool ca(int z){
13     if(z&(z<<1))return 0;
14     if(z&(z<<2))return 0;
15     return 1;
16 }
17 int main(){
18     memset(f,-1,sizeof(f));
19     scanf("%d%d",&n,&m);ma=1<<m;
20     for(int i=1;i<=n;i++){
21         scanf("%s",&ch);
22         for(int j=0;j<m;j++)
23             if(ch[j]==H) a[i]|=(1<<(m-j-1));
24     }
25     for(int i=1;i<ma;i++){
26         if(ca(i)){
27             ke[++tot]=i;
28             for(int j=0;j<m;j++)
29                 if(i&(1<<j))d[tot]++;
30         }
31     }
32     for(int i=0;i<=tot;i++){
33         int w=a[1]&ke[i];
34         if(!w)
35             f[1][0][i]=d[i];
36     }int ans=0;
37     for(int i=2;i<=n;i++){
38         for(int w=0;w<=tot;w++){
39             for(int j=0;j<=tot;j++){
40                 if(ke[w]&ke[j])continue;
41                 if(ke[j]&a[i])continue;
42                 for(int j1=0;j1<=tot;j1++){
43                     if((ke[j1]&ke[j])||(ke[j1]&ke[w]))continue;
44                     f[i][w][j]=max(f[i-1][j1][w],f[i][w][j]);
45                 }
46                 f[i][w][j]+=d[j];
47                 ans=max(f[i][w][j],ans);
48             }
49         }
50     }printf("%d\n",ans);
51     return 0;
52 }
View Code

 

以上是关于JZYZOJ1390noi2001炮兵阵地 状压DP的主要内容,如果未能解决你的问题,请参考以下文章

P2704 [NOI2001] 炮兵阵地(状压dp)

状压DP NOI2001 炮兵阵地

(状压dp)NOI 2001(POJ 1185) 炮兵阵地

NOI2001炮兵阵地

[NOI2001]炮兵阵地

NOI2001炮兵阵地