炮兵阵地 解题报告
状态压缩例题；二进制、位运算的理解与应用
版权声明：欢迎大家指正错误，有不同观点的欢迎评论，共同进步 https://blog.csdn.net/Sirius_han/article/details/79959342
G - 炮兵阵地
题目链接：POJ - 1185POJ-3254的升级版；
题意：N*M的地图标有H,G， H不能建炮兵阵地， P可以，炮兵阵地之间至少要间距两个格子（POJ-3254是不能相邻）；若用1表示有炮兵阵地， 0表示无， 则1001001可行， 1100和1010不行；
思路：POJ-3254的每行只被上一行影响， 这道题的每行被前两行影响；
#include <iostream>
#include <algorithm>
#include <string.h>
using namespace std;
int n, m, top, dp[105][100][100], fild[105], state[100], cnt[100];//打表可知， 符合条件的状态只有60多个；
bool ok(int x){
if(x&(x<<1)||x&(x<<2)) return false;
return true;
}
int Count(int x){
int ans=0;
while(x){
if(x&1) ans++;
x>>=1;
}
return ans;
}
bool fit(int x, int i){
if(x&fild[i]) return false;
return true;
}
void init(){
top=0;
for(int i=0; i<(1<<m); i++){
if(ok(i)){
state[top]=i;
cnt[top++]=Count(i);
}
}
}
int main(){
char ch;
cin >>n >> m;
init();
for(int i=1; i<=n; i++){
fild[i]=0;
for(int j=1; j<=m; j++){
cin >>ch;
if(ch=='H') fild[i]+=(1<<(m-j));
}
}
memset(dp, 0, sizeof(dp));
for(int i=0; i<top; i++){
if(fit(state[i], 1)) dp[1][i][0]=cnt[i];
}
for(int i=0; i<top; i++){
if(!fit(state[i], 2)) continue;
for(int j=0; j<top; j++){
if(!fit(state[j], 1)) continue;
if(state[i]&state[j]) continue;
dp[2][i][j]=max(dp[2][i][j], dp[1][j][0]+cnt[i]);
}
}
for(int i=3; i<=n; i++){
for(int j=0; j<top; j++){
if(!fit(state[j], i)) continue;
for(int k=0; k<top; k++){
if(!fit(state[k], i-1)) continue;
if(state[j]&state[k]) continue;
for(int l=0; l<top; l++){
if(!fit(state[l], i-2)) continue;
if(state[j]&state[l]) continue;
if(state[k]&state[l]) continue;
dp[i][j][k]=max(dp[i][j][k], dp[i-1][k][l]+cnt[j]);
}
}
}
}
int ans=0;
for(int i=0; i<top; i++){
for(int j=0; j<top; j++){
ans=max(ans, dp[n][i][j]);
}
}
cout <<ans <<endl;
return 0;
} 
