QAQ
很明显是一道深搜题目
问题是确定R与H的初始范围是个啥子
由于是整数,那么第一层的蛋糕半径和高度起码就是m(很明显啊)
然后第一层蛋糕半径最大是个(sqrt(v/m))这个也很明显,最大的半径是在最小的高度处取得的
同理有Hmax=v/(m*m)
具体范围确定了,最基础的搜索板子打出来
仔细一看,MMP只过了30的数据,这咋办,于是开始了漫长的剪枝之路
1:可行性剪枝:当当前的表面积已经超过当前的最优解时,直接返回(都知道吧)
2:可行性剪枝:如果我们知道当前剩余蛋糕体积,那么就可以估算一下剩余的表面积,于是有剩余的表面积加上当前的表面积超出了最优解,不进行搜索,也就是s+2*v/R+1< ans
3:最优剪枝:如果当前剩余的蛋糕材料太多,即使最大化剩下的蛋糕体积也做不完,那就不用搜索了,
于是我们预处理一个三维数组表示在I层半径为R高度为H的蛋糕最多剩下多少的体积
这样剪枝,就可以AC了!
一些小细节自己领悟吧!
代码如下
#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring>
#include <cmath>
using namespace std;
int ans=0x7fffffff;
int n,m;
int minx[99999];
int maxx[99][99][99];
void dfs(int x,int r,int h,int v,int s)
{
if(ans<=s) return;
if(x==m)
{
if(v==0) ans=min(ans,s);
return;
}
//if(v<minx[x]) return;
if(v>maxx[x][r][h]) return;
for(int i=r-1;i>=m-x;i--)
for(int j=min(v/(i*i),h-1);j>=m-x;j--)
if(s+2*v/i<ans) dfs(x+1,i,j,v-(i*i*j),s+2*i*j);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
minx[i]=i*i*i+minx[i-1];
for(int r=1;r<=sqrt(n);r++)
for(int h=1;h<=n/(r*r);h++)
{
for(int i=m-1;i>=1;i--)
maxx[i][r][h]=maxx[i+1][r][h]+(r-m+i)*(r-m+i)*(h-m+i);
}
for(int i=m;i<=sqrt(double(n/m));i++)
for(int j=m;j<=(n/(m*m));j++)
dfs(1,i,j,n-i*i*j,2*i*j+i*i);
int t=0x7fffffff;
if(ans!=t) printf("%d",ans);
else printf("0");
}