链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=2167
题目大意:
给你一个N*N( 3<=N<=15)个矩阵,要你选择若干个数(每个数的范围在[10,99]),使得最后所选的数总和最大。选数的规则是如果选了某个数,那么它的八个相邻方向的数都不能选。
这道题其实可以说是状压的模版题了,但是第一眼看到这题还是一脸懵逼,因为之前没学过状压嘛。什么也不懂。后来scy给我们讲了之后,才艰难的码了出来。
状态压缩到底是个什么玩意儿呢,让我们来看看这个算法到底适用于那种题目吧:
从状态压缩的特点来看,这个算法适用的题目符合以下的条件:
1.解法需要保存一定的状态数据(表示一种状态的一个数据值),每个状态数据通常情况下是可以通过2进制来表示的。这就要求状态数据的每个单元只有两种状态,比如说棋盘上的格子,放棋子或者不放,或者是硬币的正反两面。这样用0或者1来表示状态数据的每个单元,而整个状态数据就是一个一串0和1组成的二进制数。
2.解法需要将状态数据实现为一个基本数据类型,比如int,long等等,即所谓的状态压缩。状态压缩的目的一方面是缩小了数据存储的空间,另一方面是在状态对比和状态整体处理时能够提高效率。这样就要求状态数据中的单元个数不能太大,比如用int来表示一个状态的时候,状态的单元个数不能超过32(32位的机器)。 引用自 Avril的博客
而这道题很明显都符合上面的两个条件,所以,我们采取状压来做这道题。我们把每一行的状态都变成二进制,如果你要取这一个格子里的数,那么就用1来表示,若不取,就用0来表示。
附代码:
#include<cstdio>
#include<sstream>
#include<cstring>
using namespace std;
int f[16][(1<<15)+10];
int a[110][110],v[1<<15],vn,bit[110];
//v记录合法的方案,vn记录合法的方案数
//a记录输入的矩阵
//f[i][p]表示当前第i行采取的合法的方案p可在这一行获取的最大和
char s[1024];
int n,maxx;
int main()
{
bit[1]=1;for(int i=2;i<=15;i++)bit[i]=(bit[i-1]<<1);//先存下2的1~15次方,方便以后运用
while(scanf("%d",&a[1][1])!=EOF)
{
char ss;
scanf("%c",&ss);
n=1;
while(ss!='\n')
{
n++;
scanf("%d",&a[1][n]);
scanf("%c",&ss);
}
for(int i=2;i<=n;i++)for(int j=1;j<=n;j++)scanf("%d",&a[i][j]);
getchar();//清除空格与换行,scy说这句不要也罢,你们看着办
memset(f,0,sizeof(f));
vn=0;
maxx=(1<<n)-1;//maxx表示最多有maxx种方案(包括不合法方案)
for(int x=0;x<=maxx;x++)
{
if(((x<<1)& x) ==0)
{
for(int j=1;j<=n;j++)
{
if((bit[j]& x) ==bit[j])
f[1][x]+=a[1][j];
}
v[++vn]=x;
}
}
for(int i=2;i<=n;i++)
{
for(int p=1;p<=vn;p++)
{
int tt=0;for(int j=1;j<=n;j++)if(bit[j] & v[p])tt+=a[i][j];
for(int q=1;q<=vn;q++)
{
if( (v[p] & v[q])==0 && ((v[q]<<1) & v[p])==0 && ((v[q]>>1) & v[p])==0)
{
f[i][v[p]]=max(f[i][v[p]],tt+f[i-1][v[q]]);
}
}
}
}
int ans=0;for(int p=1;p<=vn;p++)if(ans<f[n][v[p]])ans=f[n][v[p]];
printf("%d\n",ans);
}
return 0;
}