Description
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,找出蛋糕的制作方案(适当的Ri和Hi的值),使S最小。
(除Q外,以上所有数据皆为正整数)
Input
有两行,第一行为N(N <= 10000),表示待制作的蛋糕的体积为Nπ;第二行为M(M <= 20),表示蛋糕的层数为M。
Output
仅一行,是一个正整数S(若无解则S = 0)。
Sample Input
100
2
Sample Output
68
Hint
圆柱公式
体积V = πR 2H
侧面积A’ = 2πRH
底面积A = πR 2
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
int inf=0x3f3f3f3f;
int m,n,s;
int Maxh,Maxr,Minv[30],Mins[30];
void dfs(int sums,int sumv,int r,int h,int depth)//从下往上搜索
{
if(depth==0)
{
if(sumv==n&&sums<s)
s=sums;
return;
}
if(sums+Mins[depth-1]>s||sumv+Minv[depth-1]>n||2*(n-sumv)/r+sums>=s)//如果当前所算之和与未算部分的可能最小值之和大于目标,则没有必要再计算下去
return;
for(int i=r-1; i>=depth; i--)
{
if(depth==m)
sums=i*i;//初始值
int hh=min((n-sumv-Minv[depth-1])/i/i,h-1);//当前层的最大高度
for(int j=hh; j>=depth; j--)
dfs(sums+2*i*j,sumv+i*i*j,i,j,depth-1);
}
}
int main()
{
s=inf;
scanf("%d%d",&n,&m);
Maxh=n,Maxr=sqrt(n);//最大高度及最大半径
Minv[0]=0,Mins[0]=0;
for(int i=1; i<=m; i++)
{
Minv[i]=Minv[i-1]+i*i*i;//从上往下依次编号,此即到此层的最小总体积
Mins[i]=Mins[i-1]+i*i;
}
dfs(0,0,Maxr,Maxh,m);//从下往上搜索
if(s!=inf)
printf("%d\n",s);
else
printf("0\n");
return 0;
}
Reflect:
1,通过从上到下预置体积和面积之和,为剪枝提供依据;
2,预置是从上到下的,而搜索是从下到上的(很关键);
3,对于每一层的半径和高的搜索范围!
4,最大半径和最大高度只需将H和R分别取一即可;