思路:说一下最重要的剪枝,如果当前已经使用了v的体积,为了让剩下的表面积最小,最好的办法就是让R尽量大,因为V = πR 2H,A' = 2πRH,A' = V / R * 2 ,最大的R一定是取当前的r,所以得到剩下最小表面积就是(n-v)/r*2,如果(n-v)/r*2 + s >= ans,ans是当前得到的最小表面积,直接剪枝。
AC代码
#include <cstdio>
#include <cmath>
#include <cctype>
#include <algorithm>
#include <cstring>
#include <utility>
#include <string>
#include <iostream>
#include <map>
#include <set>
#include <vector>
#include <queue>
#include <stack>
using namespace std;
#pragma comment(linker, "/STACK:1024000000,1024000000")
#define eps 1e-10
#define inf 0x3f3f3f3f
#define PI pair<int, int>
typedef long long LL;
const int maxn = 1e4 + 5;
int ans, n, m;
void dfs(int r, int h, int lev, int v, int s) { //v-已经使用的体积 s-已经出现的表面积
if(v > n || (lev == m && v != n)) return;
if(lev == m && v == n) {
ans = min(ans, s);
return;
}
if(s + (n-v)/r*2 >= ans) {
return; //剪枝
}
for(int i = r-1; i >= m-lev; --i)
for(int j = h-1; j >= m-lev; --j) {
dfs(i, j, lev+1, v + i*i*j, s + 2*i*j);
}
}
int main() {
while(scanf("%d%d", &n, &m) == 2) {
ans = inf;
//枚举最底层
for(int i = 25; i > 0; --i)
for(int j = 25; j > 0; --j) {
if(n-i*i*j >= 0) {
dfs(i, j, 1, i*i*j, i*i + 2*i*j);
}
}
if(ans == inf) printf("-1\n");
else printf("%d\n", ans);
}
return 0;
}
如有不当之处欢迎指出!