POJ 1190 生日蛋糕

http://poj.org/problem?id=1190

这道题 确实写完不知道怎么剪枝了 对代码要求很高

仿一波代码 

剪枝的地方

1、 用mv[]记录最小体积  因为越到上面越小 那么最小的情况就是最高一层 r = 1 ， h = 1 依次递增 但是不确定层数， 就像想一个倒置的蛋糕 这样每次都去掉最大的一块

　　　　　　　　　　　　保证mv[]记录的总是 这个层数 所能达到的最小体积 mv[i] = mv[i-1] + r*r*h//很巧妙

　　　　　　　　　　　　这样 因为总有 r[i] < r[i-1] , h[i] < h[i-1] 底层的蛋糕 > v / 2 ---->>> v > mv[dep-1] 比 这层蛋糕以上所有层的体积都大

2、当前已得表面积 tmp > ans 时 剪枝

3、2*v/R+tmp >= ans ; v/R = H*R ---->>>因为再下一层肯定比这层大 那么如果再加本层的面积 >= ans 也要减去

4、int Hm = min(H-1, (v-mv[dep-1])/ri/ri);  H[i] 满足小于 H[i-1] 并且 最大不超过（体积不超过v-mv[dep-1]） (v-mv[dep-1])/ri/ri;//缩小H的范围

#include <iostream>
#include <stdio.h>
#define INF 0x3f3f3f3f
using namespace std;

int n, m, N, M, ans, tmp;
int mv[28];//记录最小体积

void dfs(int v, int dep, int R, int H)
{
    //从最上层到最下层搜索
    if (dep == 0)
    {
        if (v == 0 && tmp < ans) ans = tmp;
        return ;
    }
    if (v < mv[dep-1] || tmp >= ans || 2*v/R+tmp >= ans) return ;//剪枝
    for(int ri = R-1; ri >= dep; ri--)
    {
        int Hm = min(H-1, (v-mv[dep-1])/ri/ri);
        for (int hi = Hm; hi >= dep; hi--)
        {
            if (v - ri*ri*hi >= 0)
            {
                if (dep == m) tmp = ri*ri;
                tmp += 2*ri*hi;
                dfs(v-ri*ri*hi, dep-1, ri, hi);
                tmp -= 2*ri*hi;
                if (dep == m) tmp = 0;
            }
        }
    }
}
int main()
{
    scanf("%d%d", &n, &m); mv[0] = 0;
    for (int i = 1; i <= m; i++) mv[i] = mv[i-1]+i*i*i;
    ans = INF;
    tmp = 0;
    dfs(n, m, n+1, n+1);
    if (ans == INF) printf("0
");
    else printf("%d
", ans);
    return 0;
}