poj 1190 dfs+神剪枝

生日蛋糕
Time Limit: 1000MS Memory Limit: 10000K
Total Submissions: 14120 Accepted: 4971

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

Source



给定蛋糕的体积,求解蛋糕的最小面积。其中有较高层的r和h都要小于较低一层的限制。

剪枝参考http://www.cnblogs.com/rainydays/p/3523162.html

共有三种


按着上述的剪枝来做了,却还是超时严重甚至本机运行都成很大问题

发现是最主要的两步剪枝没做好,原先TLE的做法在下面代码的注释里,

原先的剪枝是在本层剪上一层的枝,本以为不成问题,估计是状态空间真的太大了

所以以后剪枝还是要就剪本层递归的



#include <iostream>
#include <string.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <numeric>
#include <functional>
#define maxn 5005

using namespace std;
int n,m,s;
int check2[25],check3[25];

void calcu()
{
    memset(check2,0,sizeof check2);
    memset(check3,0,sizeof check3);

    check2[1]=2;
    check3[1]=1;
    for(int i=2;i<25;i+=1){
        check2[i]=check2[i-1]+2*i*i;
        check3[i]=check3[i-1]+i*i*i;
        cout<<check3[i]<<" ";
    }
    cout<<"\n";
}

void dfs(int layer,int r,int h,int v,int nows)
{
    if(layer==m){
        if(v+r*r*h==n&&nows+2*r*h<s){
            s=nows+2*r*h;
        }
        return ;
    }
    else if(n-v-r*r*h<check2[m-layer]) return ; //n-v<check2[m-layer+1]
    else if(s-nows-2*r*h<check3[m-layer]) return ; //s-nows<check3[m-layer+1]
    else if(2*(n-v)/r>s-nows) return ;
    else{
        for(int i=m-layer;i<r;i+=1){
            for(int j=m-layer;j<h;j+=1){
                dfs(layer+1,i,j,v+r*r*h,nows+2*r*h);
            }
        }
    }
}

int main()
{
    s=99999999;
    calcu();
    scanf("%d %d",&n,&m);
    for(int i=(int)sqrt(n)+1;i>=1;i-=1){
        for(int j=n;j>=1;j-=1){
            dfs(1,i,j,0,i*i);
        }
    }
    if(s==99999999) printf("0\n",s);
    else printf("%d\n",s);
    return 0;
}



//s 57 85
//i 3 2
//j 8 7







  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值