网上很多解题报告,基本一个样,我就不画蛇添足了。
推一下第三个剪枝:如题, n为总体积, sumv为从下向上前k-1层的总体积,h为高,r为半径,ans为已经获得的最优的方案
易得:n - sumv = h[ k ] * ( r[ k ] ^ 2 ) + h[ k + 1 ] * ( r[ k + 1 ] ^ 2 )…… + h[ m ] * ( r[ m ] ^ 2 )
而剩余部分的表面积满足:left = 2 * ( r[ k ] * h[ k ] + r[ k + 1 ] * h[ k + 1 ] ……+ r[ m ] * h[ m ] ) > 2 * ( n - sumv ) / r[ k - 1 ]
当ans - sums > left 时,才有可能获得更优的方法,所以 ans - sums > left > 2 * ( n - sumv ) / r[ k - 1 ] 即 2 * ( n - sumv ) / r[ k - 1 ] + sums < ans,
所以当2 * ( n - sumv ) / r[ k - 1 ] + sums >= ans 时也可以进行剪枝.
把代码再贴一遍,留个纪念:
#include <iostream>
#include <cmath>
using namespace std;
#define MIN(a, b) a < b ? a : b
#define M 22
int mins[M], minv[M], n, m, ans;
void dfs(int sumv, int sums, int now, int r, int h)
{
if(now == 0)
{
if(sumv == n && sums < ans)
ans = sums;
return;
}
if(sumv + minv[now - 1] > n || sums + mins[now - 1] > ans || sums + 2 * (n - sumv) / r >= ans)
return;
int i, j, realh;
for(i = r - 1; i >= now; i--)
{
if(now == m)
sums = i * i;
realh = MIN((n - sumv - minv[now - 1]) / (i * i), h - 1);
for(j = realh; j >= now; j--)
dfs(sumv + i * i * j, sums + 2 * i * j, now - 1, i, j);
}
}
int main()
{
int i, maxh, maxr;
scanf("%d%d", &n, &m);
mins[0] = minv[0] = 0;
for(i = 1; i <= m; i++)
{
minv[i] = minv[i-1] + i * i * i;
mins[i] = mins[i-1] + 2 * i * i;
}
maxh = (n - minv[m - 1]) / (m * m);
maxr = sqrt(1.0 *(n - minv[m - 1]) / m);
ans = INT_MAX;
dfs(0, 0, m, maxr + 2, maxh + 2);
if(ans == INT_MAX)
ans = 0;
printf("%d\n", ans);
return 0;
}