poj1190 生日蛋糕 （搜索剪枝）

网上很多解题报告，基本一个样，我就不画蛇添足了。

推一下第三个剪枝：如题, 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;
}