《GMOJ-Senior-2679 跨时代》题解

题目大意

给出$n$根栏杆的长度，让你从这$n$根栏杆中选出一些栏杆围成一个矩形（必须要刚好围成一个矩形，即不能出现多余的边长，且不能切断栏杆，但所给栏杆不一定要全部用上），求围成矩形的最大面积。
对于$30\%$的数据，$1\le n\le 10$；
对于$100\%$的数据，$1\le n\le 16$,$1\le \text{栏杆的长度}\le 15$。

分析

这是一道状态压缩动态规划（状压$DP$）题。

30分

对于$30\%$的数据，我们可以用暴力求出。设围成的矩形的四条边分别为$e_1,e_2,e_3,e_4$，那么在$DFS$时枚举每条边的每种情况（加入$e_1$或$e_2$或$e_3$或$e_4$），最后判断矩形的长宽是否相等，如果相等就更新答案。代码如下：

#include<cstdio>
int max(int x,int y)//最大值函数
{
	if(x>y)
	{
		return x;
	}
	return y;
}
int n,a[22];
int dfs(int s/*枚举的栏杆的编号*/,int e1/*矩形的第一条长*/,int e2/*矩形的第二条长*/,int e3/*矩形的第一条宽*/,int e4/*矩形的第二条宽*/)
{
	if(s==n+1)//栏杆枚举完了
	{
		if(e1==e2&&e3==e4)//判断矩形是否合法（长、宽是否相等）
		{
			return e1*e3;
		}
		return 0;
	}
	return max(max(max(dfs(s+1,e1+a[s],e2,e3,e4),//加入e1
					   dfs(s+1,e1,e2+a[s],e3,e4)),//加入e2
				   max(dfs(s+1,e1,e2,e3+a[s],e4),//加入e3
					   dfs(s+1,e1,e2,e3,e4+a[s]))),//加入e4
					   dfs(s+1,e1,e2,e3,e4));//不选
}
int main()
{
	scanf("%d",&n);//读入n
	for(int i=1;i<=n;i++)//读入栏杆长度
	{
		scanf("%d",&a[i]);
	}
	int t=dfs(1,0,0,0,0);//DFS
	if(t==0)//如果无解
	{
		printf("No Solution");
	}
	else
	{
		printf("%d",t);//输出答案
	}
	return 0;
}

这样的时间复杂度为$O(5^n)$，$100\%$的数据肯定过不了，于是我们考虑对其优化。

100分

因为我们知道矩形的长、宽相等，所以我们设它的长为$a$，宽为$b$，深搜时每根栏杆的情况就由$5$种降到了$3$种：加入$a$，加入$b$，不选。最后我们再判断$a$，$b$是否可以分成相等的两部分，如果可以更新答案即可。
如何判断$a$，$b$是否可以分成相等的两部分呢？我们用$f$数组来判断。$f_i$（$i$为一个$01$串）表示对于$i$这种选择栏杆的情况中，选出的栏杆能否分成相等的两段。因为$i$是一个$01$串，所以我们可以把它压缩成一个十进制数储存。
对于每个$i$，我们先求出选出的栏杆的长度总和$sum$。如果$sum$是奇数，那么显然$f_i=false$。否则我们用$01$不带权背包$g$辅助处理。我们把选出的栏杆逐个放入$g$中，然后判断$g_{sum÷2}$（一半）能否被组合，若能，表示$sum-sum÷2$（即另一半）也能被组合，$f_i=true$；否则，$f_i=false$。
根据以上思路，我们可以打出如下代码：

#include<cstdio>
bool f[1<<16],g[242];int n,l[22];
int max(int x,int y)//最大值函数
{
	if(x>y)
	{
		return x;
	}
	return y;
}
int dfs(int s/*枚举的栏杆的编号*/,int a/*长的选择情况*/,int suma/*两条长的总和*/,int b/*宽的选择情况*/,int sumb/*两条宽的总和*/)
{
	if(s==n+1)//栏杆枚举完了
	{
		if(f[a]&&f[b])//判断矩形是否合法（长、宽能否被分成相等的两部分）
		{
			return (suma/2)*(sumb/2);//注意这是两条长、宽的总和，乘的时候记得分别除以2
		}
		return 0;
	}
	return max(max(dfs(s+1,a+(1<<(s-1)),suma+l[s],b,sumb)/*加入长*/,dfs(s+1,a,suma,b+(1<<(s-1)),sumb+l[s]))/*加入宽*/,dfs(s+1,a,suma,b,sumb)/*不选*/);
}
int main()
{
	scanf("%d",&n);//读入n
	for(int i=1;i<=n;i++)//读入读入栏杆长度
	{
		scanf("%d",&l[i]);
	}
	int m=1<<n//求出栏杆选择的情况总数
	for(int i=0;i<m;i++)//枚举每种情况
	{
		int sum=0;//求出选出的栏杆长度总和
		for(int j=1;j<=n;j++)//枚举每一条边
		{
			if((i&(1<<(j-1)))!=0)//判断是否被选中
			{
				sum+=l[j];//累加长度
			}
		}
		if(sum%2!=0)//如果长度除以2不是整数（无法平均分成相等的两个整数部分）
		{
			f[i]=0;
		}
		else
		{
			g[0]=1;//初始化01不带权背包
			for(int j=1;j<=sum;j++)
			{
				g[j]=0;
			}
			for(int j=1;j<=n;j++)//枚举每一条边
			{
				if((i&(1<<(j-1)))!=0)//判断是否被选中
				{
					for(int k=sum;k>=l[j];k--)//放入背包
					{
						g[k]|=g[k-l[j]];
					}
				}
			}
			f[i]=g[sum/2];//判断能否被分成相等的两部分
		}
	}
	int t=dfs(1,0,0,0,0);//DFS
	if(t==0)//如果无解
	{
		printf("No Solution");
	}
	else
	{
		printf("%d",t);//输出答案
	}
	return 0;
}

总结

在做像这样的状压$DP$题时，可以先打出暴力，然后再尝试对其优化。