POJ 1190-生日蛋糕（dfs 深度搜索 优化减枝）

生日蛋糕 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 21462   Accepted: 7624 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 = πR 2H 侧面积A' = 2πRH 底面积A = πR 2 Source Noi 99

解题思路：一般的搜索题都是在图上找路的问题，表面上看这一题与在图上找并没有什么关系，其实我们把每层蛋糕的状态当作一个节点，它连着上一个状态（下层）和下一个状态（上层），这样一来不就会形成了一幅图，接下来我们就去找到最后那个 面积最小的节点啦！我们从初始状态出发（剩余n层要搭建）一直搜到0层，不过这样递归的搜下去肯定会TLE，所以我们需要对它 进行剪枝，剪枝包括：
最优性剪枝（就是当我们走到达某个节点的时候，发现就算继续走到终点，答案一定不会比目前发现的最 好答案优越，我们这时就需要放弃这种走法了）

可行性剪枝（就是当我们走到某个节点的时候，我们发现目前的资源不足我们到达终点，所以我们必需马上放弃这种走法）

我们这一题这两种剪枝都会用到，接下来就让我们直接来看代码吧！

AC代码：

#include<stdio.h>
#include<iostream>
#include<math.h>
const int INF = 1e9;
using namespace std;

int N, M, R, H, minerae, erae, s = 0;

int minv(int n) //用来求从此状态搭蛋糕需要的最小体积
{
    if(n == 0)
        return 0;
    else return n*n*n + minv(n - 1);
}

void Dfs(int V, int n, int r, int h) //要用n层去凑体积v，最底层半径不能超过r，高度不能超过h
{
    if(n == 0)
    {
        if(V) return; //其实这里也相当于一个可行性剪枝，因为我们的蛋糕不能有剩余
        else {
            minerae = min(minerae, erae); return;//求出最小表面积放入minerae
        } 
    }
    if(erae + 2 * V / r > minerae) return; //最优性剪枝：最小面积剪枝（剩余体积能形成的最小侧面积）
    if(V - minv(n) < 0) return; //可行性剪枝：最小体积剪枝
    for(int rr = r;rr >= n; rr--)
    {
        for(int hh = h; hh >= n; hh--)
        {
            if(n == M) //一开始我们的面积为最底层蛋糕的上表面积
                erae = rr * rr;
            erae += 2 * rr * hh; //加上每层侧面积
            Dfs(V - rr * rr * hh, n - 1, rr - 1,hh - 1);
            erae -= 2 * rr * hh; //这时我们考虑其他情况，因此我们需要把刚才加上的面积减掉
        }
    }
}

int main()
{
    while(~scanf("%d%d", &N, &M))
    {
        erae = 0;
        minerae = INF;
        R = (int)sqrt((double)N - minv(M - 1)); //最下层的最大半径
        H = R; //最大高度，可以是其他值（准确来说是(N-min(M-1))/(M^2)），自己稍微估计下
        Dfs(N, M , R, H);
        if(minerae == INF) printf("0\n");
        else printf("%d\n",minerae);
    }
    return 0;
}