一本通例题-生日蛋糕——题解<超强深搜剪枝，从无限到有限>

题目传送

显然是道深搜题。由于蛋糕上表面在最底层的半径确认后就确认了，所以搜索时的面积着重看侧面积。

找维度/搜索面临状态/对象：当前体积v，当前外表面面积s，各层的半径r[]，各层的高度h[]。

可行性剪枝考虑/找限制、上下界：

　　　1、考虑当前：当前体积v一定小于总体积N；第i层的半径和高度一定比上一层小（从下往上数层数），同时每次层的高度和半径都>=1(都是正整数)。

　　　2、更近一步，考虑未来：预处理出蛋糕制作到第i层之后再制作的蛋糕体积最小的情况，如果当前体积+这种情况>N，显然不能做成蛋糕;顶层的h、r要大于等于一，而下面每层都要比上面的大，所以h[i]、r[i]>=m-i+1

最优性剪枝考虑：

　　　1、考虑当前：当前面积s小于已搜到的答案ans（否则不会更有，回溯）。

　　　2、更近一步，考虑未来：预处理出蛋糕制作到第i层之后再制作的蛋糕面积最小的情况，如果当前面积积+这种情况>=ans，显然不会更优，回溯。

联系各个维度加强剪枝（简单粗暴的概括：尝试各维度/状态互相表示、将各自的边界融合在一起）：

　1（可行性）：由N-v>=πr²h ，r>=1，h>=1得r<=sqrt(N-v)，h<=(N-v)/r/r（因为最后的Q把π“全包了”，所以可以无视π）（先枚举r，再枚举h）

　　　  2（最优性）：利用h与r数组，dep+1到m层的体积可以表示成n-v=∑（k=dep+1,m）h[k]*r[k]*r[k]，表面积（不算底面积，因为s先前已经算上了）可以表示成2∑（k=dep+1,m）h[k]*r[k]。因为2∑（k=dep+1,m）h[k]*r[k]=2/r[dep]*∑（k=dep+1,m）h[k]*r[k]*r[dep]>=2/r[dep]*∑（k=dep+1,m）h[k]*r[k]*r[k]=2*(n-v)/r[dep], 所以当s+2*(n-v)/r[dep]>=ans时就回溯。

考虑搜索顺序：为了使搜索树较浅的地方分支少，可从底层向顶层搜索。

上AC注释代码:

 #include<iostream>
 #include<cstdio>
 #include<cmath>

 using namespace std;

 int n,m,r[],h[],minv[],mins[];//层数从下往上数 ,mins仅为侧面积
 int s,v,ans=0x7fffffff;//目前面积、体积，最终答案 

 void dfs(int k)//要干第k层
 {
     for(int i=min(r[k-]-,(int)sqrt(n-v));i>=m-k+;i--)//半径r （融合多种剪枝）
         for(int j=min(h[k-]-,(int)((n-v)/(double)i/i));j>=m-k+;j--)//高度h （融合多种剪枝）
         {
             if(v+minv[k]>n||s+mins[k]+r[]*r[]>=ans||*(n-v)/(double)r[k-]+s+r[]*r[]>=ans) return;//可行性&最优性剪枝
             v+=i*i*j;
             s+=*i*j;
             r[k]=i;
             h[k]=j;
             if(k!=m&&v!=n) dfs(k+);
             if(k==m&&v==n) ans=min(ans,s+r[]*r[]);
             v-=i*i*j;
             s-=*i*j;
         }
 }

 int main()
 {
     cin>>n>>m;
     r[]=0x7fffffff;h[]=0x7fffffff;
     for(int i=m;i>=;i--)//预处理
     {
         minv[i]=minv[i+]+(m-i+)*(m-i+);
         mins[i]=mins[i+]+(m-i+)*;
     }
     dfs();
     if(ans==0x7fffffff) ans=;//没有解输出0
     cout<<ans;
     return ;
 }

小总结：仔细分析性质、找状态维度、分析边界、确认好搜索顺序。