帅气的HYC切蛋糕

描述

我们的帅气的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

数学渣渣表示很无奈。。。。。。。。。。。。。。。。。。。。。。。。。。。看代码~~

1. #include <stdio.h>
2. #include <stdlib.h>
3. #include <math.h>
4. #define PI acosf(-1)
5. const double E = 1e-8;
6. int isEqual (double a, const double b) {
7.     return (a - b >= -E && a - b <= E);
8. }
9. int main() 
   {
10.     int R, H;
11.     double r;
12.     while(~scanf("%d%d",&R,&H)) {
13.         scanf("%lf",&r);
14.         double xiaxian = 0, shangxian = PI;
15.         double i;//i为角度(i<Pi)
16.         double zuo, you;
17.         you = (2 * PI) / (1 / r + 1);
18.         do{
19.             i = (xiaxian + shangxian) / 2;
20.             zuo = i - sin(i);
21.             if(zuo < you) xiaxian = i;
22.             else if(zuo > you) shangxian = i;
23.         }while(!isEqual(zuo,you));
24.         double S = 2 * R * sin(i / 2) * H;
25.         printf("%.2f\n",S);
26.     }
27.     return 0;
28. }
29. //注：PI 可以预处理成acosf(-1) ；//头文件为<math.h>
30. //definde 的内容通常大写
31. //常量通常不用define，而用const ...定义