POJ1190 生日蛋糕(DFS)

N久以前写的一个迭代加深和神奇剪枝的搜索，如今我竟然写不出来，悲哀啊！
出答案很简单，但是要AC要加各种神奇的优化，我都写在代码的注释中了啦。

#include<cstdio>#include<cmath>using namespace std;int best = 0xffffff,n,m,minv[25];void dfs(int i,int ri,int hi,int si,int vi)//剩余层数，上一层半径，上一层高度，当前表面积，剩余蛋糕体积{    if(i == 0)    {        if(si < best&&vi == 0) best = si;        return;    }    if(vi < minv[i]) return;//比至少剩的体积还小，做不成蛋糕    if(vi/(2*ri)+si > best) return;//剩余体积/(2*上一层半径) = 可能的最小侧面积，加上已经有的面积，还比最优解大，那么没有更优的解    for(int j = ri-1; j >= i; j--)    {        if(i == m) si = j*j;//最下一层的半径决定了上表面的面积        for(int k = hi-1; k >= i; k--)            dfs(i-1,j,k,si+2*j*k,vi-j*j*k);    }}int main(){    scanf("%d%d",&n,&m);    for(int i = 1; i <= m; i++)//至少剩余的体积，从上往下第i层的半径高度至少为i    {           minv[i] = minv[i-1] + i*i*i;    }    dfs(m,sqrt(n/m),n/(m*m),0,n);//体积/最小高度=最大底面积，开个根号就是最大地面半；体积/(最小底面积半径^2)=地面最大高度    if(best == 0xffffff)        printf("0\n");    else        printf("%d\n",best);}