POJ 1190 生日蛋糕

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 = πR2H
侧面积A’ = 2πRH
底面积A = πR2

---

实在是想不到，先膜poj某神犇。
从最底层dfs，说一下剪枝技巧

1. sums + mins[deep]> best 表示以前的到的deep + 1层到 m 层的表面积加上从顶层到deep层的最小表面积如果都大于了已经得到的best,那么1到deep层是无论半径和高度取何值都是无效的
   sumv + minv[deep] > n同理

2. 2 * (n - sumv) / r + sums >= best 这是该题的精髓，如果没有的话会造成超时，是为了把sumv和sums联系起来，原因如下:
   假设能够得到best时(为什么这样假设呢，因为如果得不到的话那么就已经被第一个剪枝滤去了,所以在第三个剪枝验证时表示已经通过了第一个剪枝的要求)，
   n - sumv = h[1] * r[1] * r[1] + … + h[deep] * r[deep] * r[deep] < h[1] * r[1] * r + … + h[deep] * r[deep] * r （因为r是deep + 1层的半径）
   其中h[1]…h[deep]表示在函数的形参情况下，1到deep层应该取得h值，r[1]同理
   两边同时处以r 再乘以2得 2 * (n - sumv) / r < 2 * (h[1] * r[1] + … + h[deep] * r[deep])
   2 * (n - sumv) / r < best - sums
   2 * (n - sumv) / r + sums < best 成立 ，则可得剪枝条件

神一样的剪枝，贴ac代码

---

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,m,ans=1e8,mnv[31],mns[31];
void dfs(int deep,int lastr,int lasth,int sums,int sumv)
//deep层，上一层半径，上一层高度，体积和，面积和
{
    if(deep==0)
    {
        if(sumv==n&&sums<ans)
            ans=sums;
        return ;
    }
    if(mns[deep]+sums>ans||mnv[deep]+sumv>n||(2*(n-sumv)/lastr+sums>=ans))
        return ;
    for(int i=lastr-1;i>=deep;i--)//枚举半径 
    {
        if(deep==m)
            sums=i*i;
        int mxh=min(lasth-1,(n-mnv[deep-1]-sumv)/(i*i));
        for(int j=mxh;j>=deep;j--)
            dfs(deep-1,i,j,sums+2*i*j,sumv+i*i*j);
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        mnv[i]=mnv[i-1]+i*i*i;
        mns[i]=mns[i-1]+2*i*i;
    }
    dfs(m,sqrt(n)+1,n+1,0,0);
    if(ans==1e8)
        ans=0;
    printf("%d\n",ans);
    return 0;
}