【搜索】【NOI1999】生日蛋糕

生日蛋糕

问题描述

    7月17日是Mr.W的生日，ACM-THU为此要制作一个体积为Nπ的M层生日蛋糕，每层都是一个圆柱体。

设从下往上数第i(1<=i<=M)层蛋糕是半径为Ri,高度为Hi的圆柱。当i<m< span="">时，要求Ri>Ri+1且Hi>Hi+1。

由于要在蛋糕上抹奶油，为尽可能节约经费，我们希望蛋糕外表面（最下一层的下底面除外）的面积Q最小。

令Q= Sπ

请编程对给出的N和M，找出蛋糕的制作方案（适当的Ri和Hi的值），使S最小。

（除Q外，以上所有数据皆为正整数）

【输入格式】

    有两行，第一行为N（N<=10000），表示待制作的蛋糕的体积为Nπ；第二行为M(M<=20)，表示蛋糕的层数为M。

【输出格式】

    仅一行，是一个正整数S（若无解则S=0）。

【输入输出样例】

输入文件

100

2

输出文件

68

附：圆柱公式

体积V=πR2H

侧面积A’=2πRH

底面积A=πR2

    这是一道不错的搜索题目，比较考验搜索中剪枝的能力，数据范围也需要一些比较神奇的剪枝，个人想到了两个，比较好实现。应该还可以更快一些的，不过这样就可以过了。

Code：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int n,m,ms=0x7fffffff;
bool f=false;
void dfs(int x,int r,int h,int v,int s)
{
	if (v+r*r*h*(m-x)<n || v+m-x>n) return; //可行性剪枝：如果体积超过了规定值，或是接下来无论如何都会使体积超标，就返回 
	if (s+2*(n-v)/r>ms || (s+2*(n-v)/r==ms && x<m)) return; //最优化剪枝：如果接下来无论怎么放都不会是总面积小于最优值，就返回 
	if (x==m) //完成一次查找 
	{
		if (v==n) //符合标准，更新最优解 
		{
			ms=min(ms,s);
			f=true;
		}
		return;
	}
	for (int i=m-x; i<r; i++)
	  for (int j=m-x; j<h; j++)
		dfs(x+1,i,j,v+i*i*j,s+2*i*j);
}
inline void work()
{
	scanf("%d%d",&n,&m);
	for (int i=m; i*i*m<=n; i++) //确定所有可能的初状态 
	  for (int j=m; i*i*j<=n; j++)
		dfs(1,i,j,i*i*j,2*i*j+i*i); //只需要开始加一次上表面，以后都只需要加侧面积即可 
	if (!f) printf("0\n"); //无解 
	else printf("%d\n",ms);
}
int main()
{
	freopen("cake.in","r",stdin);
	freopen("cake.out","w",stdout); 
	work();
	return 0;
}