GDOI Training T7 生日蛋糕

题目链接：传送门

题面

题目描述
	7月17日是Mr.W的生日，ACM-THU为此要制作一个体积为Nπ的M层生日蛋糕，每层都是一个圆
柱体。设从下往上数第i(1 <= i <= M)层蛋糕是半径为Ri, 高度为Hi的圆柱。当i < M时，要求Ri>
Ri+1且Hi > Hi+1。由于要在蛋糕上抹奶油，为尽可能节约经费，我们希望蛋糕外表面（最下一层
的下底面除外）的面积Q最小。 令Q = Sπ 请编程对给出的N和M，找出蛋糕的制作方案（适当的R
i和Hi的值），使S最小。 （除Q外，以上所有数据皆为正整数）

输入格式
	有两行，第一行为N（N <= 10000），表示待制作的蛋糕的体积为Nπ；第二行为M(M <= 20)，
表示蛋糕的层数为M。

输出格式
	仅一行，是一个正整数S（若无解则S = 0）。

Step 1.爆搜&几个浅显易懂的剪枝

设每层半径分别为$r_1$，$r_2$，…，$r_n$，高度分别为$h_1$，$h_2$，…，$h_n$。
首先可以得出
$N=\sum r_i^2{h}_i$
$S=2r_i^2+2\sum r_i{h}_i$

那么可以先枚举第一层半径和高度，然后依次枚举上面一层的半径和高度。
浅显易懂的剪枝1——设目前表面积为S，令ans保存最终解，若s>ans则可以直接返回。
浅显易懂的剪枝2——设目前体积为V，若V>N则可以直接返回。

Step 2.对于体积V的剪枝

简单地说：
如果剩下的层数均采用最小的r和h时最终的V仍会大于N，则返回。
如果剩下的层数均采用最大的r和h时最终的V仍会小于N，则返回。

然而，均采用最大的r和h时会遇到一个尴尬的问题：
由于$r_i$>$r_{i+1}$，$h_i$>$h_{i+1}$，所以剩下的层数均采用最大的r和h时，增加的体积无法在O(1)的时间内求出！

所以这里要对增加的体积进行放缩，贴个代码：

	if(V+(lr-1)*(lr-1)*(lh-1)*(m-deep)<n)	return;

其中lr表示上一层的半径，lh表示上一层的高度，m表示总层数，deep表示当前搜索深度。

Step 3.对于表面积S的剪枝

与体积的剪枝类似，当剩下的层数都取最小的表面积时，最终的S仍然大于ans时，返回。
这里同样要对增加的表面积进行放缩，贴个代码：

if(S+(n-V)*2/i>ans)	continue;

其中S代表当前表面积，n代表最终要达到的体积，V代表当前体积，i为这一层的半径。

Step 4.Coding!

注：略微卡了下常。

#include<stdio.h>
#include<cstring>
#include<algorithm>
#define re register int
using namespace std;
int read()
{
	re x=0,f=1;
	char ch=getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch=='-')	f=-1;
		ch=getchar();
	}
	while(ch>='0' && ch<='9')
	{
		x=10*x+ch-'0';
		ch=getchar();
	}
	return x*f;
}
const int Size=1001;
const int INF=0x3f3f3f3f;
int n,m,ans;
int sqr(int x) { return x*x; }
void dfs(int deep,int V,int S,int lr,int lh)
{
	if(S>ans || V>n)	return;
	if(V+sqr((m-deep)*(m-deep+1)>>1)>n)	return;
	if(V+(lr-1)*(lr-1)*(lh-1)*(m-deep)<n)	return;
	if(deep==m)
	{
		if(V==n && S<ans)
		{
			ans=S;
		}
		return;
	}
	for(re i=m-deep; i<lr; i++)
	{
		if(S+(((n-V)<<1)/i)>ans)	continue;
		for(re j=m-deep; j<lh; j++)
		{
			dfs(deep+1,V+i*i*j,S+(i*j<<1),i,j);
		}
	}
}
int main()
{
	n=read();
	m=read();
	ans=INF;
	for(re i=m; i<=100; i++)
	{
		for(re j=m; j<=200; j++)
		{
			dfs(1,i*i*j,i*i+2*i*j,i,j);
		}
	}
	printf("%d",ans);
	return 0;
}