题目大意:一个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
*/