搜索+减枝 AAA—————————————————————————————————————————–
已使用面积+剩余蛋糕最小面积>已求出的最小面积 return
剩余体积<已求出的最小体积 return
剩余体积>已求出的最大体积 return
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,m,S,V;
int bst=0x7f7f7f7f;
void search(int fl,int r,int h,int s,int v)
{
if(fl==0 && v==0){
if(s<bst) bst=s;return;
}
else if(fl==0) return;
for(int i=r-1;i>=fl;i--)
for(int j=h-1;j>=fl;j--)
search(fl-1,i,j,s+2*i*j,v-i*i*j);
}
int main()
{
scanf("%d%d",&n,&m);
//枚举最底层的蛋糕
for(int r=m;r<=sqrt(n);r++)//只有一层的话半径最大,h=1
for(int h=n/(r*r);h>=m;h--)//0+4^2*h4=n,h4最大
{
S=r*r+2*r*h;
V=n-r*r*h;
search(m-1,r,h,S,V);//下一层,当前层的半径,高度,已经使用的面积,剩余的体积
}
cout<<bst<<endl;
}
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,m,S,V,mins[25],minv[25];
int bst=0x7f7f7f7f;
void search(int fl,int r,int h,int s,int v)
{
if(fl==0 && v==0){
if(s<bst) bst=s;return;
}
else if(fl==0) return;
if(s+mins[fl]>bst) return; //剪枝1,最优性减枝
if(v>n)return;
if(v<minv[fl])return;//剪枝2,根据下界进行可行性剪枝
for(int i=r-1;i>=fl;i--)
for(int j=h-1;j>=fl;j--)
search(fl-1,i,j,s+2*i*j,v-i*i*j);
}
int main()
{
scanf("%d%d",&n,&m);
//枚举最底层的蛋糕
for(int i=1;i<=m;i++)
mins[i]=mins[i-1]+2*i*i,minv[i]=minv[i-1]+i*i*i;
for(int r=m;r<=sqrt(n);r++)//只有一层的话半径最大,h=1
for(int h=n/(r*r);h>=m;h--)//0+4^2*h4=n,h4最大
{
S=r*r+2*r*h;
V=n-r*r*h;
search(m-1,r,h,S,V);//下一层,当前层的半径,高度,已经使用的面积,剩余的体积
}
cout<<bst<<endl;
}
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,m,S,V,mins[25],minv[25];
int bst=0x7f7f7f7f;
void search(int fl,int r,int h,int s,int v)//本层,上一层的半径,上一层的高,已经使用的面积,剩余的体积
{
if(fl==0 && v==0){
if(s<bst) bst=s;return;
}
else if(fl==0) return;
if(s+mins[fl]>bst) return; //剪枝1,最优性减枝
if(v<minv[fl])return;//剪枝2,根据下界进行可行性剪枝
int maxv=0;
for(int i=1;i<=fl;i++)//因为最大体积,因此上一层都在是下一层基础上减1
maxv+=(r-i)*(r-i)*(h-i);//一共m层,当前fl,加上本层再往上只能一共fl层
if(v>maxv) return;
for(int i=r-1;i>=fl;i--)
for(int j=h-1;j>=fl;j--)
search(fl-1,i,j,s+2*i*j,v-i*i*j);
}
int main()
{
scanf("%d%d",&n,&m);
//枚举最底层的蛋糕
for(int i=1;i<=m;i++)
mins[i]=mins[i-1]+2*i*i,minv[i]=minv[i-1]+i*i*i;
for(int r=m;r<=sqrt(n);r++)//只有一层的话半径最大,h=1
for(int h=n/(r*r);h>=m;h--)//0+4^2*h4=n,h4最大
{
S=r*r+2*r*h;
V=n-r*r*h;
search(m-1,r,h,S,V);//下一层,当前层的半径,高度,已经使用的面积,剩余的体积
}
cout<<bst<<endl;
}