POJ 1190 生日蛋糕

题目链接：http://poj.org/problem?id=1190

生日蛋糕

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 17060		Accepted: 6074

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

Source

Noi 99

思路：DFS+剪枝。详见下面我们老师的讲义。

附上AC代码：

#include <cstdio>
#include <cmath>
//#pragma comment(linker, "/STACK:102400000, 102400000")
using namespace std;
const int maxm = 25;
const int inf = 0x3f3f3f3f;
int mins[maxm], mind[maxm];
int n, m, ans;

void init(){
    mins[0] = mind[0] = 0;
    for (int i=1; i<maxm; ++i){
        mins[i] = mins[i-1]+2*i*i;
        mind[i] = mind[i-1]+i*i*i;
    }
}

void dfs(int cur, int r, int h, int s, int d){
    if (cur == m){
        if (d==n && s<ans)
            ans = s;
        return ;
    }
    int ld = n-d;
    for (int i=r-1; i>=m-cur; --i){
        int maxh = (ld-mind[m-cur-1])/(i*i);
        if (maxh > h-1)
            maxh = h-1;
        for (int j=maxh; j>=m-cur; --j){
            int td = ld-i*i*j;
            if (td < mind[m-cur-1])
                continue;
            int ts;
            if (cur == 0)
                ts = i*i+2*i*j;
            else
                ts = s+2*i*j;
            if (ts+mins[m-cur-1]>=ans || ts+2*td/i>=ans)
                continue;
            dfs(cur+1, i, j, ts, d+i*i*j);
        }
    }
}

int main(){
    #ifdef LOCAL
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
    #endif
    init();
    while (~scanf("%d%d", &n, &m)){
        ans = inf;
        int maxh = (n-mind[m-1])/(m*m);
        int maxr = sqrt(1.0*(n-mind[m-1])/m);
        dfs(0, maxh, maxr, 0, 0);
        printf("%d\n", inf==ans ? 0 : ans);
    }
    return 0;
}