UESTC Training for Search Algorithm——H

积木

Description

madoka小朋友最近在玩一堆积木。这堆积木都是圆柱体。现在madoka想把M个积木堆成M层体积为N的物体。
为了让堆出来的这货看上去比较和谐，madoka规定堆在较下面的圆柱体的底面半径和高度都必须比上面的大。
现在madoka想把堆成物体的外表面涂色（不包括最下面一层底面），madoka想知道怎么来堆积木，才能使总涂色面积S最小
假设所有积木的底面半径和高度均为正整数，且每种积木madoka都有无限个。

Input

有多组数据。每组数据包括两个正整数n,m(n<=10000,m<=20)，m表示堆m层，n表示要求的体积为nπ（π是圆周率）。

Output

对每组数据输出一个正整数s，表示最小面积为sπ（π是圆周率）。如果无解输出0

Sample Input

100 2

Sample Output

68

 

/*算法思想：
  根据数据的范围可以计算出所有的case的高度最大为10000，半径最大为：100，然后DFS
  剪枝1：要是当前还剩的体积小于当前的状态所能得到体积的最小值，这种情况不可能有解
  剪枝2：要是当前的状态所能得到的最小的面积大于作为answer的面积，这种情况不可能有更优的解
*/

/*CODE*/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define inf 100000
#define N 100
using namespace std;
int n,m,ans;
int min(int a,int b)
{
    return a<b?a:b;
}
int min_v(int floor)  //计算还剩floor层是的最小体积
{
    return (floor+1)*(floor+1)*(floor)*(floor)/4;
}
/*DFS，还剩下的层数，还剩下的体积，上一层的半径，上一层的高度，已经有了多大的面积*/
void dfs(int left_floor,int left_v,int last_r,int last_h,int last_s)
{
    if(left_v<min_v(left_floor)) return;  //如果还剩的体积小于当前状态所能得到的最小的体积，退出
    if(last_s+2*left_v/last_r>=ans) return;  //如果当前的状态所能得到的最小面积大于作为answer的面积，退出
    if(left_floor==0)  //已经搜完了所有层
    {
        if(left_v==0) ans=last_s;  //刚好满足体积的要求
        return;
    }
    for(int i=last_r-1;i>=left_floor;i--)  //枚举这一次的半径
    {
        //搜索第一层时面积初始值应该是第一层的圆面的面积，以后搜索的时候，加上当前层的侧面积
        if(left_floor==m) last_s=i*i;
        int maxh=min(last_h-1,(left_v-min_v(left_floor-1))/i/i);
        for(int j=maxh;j>=left_floor;j--)  //枚举这一次的高度
            dfs(left_floor-1,left_v-i*i*j,i,j,last_s+2*i*j);
    }
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        ans=inf;
        dfs(m,n,100,10000,0);  //最大的第一层r为100，最大的第一层h为10000
        if(ans==inf) ans=0;  //没有找到，无解
        printf("%d\n",ans);
    }
}