P2704 [NOI2001] 炮兵阵地（状压dp）

题目传送门

题意： 给你一张N*M大小的由P和H组成的地图，在点P上你可以放置炮兵，炮兵的攻击范围如图，问你这张地图中最多放多少个炮兵从而让他们不会互相攻击。（N<=100,M<=10）

思路： 每一行的状态用二进制压缩，1为放置炮兵，0为不放置。最开始把每一行中炮兵不会互相攻击的状态全部处理出来，然后输入地图的时候把所有H看成1，方便后面枚举的时候判断某个状态是否在第i行适用（两个状态按位与为0则适用）。用f[i][j][k]表示前i行，第i行状态为j，并且第i-1行状态为k时，可以放置最大炮兵数量，状态转移方程也就很简单了。

代码：

#include<bits/stdc++.h>
#define endl '\n'
#define mp make_pair
#define pb push_back
#define ll long long
//#define int long long
#define pii pair<int,int>
#define sz(x) (int)(x).size()
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read()
{
    int x=0,f=1;
    char ch=gc();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            f=-1;
        ch=gc();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=gc();
    }
    return x*f;
}
using namespace std;
const int N=1e2+10;
const int inf=0x3f3f3f3f;
const int mod=998244353;
const double eps=1e-7;
string s[N];
int sr[N],sum[N],cnt=0;
int f[111][111][111],n,m,line[N];
void dfs(int cur,int pos,int x)
{
    if(pos >= m)
    {
        sum[++cnt] = x;
        sr[cnt] = cur;
        return ;
    }
    dfs(cur,pos+1,x);
    dfs(cur|(1<<pos),pos+3,x+1);
}
int solve()
{
    cin>>n>>m;
    dfs(0,0,0);
    for(int i=1; i<=n; i++)
    {
        cin>>s[i];
        for(int j=0;j<s[i].size();j++)
        {
            if(s[i][j]=='H')
                line[i] |= (1<<j);
        }
    }
    for(int i=1;i<=cnt;i++)
    {
        if(sr[i]&line[1])
            continue;
        for(int j=1;j<=cnt;j++)
        {
            if(sr[i]&sr[j]) continue;
            f[1][i][j] = sum[i];
        }
    }
    for(int i=2;i<=n;i++)
    {
        for(int j=1;j<=cnt;j++)
        {
            if(sr[j]&line[i]) continue;
            for(int k=1;k<=cnt;k++)
            {
                if(sr[j]&sr[k]) continue;
                if(sr[k]&line[i-1]) continue;
                for(int l=1;l<=cnt;l++)
                {
                    if(sr[j]&sr[l]) continue;
                    if(sr[k]&sr[l]) continue;
                    if(sr[l]&line[i-2]) continue;
                    f[i][j][k] = max(f[i][j][k],f[i-1][k][l] + sum[j]);
                }
            }
        }
    }
    int res=0;
    for(int i=1;i<=cnt;i++)
        for(int j=1;j<=cnt;j++)
            res = max(res,f[n][i][j]);
    cout<<res<<endl;
    return 0;
}
signed main()
{
//    int _;
//    cin>>_;
//    while(_--)
    solve();
}