POJ 1190 生日蛋糕题解

题目地址：http://poj.org/problem?id=1190

一道很有趣的搜索题……主要是剪枝……

我弄了5个剪枝：

1.当前剩余层数>=上层半径，剪掉

2.当前剩余层数>=上层高度，剪掉

3.体积<1，剪掉

4.体积<当前可能的最小体积（我用等差数列求和公式算的……也不知道对不对）

5.神奇剪枝网上看的，当前(剩余体积*2)/上层半径+当前表面积总和>当前搜出的最优解，剪掉

代码很简单，20行

 1 #include <iostream>
 2 #include <cmath>
 3 #include <climits>
 4 #define pow2(k) k*k
 5 #define v(r,h) pow2(r)*h
 6 #define min(a,b) a<b?a:b
 7 using namespace std;
 8 int n,m,s=INT_MAX,beg;
 9 void dfs(int r,int h,int c,int lv,int ans)
10 {
11     if(c>m){if(!lv)s=min(s,ans);return;}if(r<=m-c+1||h<=m-c+1||ans+lv*2/r>s||lv<1||lv<(m-c+2)*(m-c+1)/2)return;
12     for(int i=r-1;i>0;i--)for(int j=1;j<h&&v(i,j)<=lv;j++)dfs(i,j,c+1,lv-v(i,j),ans+2*i*j);
13 }
14 int main()
15 {
16     ios::sync_with_stdio(0);
17     cin>>n>>m,beg=sqrt(n)+1;
18     for(int i=beg;i>=1;i--)for(int j=1;v(i,j)<=n;j++)dfs(i,j,2,n-v(i,j),pow2(i)+2*i*j);
19     cout<<(s==INT_MAX?0:s);
20     return 0;
21 }