POJ 1190 生日蛋糕

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 = πR  ²H 
侧面积A' = 2πRH 
底面积A = πR  ² 

   一开始看到这个题的时候完全没有思路，时隔一个多月再回来做，看郭炜老师的视频就是清晰易懂。

   解题思路：确定是深度优先搜索，枚举每一层可能的高度和半径。

  确定搜索范围：不要漫无边际的搜，深搜本身就耗时，确定好边界，即底层蛋糕的最大半径和最大高度。

   搜索顺序： 从底层往上搜，而不是从上往下搜。

                       在同一层进行尝试的时候，半径和高度都是从大到小尝试。

 思路有了，接下来就是考虑如何去剪枝，老师讲了四种方法。

 剪枝1：搭建过程中发现面积已经超过了当前求得最优表面积，return。（最优性剪枝）；

剪枝2：搭建过程中预见到再往上搭建，高度已经无法安排，或者半径无法安排，即高度和半径要越来越小，当剩余的层数大于当前的高度或者半径时就无法搭建（可行性剪枝）；

剪枝3：搭建过程中发现还没有搭建的那些层的体积，按最小搭建都会超过还剩余的体积，则停止搭建（可行性剪枝 ）；

剪枝4：搭建过程中发现那些还没有搭的层的体积，最大也到不了还剩的体积，则停止搭建（可行性剪枝）。

  最难想的就是剪枝了，剪枝想到了，接下来就容易处理了。

详见Code：

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
using namespace std;
int N,M;
int minArea=1<<30;
int area=0;
int minV[30];
int minA[30];
int maxforNRH(int n,int r,int h)
{
    int v=0;
    for(int i=0; i<n; i++)
        v+=(r-i)*(r-i)*(h-i);
    return v;
}
void dfs(int v,int n,int r,int h)
//要用n层去凑体积v,最底层半径不能超过r，高度不能超过h
{
    if(n==0)
    {
        if(v) return ; //搭建到最顶层，体积有剩余，不能成功
        else
        {
            minArea=min(minArea,area);//体积无剩余，取最优解
            return ;
        }
    }
    if(v<=0) return ;
    if(area+minA[n]>minArea)  //剪枝1
        return ;
    if(n>r||n>h)     //剪枝2
        return ;
    if(minV[n]>v)    //剪枝3
        return ;
    if(maxforNRH(n,r,h)<v)  //剪枝4，这个剪枝最强没有过不了
        return ;

    for(int rr=r; rr>=n; rr--)
    {
        if(n==M)//顶面积
            area=rr*rr;
        for(int hh=h; hh>=n; hh--)
        {
            area+=2*rr*hh;
            dfs(v-rr*rr*hh,n-1,rr-1,hh-1);
            area-=2*rr*hh;//回溯
        }
    }
}
int main()
{
    minV[0]=0;
    minA[0]=0;
    //预处理，把搭建m层蛋糕所需最小的体积，表面积求出来
    for(int i=1; i<30; i++)
    {
        minV[i]=minV[i-1]+i*i*i;
        minA[i]=minA[i-1]+2*i*i;
    }
    while(~scanf("%d%d",&N,&M))
    {
        if(minV[M]>N)//连最少的都不够
        {
            printf("0\n");
        }
        else
        {
            //大致确定最底层蛋糕的的高度和半径，也就是上面说的确定搜索的范围
            int maxH=(N-minV[M-1])/(M*M) + 1;
            int maxR=sqrt(double(N-minV[M-1])/M) + 1;
            area=0;
            minArea=1<<30;
            dfs(N,M,maxR,maxH);
            if(minArea==1<<30)//不能搭建成功
                printf("0\n");
           else printf("%d\n",minArea);
        }
    }
    return 0;
}