POJ1185:炮兵阵地(状压dp)

题目:http://poj.org/problem?id=1185

大神的题解:

方法就是用DP[i][r][p]表示第i行状态为r,第i-1行状态是p时的最多个数。而这里p受到r的限制,而第i-2行状态q则受到r和p两个状态限制。状态转移方程就是:

DP[i][r][p] = MAX{DP[i-1][p][q] +num[r]}

其中,p是受到r的限制时枚举的状态,q是受到r和p共同限制时候的状态,num[r]表示状态r里面的布局炮兵所摆的个数。

这里我们可以看到就要枚举i,r,p,q,这4 个变量,i的范围是100,而其他几个则都是1<<10,复杂度颇为偏高。而实际上由于每一行里面有很多都是某些位置被其他位置影响的。比如:   1110001,  如果第一个位置放上炮兵,那么第二第三的位置都会受到影响,而一个也放不了。

解决方案就是不去管那些相互有影响的状态,把形如:

1000000    0100000    ... ...

1001000    0100001    ... ...

...

1001001

这些相互之间没有影响的状态找出来,这样所有的状态数就会减少至少于60种(我算了一下,好像是58种),这样一来就是60*60*60*100,可以过了。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#define mod 100000000
using namespace std;
int m,n,top,state[],cur[],num[],dp[][][];
char tu[][];
inline bool ok(int x)
{
if(x&(x<<)) return false;
if(x&(x<<)) return false;
return true;
}
inline void init()
{
int tol=<<n;
top=;
for(int i=; i<tol; i++)
{
if(ok(i)) state[++top]=i;
}
}
inline bool fit(int x,int k)
{
if(x&cur[k]) return false;
return true;
}
inline int icount(int x)
{
int cnt=;
while(x)
{
cnt++;
x=x&(x-);
}
return cnt;
}
int main()
{
while(scanf("%d%d",&m,&n)!=EOF)
{
init();
for(int i=; i<=m; i++)
scanf("%s",tu[i]+);
for(int i=;i<=m;i++)
{
cur[i]=;
for(int j=;j<=n;j++)
{
if(tu[i][j]=='H')
cur[i]+=(<<(n-j));
}
}
memset(dp,,sizeof(dp));
for(int i=; i<=top; i++)
{
num[i]=icount(state[i]);
if(fit(state[i],))
dp[][][i]=num[i];
}
for(int i=; i<=m; i++)
{
for(int t=; t<=top; t++)
{
if(!fit(state[t],i)) continue;//符合本行
for(int j=; j<=top; j++)
{
if(state[j]&state[t]) continue;//符合上一行
for(int k=; k<=top; k++) //符合上上行
{
if(state[k]&state[t]) continue;
if(state[k]&state[j]) continue;
dp[i][j][t]=max(dp[i][j][t],(dp[i-][k][j]+num[t]));
}
}
}
}
int maxx=-;
for(int i=; i<=m; i++)
{
for(int j=; j<=top; j++)
{
for(int k=; k<=top; k++)
{
maxx=max(maxx,dp[i][j][k]);
}
}
}
printf("%d\n",maxx);
}
return ;
}
上一篇:剑指架构师系列-spring boot的logback日志记录


下一篇:Vue Route Building the UI back-end framework