poj1190 神剪枝啊 生日蛋糕 (dfs)

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  ² 

又见神剪枝啊，不加

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;
}