qduoj 帅气的HYC切蛋糕

帅气的HYC切蛋糕

发布时间: 2015年11月1日 17:02   最后更新: 2015年12月13日 22:14   时间限制: 1000ms   内存限制: 128M

描述

我们的帅气的HYC获得了蛋糕后，当然要找小伙伴分享啦，

这个蛋糕是一个半径为R,高度为H的圆柱形蛋糕，蛋糕上面平铺一层奶油，由于这个蛋糕太大，于是他找到ltwy帮他切蛋糕。

ltwy说要考考他，如果回答正确就帮忙，否则整个蛋糕就归ltwy了。

HYC犹豫了一下，问是什么问题.

ltwy说题目很简单，将蛋糕切下来一块，分成两个部分，告诉HYC这两部分的上表面面积（平铺奶油的那一面）之比r，问切下来的那一块蛋糕的切面面积。

HYC以为是这样:

想了一下确实很简单。不就是2 *Ｒ*Ｈ吗。

没想到ltwy是这么切的:

问此时HYC心里阴影面积（切面面积）。

输入

R 和 H 都是整数（1 < H < 10, 1 < R < 1000），r 是浮点数, 0 < r < 1，以EOF结束。

输出

输出面积，保留两位小数

样例输入1  复制

663 5 0.025
830 1 0.505

样例输出1

3144.76
1602.32

二分，刚开始一直调不对，样例都过不了，，调了很长时间，最后还是请教了某k君，发现是我求圆面积公式记错了。。。。。。。。。。。。。。。。。。

#include<stdio.h>
#include<math.h>
#define exp 1e-12
double judge(double x);
double m;
double pi;
double h,R,r;
int main()
{
	pi=acos(-1.0);
	
	double le,ri,mid;
	while(scanf("%lf %lf %lf",&h,&R,&r)!=EOF)
	{
		le=0;ri=pi/2;
		while(fabs(ri-le)>=exp)
		{
			mid=(le+ri)/2;
			if(judge(mid)*judge(le)<0)
			ri=mid;
			if(judge(mid)*judge(ri)<0)
			le=mid;
		}		
		printf("%.2lf\n",sin(le)*2*h*R);
	}				
	return 0;
}
double judge(double x)
{
	return (sin(2*x) - 2*x + 2*pi*r/(r+1));
}