DFS.02-生日蛋糕

生日蛋糕

问题描述：要制作一个体积为 N·Π的 M层生日蛋糕，每层都是一个圆柱体
每层的高度和半径都比下一层至少 少1， 由于要在蛋糕上抹奶油，希望蛋糕的外表面最小
根据N和M，找到最小的外表面S

思路

由第 n 层蛋糕的状态 到第 n-1 层蛋糕的状态
枚举 n-1 层蛋糕的半径和高度
针对每个枚举结果，再去枚举 n-2 层蛋糕的状态

第一次写的/错误代码

#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
int N,M;//全局变量 体积N，层数M 
int maxR,maxH;//底层高为1的最大半径，半径为1的最大
int minarea=1<<30;
int area=0;

void dfs(int v,int r,int h,int m){
 if(v<=0)
   return;
 if(m<0)
   return;
 if(m==0)
   minarea=min(area,minarea);
 for(int rr=r;rr>=m;--rr){
   if(m==M)//m是剩余要搭的蛋糕层数 
     area=rr*rr;//相同时正好是整个蛋糕的底面积 
   for(int hh=h;hh>=m;--hh){
      area+=2*rr*hh;
      minarea=min(minarea,area);
      dfs(v-rr*rr*hh,rr-1,hh-1,m-1);
      area-=2*rr*hh;
      }   
   }
}   
int main(){
 cin>>N>>M;
 maxR=sqrt(N);//wrong！
 maxH=N;//wrong！
 cout<<maxR<<maxH<<endl;
 
 dfs(N,maxR,maxH,M);
 if(minarea== 1<<30)
       cout<<0<<endl;
 else
    cout<<minarea<<endl; 
} 

错误原因分析

1.在 maxR 和 maxH 的取值上有问题
搜索范围过大，没有考虑上层蛋糕占用的半径高度，导致结果出问题
2.优化剪枝过差，时间过长

修改方案

1.剪枝优化
1）搭建过程中 发现 area >= minarea 停止搜索
----- if( area > minarea ) return;
2）搭建过程中预判 搭完后 area>=minarea 停止搜索
3）预判再向上搭 高度无法满足，半径无法满足 停止搜索
------rr < n ; hh < n;
4）搭建过程中预判未搭的体积 >剩余的体积，停止搜索
------restv < v - rr * rr * hh
5）搭建过程中预判未搭的体积 <剩余的体积，停止搜索
------restv > v - rr * rr * hh
2.范围缩小

AC代码

#include<iostream>
#include<string.h>
#include<cmath>
using namespace std;
int N,M;
int minArea=1<<30;
int area;
int minV[30];//根据行数储存最小体积 
int minA[30];//最小表面积

int maxforNRH(int m,int r,int h)//！根据剩余的层数来计算最小体积 
{
    int v=0;//重置 
    for(int i=0; i<m; i++)
        v+=(r-i)*(r-i)*(h-i);//从下到上 
    return v;
}
void dfs(int v,int m,int r,int h)
//要用n层去凑体积v,最底层半径不能超过r，高度不能超过h
{
    if(m==0){ 
        if(v) return ; //不能欺骗客户 
        else{ 
            minArea=min(minArea,area);//体积无剩余，取最优解
            return ;
        }
    }
    
    if(v<=0) return ;//不能赔钱 
   
    if(area+minA[m]>minArea)//剪枝1：对第n层**预判**最小面积 
        return ;
   
    if(m>r||m>h)//            剪枝2：半径与高度>=对应层数 
        return ;
   
    if(minV[m]>v)//           剪枝3：对第n层**预判**不可能小于最小体积 
        return ;
   
    if(maxforNRH(m,r,h)<v)//  剪枝4：**预判**体积缺少 
        return ;
    for(int rr=r; rr>=m; rr--)
    {
        if(m==M)//重叠顶面积
            area=rr*rr;
        for(int hh=h; hh>=m; hh--)
        {
            area+=2*rr*hh;
            dfs(v-rr*rr*hh,m-1,rr-1,hh-1);
            area-=2*rr*hh;//回溯
        }
    }
}

int main()
{
    minV[0]=0;
    minA[0]=0;
    for(int i=1; i<M; i++)//每层对应的最小体积，面积 
    {
        minV[i]=minV[i-1]+i*i*i;//从上向下+1
        minA[i]=minA[i-1]+2*i*i;
    }
    while(scanf("%d%d",&N,&M))//多组输入 
    {
        if(minV[M]>N)//连最底层的都不够
            printf("0\n");
        else
        { 
            int maxH=(N-minV[M-1])/(M*M) + 1;
            //除去底层，用剩余的体积/底层的最小半径+底层的最小高度 
            int maxR=sqrt(double(N-minV[M-1])/M) + 1;
           
            area=0;
            minArea=1<<30;
           
            dfs(N,M,maxR,maxH);
            
            if(minArea==1<<30)//不能搭建成功
                printf("0\n");
           else printf("%d\n",minArea);
       }
   }
return 0;
}

这个AC代码不是我自己写的…
但是从这个答案里学到了很多
用数组记录每层的过程值与枚举的状态比较，预判的方法
希望有一天我自己也能写出来，(｡･∀･)ﾉﾞ