POJ

来源:互联网 发布:淘宝店铺等级在哪儿看 编辑:程序博客网 时间:2024/06/06 12:29
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
1002
Sample Output
68
Hint
圆柱公式 
体积V = πR 2
侧面积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;}