深搜剪枝——生日蛋糕

问题 B: 【例题2】生日蛋糕

时间限制: 1 Sec  内存限制: 128 MB
提交: 14  解决: 8
[提交][状态][讨论版][命题人:add_cy]

题目描述

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外，以上所有数据皆为正整数）

输入

有两行，第一行为N（N <= 10000），表示待制作的蛋糕的体积为N;

第二行为M(M <= 20)，表示蛋糕的层数为M。

输出

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

样例输入

100
2

样例输出

68

提示

圆柱公式

体积V = πR2H

侧面积A' = 2πRH

底面积A = πR2

从下往上搜，枚举搜索面对的状态有：正在搜索蛋糕第dep层，当前外表面面积s，当前体积v，第dep+1层的高度和半径

整个蛋糕的上表面面积之和等于最底层的圆面积，可以在第dep层直接累加到s中。这样在第dep-1层往上的搜索中，只需计算侧面积

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define inf 0x3f3f3f3f
using namespace std;
int n,m,ans=0;
void dfs(int dep,int s,int v,int lasth,int lastr)
{
   if(v==n&&dep==0)
   {
       if(s<ans&&ans)
        ans=s;
        else if(ans==0)ans=s;
        return;
   }
   if(dep!=m&&dep*(lastr-1)*(lastr-1)*(lasth-1)+v<n)//当前体积+未遍历的层的体积无法到达n
    return;
   if(m<=0||v>=n||(s>=ans&&ans))
    return;
   for(int i=lastr-1;i>=dep;i--)
   {
       for(int j=lasth-1;j>=dep;j--)
       {
           int s1=2*i*j;
           int v1=i*i*j;
           if(dep==m)s1+=i*i;
           if(s+2*(n-v)/i>=ans&&ans)continue;//当前侧面积+剩余侧面积>已经求得的最小值
           dfs(dep-1,s+s1,v+v1,j,i);
       }
   }
}
int main()
{
    cin>>n>>m;
    dfs(m,0,0,1000,100);
    cout<<ans<<endl;
    return 0;
}