ACcoders Problem 2051 题解

题意

如上图，有一个生日蛋糕，一共有 $m$ 层（层数从上到下递增），体积为 $n$（$\pi$ 忽略不计），每层的半径和高都严格单调递增，问符合这些条件的生日蛋糕表面积最小为多少（不算底面积）？

思路

如上图，我们可以观察到红色部分拼合起来正好为一个底面积（此后式子中 $\pi$ 忽略不计），所以面积即为$r_m^2$

然后我们再计算蓝色部分。侧面积为 $\sum _{i=1}^{m}2r_i{h}_i$，且根据题意中的单调性，又可以列出两条式子，分别为
$$r_i\le r_{i+1}-1,h_i\le h_{i+1}-1$$

考虑搜索，每次传入三个参数，分别为当前层数 $u$，第 $m$ 到 $u+1$ 层的总表面积 $S$，第 $m$ 到 $u+1$ 层的总体积 $V$。

我们考虑再预处理出两个数组，一个是从 $1$ 到 $u$ 层的最小侧面积和 $minS$，另一个是从 $1$ 到 $u$ 层的最小体积和 $minV$。根据刚刚的式子，每次半径最小与相邻层差1，所以从 $1$ 到 $m$ 层依次枚举，每次的半径也就是对应的层数，随后以 $i$ 为半径来套公式。

随后，我们进行剪枝操作。

剪枝 1，搜索顺序

考虑从底向顶枚举层数。枚举 $h,r$ 时应该先枚举 $r$，因为在圆柱体积公式中，$r$ 是二次幂，而 $h$ 是一次幂，显然 $r$ 贡献更大。还是根据刚刚的式子所以第 $u$ 层半径最小也就为 $u$，$h$ 也是同理。那最大值该怎么来算呢？又是根据一开始列出的式子，与相邻层数差 $1$，所以最多能取到 $r_{i+1}-1$，$h$ 也是同理。但是还没完，因为总体积最多为 $n$，所以现在体积最多还剩下 $n-V$，那么当前层也就是 $r^2\times h\le n-v$（有可能没全用完，所以是小于等于），设此时 $h$ 为它最小的时候 $1$，那么这个式子就能变成
$$r\le \sqrt{n-v},h\le \frac{n-v}{r^2}$$

这也是一种最大值可能，因为要使搜索过程更优，所以取最小值。

剪枝 2，可行性剪枝

如果当前已经有的体积加上还没被用过的体积和的最小值（就是我们之前已处理的 $minV$），已经不符合条件了，也就是大于 $n$ 了，就不可能有符合要求的答案了，直接返回

剪枝 3，最优性剪枝

如果当前已经有的表面积加上还没被用过的表面积和的最小值（就是我们之前已处理的 $minS$），不比我们的最优结果更优，再求它也就没有意义了，直接返回。

剪枝 4，聪明の剪枝

我们还可以考虑一下表面积与体积之间的关系。根据刚才的最优性剪枝，如果当前已经有的表面积加上还没被用过的表面积和的最小值不比我们的最优结果更优就返回，如果我们推出一个关系，还没被用过的表面积和的最小值大于 $x$，那么这个 $x$ 加上当前已经有的表面积还大于最优结果，那也可以再剪枝。结论为 $x=\frac{2(n-v)}{r_{u+1}}$，以下为证明（从第 $1$ 到 $u$ 层的表面积记为 $S_{1-u}$）：
$$S_{1-u}=\sum _{i=1}^{u}2r_i{h}_i=\sum _{i=1}^{u}2r_i{h}_i\times r_{u+1}\frac{1}{r_{u+1}}=\frac{2}{r_{u+1}}\sum _{i=1}^u{r}_i{h}_i{r}_{u+1}>\frac{2}{r_{u+1}}\sum _{i=1}^u{r}_i^2{h}_i$$

我们知道，$\sum _{i=1}^u{r}_i^2{h}_i$ 即为 $n-V$，又得
$$S_{1-u}>\frac{2(n-V)}{r_{u+1}}$$

所以，当 $S+\frac{2(n-V)}{r_{u+1}}>ans$ 的时候，也可以剪枝。

代码

#include<bits/stdc++.h>
using namespace std;
const int INF=0x3f3f3f3f;
int mins[1000001];
int minv[1000001];
int R[100001],H[100001];
int n,m,ans=INF;
void dfs(int u,int S,int V)
{
	if(V+minv[u]>n) return ;
	if(S+mins[u]>=ans) return ;
	if(S+2*(n-V)/R[u+1]>=ans) return ;
	if(!u)
	{
		if(V==n) 
			ans=min(ans,S);
		return ;
	}
	for(int i=min(R[u+1]-1,(int)sqrt(n-V));i>=u;i--)
		for(int j=min(H[u+1]-1,(n-V)/i/i);j>=u;j--)
		{
			int t=0;
			if(u==m) t=i*i;
			R[u]=i;H[u]=j;
			dfs(u-1,S+2*i*j+t,V+i*i*j);
		}
	return ;
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		minv[i]=minv[i-1]+i*i*i;
		mins[i]=mins[i-1]+2*i*i;
	}
	R[m+1]=H[m+1]=INF;
	dfs(m,0,0);
	if(ans==INF) puts("-1");
	else printf("%d",ans);
	return 0;
}