POJ1190-生日蛋糕

描述

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外，以上所有数据皆为正整数）

输入

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

输出

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

样例输入

100
2

样例输出

68

【代码】附带解析

#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
int N,M; 
int minS[110],minV[110];//每层蛋糕的最小面积和最小体积 
int area=0,minarea=1e5;//当前正在堆的蛋糕的面积和最小的面积 
int maxv(int n,int r,int h)
{
    int v=0;
    for(int i=0;i<n;i++)
    {
        v+=(r-i)*(r-i)*(h-i); 
    }
    return v;
}
void dfs(int v,int n,int r,int h)
{
    if(n==0)//已经堆到了最后一层 
    {
        if(v)//体积如果没有堆完 
            return;
        else
        {
            minarea=min(area,minarea);//最小的面积进行替换 
            return;
        }
    }
    if(v<=0)
    return;//体积已经被堆完了
    if(r<n||h<n)
    return;//已经小于诶一层最小的半径和高度了
    if(v<minV[n])
    return;
    if(minS[n]+area>=minarea)
    return;
    if(maxv(n,r,h)<v)//所有层的体积大于v就返回
    return;
    for(int rr=r;rr>=n;rr--)
    {
        if(n==M)
            area=rr*rr;//蛋糕的上表面面积投影到最底下就是最底层的面积
        for(int hh=h;hh>=n;hh--)//对每种情况进行搜索
        {
                area+=2*rr*hh;
                dfs(v-rr*rr*hh,n-1,rr-1,hh-1);
                area-=2*rr*hh;
         } 
    }
}
int main()
{
    int maxR,maxH;
    cin>>N>>M;//输入体积和高度
    minS[0]=0,minV[0]=0;
    for(int i=1;i<=M;i++)
    {
        minS[i]=minS[i-1]+2*i*i;
        minV[i]=minV[i-1]+i*i*i;//求出每层蛋糕所需要的最少的面积和体积    
    }    
    if(minV[M]>N)
    {
        cout<<"0"<<endl;
    }
    maxH=(N-minV[M-1])/M*M+1;
    maxR=sqrt((N-minV[M-1])/M)+1;
    dfs(N,M,maxR,maxH);//用M层去堆N体积的蛋糕
    //并且最大的半径 是maxR，最大的高度是maxH
    cout<<minarea<<endl;
    return 0;
}