hdu 2289 要二分的杯子

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2289

大意是 一个Cup,圆台形，给你它的顶部圆的半径，底部圆的半径，杯子的高度，和此时里面装的水的体积，求水的高度。

这道题好在能够锻炼算法思维，或者说形成算法思维。

这是一道二分逼近求值的应用

题目已经给定杯子的高度是0~100，故在这个范围内进行二分，但走了涉及到还有小数的问题，故可以取一个精度较高的误差值

来作为二分的条件，二分一次，以这个值作为水的高度求出水的体积与实际值比较一次，当误差不大于误差值时，便取为高。

另圆台体积公式：设上底的半径为r ,下底的半径为R ,高为h则V＝ (1/3)*pi*h*(R^2 + Rr +r^2) pi为圆周率

 #include <stdio.h>
 #include <math.h>
 double PI=acos(-1.0); //圆周率
 double haha = 0.000000001;//取误差值

 double yj(double r,double R,double h,double H)
 {
     double u = h/H*(R-r) + r; //求出水面圆的半径
     return PI/*(r*r+r*u+u*u)*h;//返回水圆台的体积
 }

 double ls(double r,double R,double H,double v)
 {
     double left,right,temp,mid;
     left = ,right = ;
     while ((right-left)>haha)//二分条件变成了与误差值相比较
     {
         mid = (left+right)/;
         temp= yj(r,R,mid,H);
         if (fabs(temp-v)<=haha)
             return mid;
         else if(temp<v)
             left=mid+haha;
         else
             right = mid-haha;
     }
     return (right+left)/;
 }

 int main()
 {
     int t;
     double r,R,H,v;
     while (~scanf("%d",&t)){
     while (t--)
     {
         scanf("%lf%lf%lf%lf",&r,&R,&H,&v);
         printf("%.6lf\n",ls(r,R,H,v));
     }
     }
     return ;
 }