codevs 1710 生日蛋糕

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

借这个题目来谈谈剪枝。
在我看来，剪枝是一种小技巧，而不是一种方法。搜索的进程可以看作是从树根出发，遍历一棵倒置的树——搜索树的过程。而所谓的剪枝，顾名思义，就是通过某种判断，避免一些不必要的遍历过程，形象的说，就是减去了搜索树中的某些“枝条”，以达到优化的目的。
剪枝的方法有两种，一种叫可行性剪枝，另一种叫最优化剪枝（“上下界剪枝”）。（记忆化搜索可以将中间状态记录下来，避免后面不必要的搜索）
剪枝要保证正确性，准确性和高效性。（不能巨斧砍大树。。）

生日蛋糕这道题的优化方法有如下几种：
1、前缀和优化。用vmin[i]和smin[i]记录前i层及其自身的体积（侧面积）最小情况（即第一层r=1,h=1，第二层r=2,h=2…）的总和，用于后面的剪枝。
2、v+r*r*h+vmin[i+1]<=n 如果当前的总体积+该层往上的所有的总体积的最小和>给定蛋糕的体积，剪掉。（可行性剪枝）
3、s+2*r*h+smin[i+1]< best 如果当前的总表面积+该层往上的所有的表面积的最小和>目前最优解,剪掉。（最优化剪枝）
4、v+r*r*h*(m-i+1)>=n 如果后面的所有层均以当前表面积做成蛋糕，如果此时的体积仍然小于蛋糕的体积，就要剪掉。（可行性剪枝）
这种剪枝方法等价于：假设剩余的体积全部以当前层的半径做成蛋糕，因为当前的半径的最大的，所以得出的表面积是最小的，根据这个“最小表面积”判断表面积是否超过最优解，超出则剪掉。
5、搜索的顺序:r从大往小循环，h从小往大循环。由侧面积s=2*r*h,体积v=r*r*h得：在s中r的增长是线性的，在v中r的增长是平方的。r从大往小，h从小往大，可以更快的到达目标（best）即最优解，之后越多不必要的搜索就会被剪枝，从而提高效率。

代码如下：

#include<cstdio>
#include<cmath>
using namespace std;
int n,m,best=0x7fffff,vmin[25],smin[25];

void dfs(int i,int ri,int hi,int s,int v)//i当前层数，ri,hi上一层的半径和高，s,v之前的面积，体积总和 
{
    for (int r=m-i+1; r<=ri-1; r++) //最小状态，第i层（从上往下）的半径和高度均为i 
        for (int h=hi-1; h>=m-i+1; h--)
        {
            if (i==1) s=r*r;
            if (s+2*r*h+smin[i+1]<best&&v+r*r*h+vmin[i+1]<=n&&v+r*r*h*(m-i+1)>=n)
            {
                if (i==m) best=s+2*r*h;
                else dfs(i+1,r,h,s+2*r*h,v+r*r*h);
            }
        }
}
int main()
{
    scanf("%d%d",&n,&m);
    vmin[0]=0; //第0层的体积和表面积都是0(前缀和技巧？) 
    smin[0]=0;
    for (int i=1; i<=m; i++)//从下往上，m层为最上面的一层 
    {
        vmin[m-i+1]=vmin[m-i+2]+i*i*i;//min[i]表示前i层及其自身的体积（侧面积）的总和，用于剪枝 
        smin[m-i+1]=smin[m-i+2]+2*i*i;//smin只记录侧面积 
    }
    dfs(1,sqrt((n-vmin[2])/m)+1,(n-vmin[2])+1,0,0);//最下面那层的上一层的半径为本层半径+1（高同理）
    printf("%d\n",best);
    return 0;
}