[POJ]poj1185 炮兵营地(状压DP)

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示：

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。

N <= 100；M <= 10。

ACM课作业，基础状压DP，实际合法状态数很少，O(n*state^3)可以接受。

#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
#define rep(i,n) for(int i=1;i<=n;i++)
#define st first
#define nd second 
#define pii pair<int,int>
#define ll long long
#define pll pair<ll,ll>
#define pb push_back
#define mp make_pair


const int N=1e2+7;
const int M=1e2+7;

int n,m;

int st[N],f[N][M][M],sum[M];

vector<int> state;

int main(){
    ios::sync_with_stdio(false);
    cin>>n>>m;
    rep(i,n){
        string s;
        cin>>s;
        rep(j,m)st[i]=st[i]*2+ (s[j-1]=='H'? 1: 0);
    }
    
    int maxstate=1<<m;
    for(int i=0;i<maxstate;i++){
        if( !(i&(i<<1)) && !(i&(i<<2)) ){
            state.pb(i); 
            int res=i;
            while(res){
                if(res&1)sum[state.size()-1]++;
                res>>=1;
            }
        }
    }
    
    memset(f,-1,sizeof(f));
    f[0][0][0]=0;
    rep(i,n){
        for(int j=0;j<state.size();j++){ 
            if(state[j]&st[i])continue; 
            for(int k=0;k<state.size();k++){
                if(state[k]&state[j])continue;
                for(int l=0;l<state.size();l++)
                    {
                        if(state[l]&state[j])continue;
                        if(f[i-1][k][l]!=-1){
                            f[i][l][j]=max(f[i-1][k][l]+sum[j],f[i][l][j]);
                        //    printf("%d %d %d
",f[i-1][k][l],sum[j],f[i][l][j]); 
                        }
                            }            
            }
            }    
    }
    int ans=0;
    for(int i=0;i<state.size();i++)
        for(int j=0;j<state.size();j++)ans=max(ans,f[n][i][j]);
        
    cout<<ans;
    
}