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 = πR2H
侧面积A’ = 2πRH
底面积A = πR2
实在是想不到,先膜poj某神犇。
从最底层dfs,说一下剪枝技巧
sums + mins[deep]> best 表示以前的到的deep + 1层到 m 层的表面积加上从顶层到deep层的最小表面积如果都大于了已经得到的best,那么1到deep层是无论半径和高度取何值都是无效的
sumv + minv[deep] > n同理2 * (n - sumv) / r + sums >= best 这是该题的精髓,如果没有的话会造成超时,是为了把sumv和sums联系起来,原因如下:
假设能够得到best时(为什么这样假设呢,因为如果得不到的话那么就已经被第一个剪枝滤去了,所以在第三个剪枝验证时表示已经通过了第一个剪枝的要求),
n - sumv = h[1] * r[1] * r[1] + … + h[deep] * r[deep] * r[deep] < h[1] * r[1] * r + … + h[deep] * r[deep] * r (因为r是deep + 1层的半径)
其中h[1]…h[deep]表示在函数的形参情况下,1到deep层应该取得h值,r[1]同理
两边同时处以r 再乘以2得 2 * (n - sumv) / r < 2 * (h[1] * r[1] + … + h[deep] * r[deep])
2 * (n - sumv) / r < best - sums
2 * (n - sumv) / r + sums < best 成立 ,则可得剪枝条件
神一样的剪枝,贴ac代码
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,m,ans=1e8,mnv[31],mns[31];
void dfs(int deep,int lastr,int lasth,int sums,int sumv)
//deep层,上一层半径,上一层高度,体积和,面积和
{
if(deep==0)
{
if(sumv==n&&sums<ans)
ans=sums;
return ;
}
if(mns[deep]+sums>ans||mnv[deep]+sumv>n||(2*(n-sumv)/lastr+sums>=ans))
return ;
for(int i=lastr-1;i>=deep;i--)//枚举半径
{
if(deep==m)
sums=i*i;
int mxh=min(lasth-1,(n-mnv[deep-1]-sumv)/(i*i));
for(int j=mxh;j>=deep;j--)
dfs(deep-1,i,j,sums+2*i*j,sumv+i*i*j);
}
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
mnv[i]=mnv[i-1]+i*i*i;
mns[i]=mns[i-1]+2*i*i;
}
dfs(m,sqrt(n)+1,n+1,0,0);
if(ans==1e8)
ans=0;
printf("%d\n",ans);
return 0;
}