POJ1190 生日蛋糕 强大的dfs剪枝!!

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 = πR2
侧面积A' = 2πRH 

底面积A = πR2 

/*by SilverN*/
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const int mxn=30;
int ms[mxn],mv[mxn];
int n,m;
int ans;
void DFS(int smm,int v,int dep,int h,int r){//表面积 体积 层数 高 半径 
    if(v>n)return;
    if(!dep){
        if(v==n && smm<ans)ans=smm;
        return;
    }
    if(v+mv[dep-1]>n || smm+ms[dep-1]>ans || (n-v)/r*2+smm>=ans)return;//当前体积+最小体积大于要求||当前面积+最小面积大于要求

 //还有就是最重要的。如果就当前体积来说,如果我们选择半径r来计算,算出来的表面积比之前算的还要大,直接return。因为后续半径

肯定是比r要大的。百分百不能满足了。
    for(int i=r-1;i>=dep;--i){//对于深度为dep的一层,r最小也要是dep,最大是r-1,自己算算就知道了
        if(dep==m)smm=i*i;// 如果这是最上面的那一层,我们要单独计算一下顶的面积。
        int mxh=min((n-v-mv[dep-1])/(i*i),h-1); //(总体积-之前用的体积-这一层的最小体积)/(当前的半径),得到这一层最大的高度
        for(int j=mxh;j>=dep;--j){//高度同理,最小也是dep,最大的就是上面求得的
            DFS(smm+2*i*j ,v+i*i*j ,dep-1,j,i);//表面积之前的加上新加一层侧面  体积  深度-1 , 高度 半径
        }
    }
    return;
}
void init(){
    for(int i=1;i<=20;i++){
        ms[i]=ms[i-1]+i*i*2;//最小侧面积 
        mv[i]=mv[i-1]+i*i*i;//最小体积 
    }
    return;
}
int main(){
    init ();
    int i,j;
    scanf("%d%d",&n,&m);
    ans=1e9;
    DFS(0,0,m,n+1,n+1);
    if(ans==1e9)ans=0;
    printf("%d\n",ans);
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值