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

Source

Noi 99

给定蛋糕的体积，求解蛋糕的最小面积。其中有较高层的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