P1731 [NOI1999] 生日蛋糕

最新推荐文章于 2024-09-11 09:12:32 发布

SecondBD.

最新推荐文章于 2024-09-11 09:12:32 发布

阅读量158

点赞数 1

文章标签：深度优先图论算法

本文链接：https://blog.csdn.net/2301_78217149/article/details/131028924

版权

[NOI1999] 生日蛋糕

题目背景

数据加强版 link

题目描述

7 月 17 日是 Mr.W 的生日，ACM-THU 为此要制作一个体积为 $N\pi$ 的 $M$ 层生日蛋糕，每层都是一个圆柱体。

设从下往上数第 $i$ （ $\leq i \leq M$ ）层蛋糕是半径为 $R_i$ ，高度为 $H_i$ 的圆柱。当 $\lt M$ 时，要求 $R_i \gt R_{i+1}$ 且 $H_i \gt H_{i+1}$ 。

由于要在蛋糕上抹奶油，为尽可能节约经费，我们希望蛋糕外表面（最下一层的下底面除外）的面积 $Q$ 最小。

请编程对给出的 $N$ 和 $M$ ，找出蛋糕的制作方案（适当的 $R_i$ 和 $H_i$ 的值），使 $S=\dfrac{Q}{\pi}$ 最小。

（除 $Q$ 外，以上所有数据皆为正整数）

输入格式

第一行为一个整数 $N$ （ $\leq 2 \times 10^4$ ），表示待制作的蛋糕的体积为 $N\pi$ 。

第二行为 $M$ （ $\leq 15$ ），表示蛋糕的层数为 $M$ 。

输出格式

输出一个整数 $S$ ，若无解，输出 $0$ 。

样例 #1

样例输入 #1

100
2

样例输出 #1

暴力代码

#include<bits/stdc++.h>
using namespace std;
int r[25],h[25],N,M,minn=0x3f3f3f,s_[25],v_[25];
void dfs(int dep,int s,int v,int ar,int ah)
{
	if (v+v_[dep]>N) return;
	if (s+s_[dep]>minn) return;
	if (s>minn) return;
	if(dep==0)
	{
		if(v==N) minn=min(s,minn);
		return;
	}
	for(int i=ar-1;i>=dep;i--)
		for(int j=ah-1;j>=dep;j--)
		{
			int stmp=2*i*j,vtmp=i*i*j;
			if(M==dep) dfs(dep-1,s+stmp+i*i,v+vtmp,i,j);
			else dfs(dep-1,s+stmp,v+vtmp,i,j);
		}
}
int main()
{
	cin>>N>>M;
	r[M+1]=sqrt(N)+1;
	h[M+1]=N*r[M+1]*r[M+1];
	dfs(M,0,0,r[M+1],h[M+1]);
	cout<<minn;
	return 0;
}//没有任何剪枝

在这里插入图片描述
//九个TLE

AC代码

#include<bits/stdc++.h>
using namespace std;
int R,H,N,M,minn=0x3f3f3f,s_[25],v_[25];
void dfs(int dep,int s,int v,int ar,int ah)//ar记录上一次选的半径，ah记录上一次选的高度 
{
	if (v+v_[dep]>N) return;//体积比N大，剪枝 
	if (s+s_[dep]>minn) return;//表面积比s大，剪枝 
	if (2*(N-v)/ar+s>=minn) return;
	if(dep==0)//到达最上层 
	{
		if(v==N) minn=min(s,minn);//求最小面积 
		return;
	}
	int tmp=sqrt(N-v);
	int imaxd=min(ar-1,tmp);//下一层的一定比上一层的半径小 
	int t_jmaxd=ah-1;//下一层的高度一定比上一层的小 
	for(int i=imaxd;i>=dep;i--)
	{
		int jmaxd=min(t_jmaxd,(N-v)/(i*i));
		for(int j=jmaxd;j>=dep;j--)
		{
			int stmp=2*i*j,vtmp=i*i*j;
			if(M==dep) dfs(dep-1,s+stmp+i*i,v+vtmp,i,j);//如果是最底层，多加一个底面积 
			else dfs(dep-1,s+stmp,v+vtmp,i,j); 
		}
	}
	return;
}
int main()
{
	cin>>N>>M;
	R=sqrt(N)+1;
	H=N*R*R;
	dfs(M,0,0,R,H);//从最底层向上搜 
	if (minn==0x3f3f3f) cout<<-1;没有最小的情况，输出-1
	else cout<<minn;//输出最小面积
	return 0;
}