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外,以上所有数据皆为正整数)
有两行,第一行为N(N <= 10000),表示待制作的蛋糕的体积为Nπ;第二行为M(M <= 20),表示蛋糕的层数为M。
仅一行,是一个正整数S(若无解则S = 0)。
100 2
68
圆柱公式
体积V = πR 2H
侧面积A' = 2πRH
体积V = πR 2H
侧面积A' = 2πRH
底面积A = πR 2
深搜各层半径和高,外加剪枝。如果确定了一个底层的半径和高度,蛋糕层数,那么往上依次半径-1,高度-1 这样可以算一个最大体积,剩下的体积就不能比它再大了。另外,如果当前表面积已经超过了我们保留的最小值,那么也就不用继续了。
#include<cstdio>
#include<queue>
#include<cmath>
#include<map>
#include<climits>
using namespace std;
int getMaxV(int r,int h,int level)
{
if(level <= 0 || r <= 0 || h <= 0) return 0;
return getMaxV(r-1,h-1,level-1) + r*r*h;
}
int ans = INT_MAX;
int V,L;
void dfs(int lr,int lh,int s,int v,int level)
{
if(s >= ans) return;
if(v < 0 || v > V) return;
if(level == 0)
{
if(v == 0)
ans = s;
return;
}
for(int r = lr-1; r >= level; r--)
for(int h = lh-1; h >= level; h--)
{
if(getMaxV(r,h,level) < v) continue;
int s_ ;
if(level == L) s_ = r*r+2*r*h;
else s_ = s+2*r*h;
int v_ = v - r*r*h;
dfs(r,h,s_,v_,level-1);
}
}
int main()
{
int v,level;
scanf("%d%d",&v,&level);
V = v,L = level;
dfs(sqrt(v),v,0,v,level);
if(ans == INT_MAX)
printf("0\n");
else printf("%d\n",ans);
return 0;
}