POJ1190-生日蛋糕

生日蛋糕

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 19089		Accepted: 6794

Description

7月17日是Mr.W的生日，ACM-THU为此要制作一个体积为Nπ的M层生日蛋糕，每层都是一个圆柱体。 
设从下往上数第i(1 <= i <= M)层蛋糕是半径为Ri, 高度为Hi的圆柱。当i < M时，要求Ri > Ri+1且Hi > Hi+1。 
由于要在蛋糕上抹奶油，为尽可能节约经费，我们希望蛋糕外表面（最下一层的下底面除外）的面积Q最小。 
令Q = Sπ 
请编程对给出的N和M，找出蛋糕的制作方案（适当的Ri和Hi的值），使S最小。 
（除Q外，以上所有数据皆为正整数） 

Input

有两行，第一行为N（N <= 10000），表示待制作的蛋糕的体积为Nπ；第二行为M(M <= 20)，表示蛋糕的层数为M。

Output

仅一行，是一个正整数S（若无解则S = 0）。

Sample Input

100
2

Sample Output

68

Hint

圆柱公式 
体积V = πR ²H 
侧面积A' = 2πRH 
底面积A = πR ² 

Source

Noi 99

解题思路：dfs，先进行预处理，处理出从顶层到第i层一共需要的最少的体积和表面积，然后从底层网上推

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#include <queue>
#include <vector>
#include <set>
#include <stack>
#include <map>
#include <climits>
#include <bitset>

using namespace std;

#define LL long long
const int INF=0x3f3f3f3f;

int n,m;
int minv[21],mins[21];
int ans;

void dfs(int v,int s,int l,int r,int h)
{
    if(l==0)
    {
        if(v==n&&s<ans) ans=s;
        return ;
    }
    if(v+minv[l]>n||s+mins[l]>ans||2*(n-v)/r+s>=ans) return ;
    int hh;
    for(int i=r-1;i>=l;i--)
    {
        if(l==m) s=i*i;
        hh=min((n-v-minv[l-1])/(i*i),h-1);
        for(int j=hh;j>=l;j--)
            dfs(v+i*i*j,s+2*i*j,l-1,i,j);
    }
}

int main()
{
    minv[0]=0;
    mins[0]=0;
    for(int i=1;i<=20;i++)
    {
        minv[i]=minv[i-1]+i*i*i;
        mins[i]=mins[i-1]+2*i*i;
    }
    while(~scanf("%d%d",&n,&m))
    {
        ans=INF;
        dfs(0,0,m,n+1,n+1);
        if(ans==INF) printf("0\n");
        else printf("%d\n",ans);
    }
    return 0;
}