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外,以上所有数据皆为正整数)
设从下往上数第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 = πR2H
侧面积A' = 2πRH
底面积A = πR2
体积V = πR2H
侧面积A' = 2πRH
底面积A = πR2
非常多的神奇的剪枝,先求出前缀和数组便于剪枝。
几处剪枝 1、当当前体积+剩下几层的最小体积>总体积n 不符合条件,剪枝。
2、当前表面积+剩下几层最小表面积>ans 没必要进行下去,剪枝。
3、(ans-当前表面积)<=(总体积-当前体积)/r*2 大于最小可能的剩余表面积,没必要进行,一定不是最优解 体积V = πR2H r和h对表面积的影响都是呈线性关系,r呈指数倍关系对V产生影响。所以说,(总体积-当前体积)/r*2为最小可能剩余面积,不是最大可能剩余面积
4、下一次搜索的h可能的最大值,在h-1和(总体积-当前体积-剩下几层的最小体积)/(i*i) 选择比较小的。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
int n,m,ans;
int Vmin[25],Smin[25];
void init()//自下而上,一层为最上面的一层
{
for(int i=1;i<=20;i++)
{
Vmin[i]=Vmin[i-1]+i*i*i;//Vmin[1]表示最上层的最小体积,为1;
Smin[i]=Smin[i-1]+i*i*2;
}
}
void dfs(int v,int s,int cnt,int r,int h)
{
if(v>n) return;
if(!cnt)
{
if(v==n&&s<ans)
ans=s;
return;
}
if(v+Vmin[cnt-1]>n||s+Smin[cnt-1]>ans||(n-v)/r*2+s>=ans) return;
for(int i=r-1;i>=cnt;--i)
{
if(cnt==m) s=i*i;
int Max=min((n-v-Vmin[cnt-1])/(i*i),h-1);//下一次要搜索的高的最大值
for(int j=Max;j>=cnt;--j)
{
dfs(v+i*i*j,s+2*i*j,cnt-1,i,j);
}
}
}
int main()
{
scanf("%d%d",&n,&m);
init();
ans=0x3f3f3f3f;
dfs(0,0,m,n,n);//从最下层开始搜索
if(ans==0x3f3f3f3f)
ans=0;
printf("%d",ans);
return 0;
}