poj1185:炮兵阵地(状压dp)

时间:2023-03-09 20:48:37
poj1185:炮兵阵地(状压dp)

也算是比较基础的状压dp了,跟做过的第二道比较又稍微复杂了一点

需要记录之前两行的状态。。

统计结果也稍有不同

另外还学习了一个得到一个整数二进制位 1 的个数的位运算方法

详见代码:

#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
#include<string>
#include<ctype.h>
using namespace std;
#define MAXN 10000
int s[];
int num[];
int cant[];
int dp[][][];
int ns;
int n,m;
int ct(int x)
{
int res=;
while(x)
{
res++;
x&=(x-);
}
return res;
}
void setstate()
{
memset(num,,sizeof(num));
for(int i=;i<<<;i++)
{
if((i&(i<<))==&&((i&(i<<))==))
{
s[ns]=i;
num[ns++]=ct(i);
}
}
}
void ini()
{
memset(cant,,sizeof(cant));
char s[];
for(int i=;i<=n;i++)
{
scanf("%s",s);
for(int j=;j<m;j++)
{
if(s[j]=='H')
{
cant[i]|=(<<j);
}
}
}
}
void solve()
{
memset(dp,,sizeof(dp));
for(int i=;i<ns&&s[i]<(<<m);i++)
{
if((s[i]&cant[])==)
dp[][][i]=num[i];
}
for(int i=;i<=n;i++)
{
for(int t=;t<ns&&s[t]<(<<m);t++)
{
if(s[t]&cant[i])
continue;
for(int j=;j<ns&&s[j]<(<<m);j++)
{
if(s[j]&cant[i-])
continue;
for(int k=;k<ns&&s[j]<(<<m);k++)
{
if(s[k]&cant[i-])
continue;
if(s[j]&s[k])
continue;
if(s[t]&cant[i])
continue;
if((s[t]&s[j])||(s[t]&s[k]))
continue;
dp[i][k][t]=max(dp[i][k][t],dp[i-][j][k]+num[t]);
}
}
}
}
int ans=;
for(int i=;i<=n;i++)
{
for(int j=;j<ns&&s[j]<(<<m);j++)
{
for(int k=;k<ns&&s[k]<(<<m);k++)
ans=max(ans,dp[i][j][k]);
}
}
printf("%d\n",ans);
}
int main()
{
ns=;
setstate();
while(scanf("%d%d",&n,&m)!=EOF)
{
ini();
solve();
}
return ;
}