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外,以上所有数据皆为正整数)
设从下往上数第i(1 <= i <= M)层蛋糕是半径为Ri, 高度为Hi的圆柱。当i < M时,要求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 = πR 2H
侧面积A' = 2πRH
底面积A = πR 2
体积V = πR 2H
侧面积A' = 2πRH
底面积A = πR 2
又见神剪枝啊,不加
2*v/R+tmp>=ans
这句话的话,是无法避免超时的
它的意思是:当前解+剩余资源的最优解>=最优解时,剪枝,有点像a*的估价函数
#include <iostream>
#include <cmath>
using namespace std;
int R,H,n,m,ans,mv[21],tmp;
void dfs(int v,int dep,int R,int H) {
if (dep==0) {
if (v==0 && tmp<ans) ans=tmp;
return;
}
if (v-mv[dep-1]<0 || tmp>=ans || 2*v/R+tmp>=ans) return;
for (int r=R-1;r>=dep;r--) {
int Hm = min(H-1,(v-mv[dep-1])/r/r);
for (int h=Hm;h>=dep;h--)
if ((v-r*r*h)>=0) {
if (dep==m) tmp=r*r;
tmp+=2*r*h;
dfs(v-r*r*h,dep-1,r,h);
tmp-=2*r*h;
if (dep==m) tmp=0;
}
}
}
int main() {
memset(mv,0,sizeof(mv));
for (int i=1;i<=m;i++)
mv[i]=mv[i-1]+i*i*i;
while (~scanf("%d%d",&n,&m)) {
ans=10000000;
dfs(n,m,n+1,n+1);
if (ans==10000000) printf("0\n");
else printf("%d\n",ans);
}
return 0;
}