poj1190生日蛋糕（dfs+剪枝）

->题目请戳这里<-

题目大意：一个m层的蛋糕，每层蛋糕均为圆柱体，总体积为n*PI，从下往上的层数满足每一层的半径和高度逐层严格递减。蛋糕最小表面积（不算底面）。（其中m、n以及每层的半径和高度均为整数，设最小表面积为s*PI，输出s）。如果不存在解，输出0。

题目分析：因为所有数据都是整数，层数不超过20层，每层的半径和高度有是单调递减的，因此可以用dfs遍历所有解，求最小的。用dfs逐层枚举该层的r，h。注意剪枝。

剪枝1：我们从m层往1层搜，即自下而上地搜，当搜到第i层的时候，考虑到最上一层最小的r和h均为1，因为都是整数，而每层半径和高度单调增加的，所以第i层的半径和高度的下界就是i，上界复杂些，对于半径而言，设搜到第I层的时候剩下体积为curv，用极限法，假设剩下的i层是一个圆柱，那么高度最小是1，那么最大半径就是sqrt（curv），而考虑到每层半径单调递减，所以第i层半径还应小于第i+1层半径减1，综上：第i层半径为第i+1层半径减1和sqrt（curv）之间的最小值，同理，第i层的高度上界也要考虑2个值，取第i+1层高度减1和curv之间的最小值。

剪枝2：搜到第i层的时候，剩下的i-1层每层的体积的最小值是一定的，我们从上往下数，最上面一层最小体积为1*1*1，第二层为2*2*2，第三层为3*3*3...第i层为i * i * i ，所以当我们搜到第i层的时候发现剩余体积curv<=∑i*i*i时，可直接返回。

剪枝3：这个剪枝很关键，也很不好想。利用不等式放缩：搜到第i层时第i + 1层半径为lastr，剩下体积(第1层到第i层体积)curv = r1 * r1 * h1 + r2 * r2 * h2 + ... + ri * ri * hi < r1 * lastr * h1 + r2 * lastr * h2 + ... + ri * lastr * hi = lastr*(r1*h1 + r2*h2 + ... + ri*hi)  == >> 2*curv < lastr * 2 * (r1*h1 + r2*h2 + ... + ri*hi)，而2 * (r1 * h1 + ... + ri * hi)就是剩下的i层的表面积，等于我们这次搜索求出的总表面积-第m层到i+1层的表面积，假设当前我们已经求了一个最小表面积ans了，那么我们现在求的剩下的表面积如果比ans小，那么继续往下搜，否则没必要往下搜了，返回。

即当2*curv < (ans - curs)*lastr的时候，我们才继续往下搜，这样才可能求得比当前ans更小的表面积。

详情请见代码：

#include <iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;

int n,m;
int ans;
const int INF = 0x3f3f3f3f;

void dfs(int curm,int curv,int curs,int lastr,int lasth)
{
    if(curv < 0 || curm < 0)
        return;
    if(curs > ans)
        return;
    if(curv < curm * curm * (curm + 1) * (curm + 1) / 4)
        return;
    //printf("%d\n",lastr * (curs + ans));
    if(2 * curv / lastr >= (ans - curs))
        return;
    //printf("hoho\n");
    if(curm == 0 && curv == 0)
    {
        if(curs < ans)
            ans = curs;
		//printf("curm:%d   curv:%d   curs:%d  lastr:%d  lasth:%d\n",curm,curv,curs,lastr,lasth);
        return;
    }
    int i,j;
    int mxr = (int)sqrt(curv * 1.0 ) > lastr - 1?lastr - 1:(int)sqrt(curv * 1.0 );
    int mnr = curm;
    int mxh = curv > lasth - 1?lasth - 1:curv;
    int mnh = curm;
    for(i = mxr;i >= mnr;i --)//i = mnr;i < mxr;i ++)//
    {
        for(j = mxh;j >= mnh;j --)//j = mnh;j < mxh;j ++)//
        {
            int tmpv = i * i * j;
            int tmps = curs;
            if(curm == m)
                tmps += i * i;
            dfs(curm - 1,curv - tmpv,tmps + 2 * i * j,i,j);
        }
    }
}

int main()
{
    while(scanf("%d%d",&n,&m) != EOF)
    {
        ans = INF;
        dfs(m,n,0,INF,INF);
        if(ans == INF)
            printf("0\n");
        else
            printf("%d\n",ans);
    }
    return 0;
}
//144K	32MS
/*
10000 8
1532
*/