JZYZOJ1390【noi2001】炮兵阵地 状压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 }
