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/3*(r*r+r*u+u*u)*h;//返回水圆台的体积
}

double ls(double r,double R,double H,double v)
{
    double left,right,temp,mid;
    left = 0,right = 100;
    while ((right-left)>haha)//二分条件变成了与误差值相比较
    {
        mid = (left+right)/2;
        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)/2;
}

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 0;
}