要制作一个体积为(Nπ)的(M)层的蛋糕,要求从下往上每一层的半径和高度都是递减的,求一种制作方案使得蛋糕外表面(最下一层的下底面除外)的面积Q最小.令(Q=Sπ),输出最小的S.
反正每个公式里面都有(π),我们干脆把它约掉就好了.然后从下往上开始搜索第x层的(r[x])和(h[x]).
剪枝一:每一层(r[x])上界是(r[x-1]-1)(满足递减),下界是剩余层数(y=m+1-x)(保证第m层的r最小能够取到1),(h[x])上下界同理可得.
剪枝二:从大到小枚举(r[x])和(h[x]),保证状态数能够减少.
剪枝三:如果当前表面积(snow)+底面积(r[1]*r[1])大于等于当前的最小ans,直接return掉.
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<set>
#define ll long long
using namespace std;
inline int read(){
int x=0,o=1;char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')o=-1,ch=getchar();
while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
return x*o;
}
const int N=20;
int n,m,ans,r[N],h[N];
inline void dfs(int x,int y,int vnow,int snow){
if(snow+r[1]*r[1]>=ans)return;
if(vnow==0&&x==m+1){
snow+=r[1]*r[1];
if(snow<ans)ans=snow;
return;
}
if(x==m+1)return;
if(vnow-r[x-1]*r[x-1]*h[x-1]*y>0)return;//无法凑出体积n的蛋糕了
for(int i=r[x-1]-1;i>=y;--i)
for(int j=h[x-1]-1;j>=y;--j){
if(vnow-i*i*j>=0){
r[x]=i;h[x]=j;
dfs(x+1,y-1,vnow-i*i*j,snow+2*i*j);
h[x]=0;r[x]=0;
}
}
}
int main(){
n=read();m=read();
r[0]=(int)sqrt(n);h[0]=(int)sqrt(n);
ans=1e9;dfs(1,m,n,0);
if(ans==1e9)puts("0");
else printf("%d
",ans);
return 0;
}