bailian.openjudge 1190:生日蛋糕

http://bailian.openjudge.cn/practice/1190/

描述

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外,以上所有数据皆为正整数)

输入

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

输出

仅一行,是一个正整数S(若无解则S = 0)。

样例输入

100
2

样例输出

68

提示

圆柱公式
体积 V = π R 2 H V = πR^2H V=πR2H
侧面积 A ′ = 2 π R H A' = 2πRH A=2πRH
底面积 A = π R 2 A = πR^2 A=πR2

代码

#include<algorithm>
#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;

int N, M;
int best;
int leftMinS[22],leftMinV[22];

void dfs(int now,int sumS,int sumV,int r,int h)
{
    if(now==0)
    {
        if(sumV==N && best>sumS)
            best=sumS;
        return;
    }

    if(sumS+leftMinS[now]>best || sumV+leftMinV[now]>N)
        return;

    /*    
     *  2*r*h=S
     *  r*r*h=V
     *  => V*2/r = S
     *  S + sumS >= best => return;   
    */
    if((2*(N-sumV)/r+sumS)>=best)
        return;

    for(int i=r-1; i>=now; i--)
    {
        if(now==M)
            sumS=i*i;
        int maxH=min((N-leftMinV[now-1]-sumV)/(i*i),h-1);
        for(int j=maxH; j>=now; j--)
            dfs(now-1,sumS+i*j*2,sumV+i*i*j,i,j);
    }
}


int main()
{
    leftMinS[0]=0;
    leftMinV[0]=0;
    for(int i=1; i<22; i++)
    {
        leftMinS[i]=leftMinS[i-1]+i*i*2;
        leftMinV[i]=leftMinV[i-1]+i*i*i;
    }
    
    scanf("%d", &N);getchar();
    scanf("%d", &M);getchar();
    
    best=1<<30;
    dfs(M,0,0,100,10000);
    if(best==1<<30)
		best=0;
	
    printf("%d\n", best);
    
    return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值