Openjudge 6048泰国佛塔(dfs)

描述

泰囧，大家应该都看过。不知道大家有没有注意到宝宝虔诚祈福的寺庙里的金灿灿的佛塔。      

众所周知，泰国佛塔外面是贴金的，而且金箔过一段时间需要修补。泰国经济自上世纪末就不景气，贴金的开支对于寺庙来说也是个比较大的负担。          

这里请同学来帮助寺庙设计下佛塔的模型，以使佛塔的外表面（最下层的下底面除外）的面积最小。从而节省金箔开支。      

这里抽象佛塔形状为，一层层的圆柱体，自底向上直径递减。制作一个体积为Nπ的M层佛塔。设从下往上数第i(1 <= i <= M)层佛塔是半径为Ri, 高度为Hi的圆柱。当i < M时，要求R_i> R_i+1且H_i> H_i+1。令外表面面积为Q=Sπ。      

请编程，对于给出的N和M，找出佛塔的制作方案（适当的R_i和H_i的值），使S最小。（除Q外，以上所有数据皆为正整数）

输入

包括若干行。每行包括两个整数N（N <= 10000）、M(M <= 20)，代表一组测试数据。分别表示待建造的佛塔的体积为Nπ,和佛塔的层数为M。

输出

对于每组测试数据，输出占一行：一个正整数S（若无解则S = 0）。

样例输入

100 227 15

样例输出

680

跟切蛋糕类似～ 

AC代码：

#include <iostream>
#include <cmath>
#include "cstdio"
#include "cstring"
#include "algorithm"
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;
}