题目链接:传送门
题面
题目描述
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,找出蛋糕的制作方案(适当的R
i和Hi的值),使S最小。 (除Q外,以上所有数据皆为正整数)
输入格式
有两行,第一行为N(N <= 10000),表示待制作的蛋糕的体积为Nπ;第二行为M(M <= 20),
表示蛋糕的层数为M。
输出格式
仅一行,是一个正整数S(若无解则S = 0)。
Step 1.爆搜&几个浅显易懂的剪枝
设每层半径分别为
r
1
r_1
r1,
r
2
r_2
r2,…,
r
n
r_n
rn,高度分别为
h
1
h_1
h1,
h
2
h_2
h2,…,
h
n
h_n
hn。
首先可以得出
N
=
∑
r
i
2
h
i
N=\sum r_i^2h_i
N=∑ri2hi
S
=
2
r
i
2
+
2
∑
r
i
h
i
S=\ 2r_i^2 + 2\sum r_ih_i
S= 2ri2+2∑rihi
那么可以先枚举第一层半径和高度,然后依次枚举上面一层的半径和高度。
浅显易懂的剪枝1——设目前表面积为S,令ans保存最终解,若s>ans则可以直接返回。
浅显易懂的剪枝2——设目前体积为V,若V>N则可以直接返回。
Step 2.对于体积V的剪枝
简单地说:
如果剩下的层数均采用最小的r和h时最终的V仍会大于N,则返回。
如果剩下的层数均采用最大的r和h时最终的V仍会小于N,则返回。
然而,均采用最大的r和h时会遇到一个尴尬的问题:
由于
r
i
r_i
ri>
r
i
+
1
r_{i+1}
ri+1,
h
i
h_i
hi>
h
i
+
1
h_{i+1}
hi+1,所以剩下的层数均采用最大的r和h时,增加的体积无法在O(1)的时间内求出!
所以这里要对增加的体积进行放缩,贴个代码:
if(V+(lr-1)*(lr-1)*(lh-1)*(m-deep)<n) return;
其中lr表示上一层的半径,lh表示上一层的高度,m表示总层数,deep表示当前搜索深度。
Step 3.对于表面积S的剪枝
与体积的剪枝类似,当剩下的层数都取最小的表面积时,最终的S仍然大于ans时,返回。
这里同样要对增加的表面积进行放缩,贴个代码:
if(S+(n-V)*2/i>ans) continue;
其中S代表当前表面积,n代表最终要达到的体积,V代表当前体积,i为这一层的半径。
Step 4.Coding!
注:略微卡了下常。
#include<stdio.h>
#include<cstring>
#include<algorithm>
#define re register int
using namespace std;
int read()
{
re x=0,f=1;
char ch=getchar();
while(ch<'0' || ch>'9')
{
if(ch=='-') f=-1;
ch=getchar();
}
while(ch>='0' && ch<='9')
{
x=10*x+ch-'0';
ch=getchar();
}
return x*f;
}
const int Size=1001;
const int INF=0x3f3f3f3f;
int n,m,ans;
int sqr(int x) { return x*x; }
void dfs(int deep,int V,int S,int lr,int lh)
{
if(S>ans || V>n) return;
if(V+sqr((m-deep)*(m-deep+1)>>1)>n) return;
if(V+(lr-1)*(lr-1)*(lh-1)*(m-deep)<n) return;
if(deep==m)
{
if(V==n && S<ans)
{
ans=S;
}
return;
}
for(re i=m-deep; i<lr; i++)
{
if(S+(((n-V)<<1)/i)>ans) continue;
for(re j=m-deep; j<lh; j++)
{
dfs(deep+1,V+i*i*j,S+(i*j<<1),i,j);
}
}
}
int main()
{
n=read();
m=read();
ans=INF;
for(re i=m; i<=100; i++)
{
for(re j=m; j<=200; j++)
{
dfs(1,i*i*j,i*i+2*i*j,i,j);
}
}
printf("%d",ans);
return 0;
}