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
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <map>
#include <queue>
#include <cmath>
#include <vector>
using namespace std;
int n, m, mins = 2100000000;//选一个足够大的数
void dfs(int r, int h, int s, int v, int d)//r为现在这一层的半径,h为现在这一层的高度,s为现在已有蛋糕的表面积,v为现在已有蛋糕的体积,d为现在层数
{
if (d == m) {
if (v == n) {
mins = s;
}
return;
}
if (s + 2 * (n - v) / r >= mins) return;//如果现有面积加后续的最小面积大于mins,剪掉
if (v + (r - 1)* (r - 1) * (h - 1) * (m - d) < n) return;//如果现有体积加后续最大体积小于所需体积,剪掉
if (v + m - d > n) return; //如果现有体积加后续最小体积大于所需体积,剪掉
for (int i = r - 1; i >= m - d; i--) {
for (int j = h - 1; j >= m - d; j--) {
if (v + i * i * j <= n && s + 2 * i * j < mins)//如果加上这一层体积小于等于所需体积,面积小于现有最小面积
dfs (i, j, s + 2 * i * j, v + i * i * j, d + 1);
}
}
}
int main()
{
#ifndef ONLINE_JUDGE
freopen ("in.txt", "r", stdin);
#endif // ONLINE_JUDGE
//printf ("%d\n", mins);
scanf ("%d%d", &n, &m);
for (int i = m; i * i * m <= n; i++)//从最下面一层开始枚举,先枚举r
for (int j = m; i * i * j <= n; j++) {//再枚举h
if(i * i + 2 * i * j < mins)//如果这一层的表面积小于现在已有的最小表面积
dfs (i, j, i * i + 2 * i * j, i * i * j, 1);
}
printf ("%d\n", mins);
return 0;
}