HDU2298

题意：圆台，给出上底半径和下底半径，高和里面的水的体积，求水的高度，下底比上底小。

圆台体积

任意给出水高h可求水面半径r1 = h/H*(R-r)+r; H是 圆台高度，r是下底半径，R是上底半径。

二分水高，注意精度。

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cmath>
#define eps 1e-9
#define pi acos(-1)
using namespace std;

double jg(double r,double R,double h,double H)
{
    double r1 = h/H*(R-r)+r;
    return pi/3*(r*r+r*r1+r1*r1)*h;
}

int main()
{
    int n;cin>>n;
    double r,R,H,V;
    while(n--)
    {
        cin>>r>>R>>H>>V;
        double l = 0,t = H;
        while(fabs(l-t)>eps)
        {
        	//cout<<l<<" "<<r<<endl;
            double mid = (l + t) / 2 ;
            double v = jg(r,R,mid,H);
            if(fabs(V-v)<eps)
            {
                //printf("%.6lf\n",r);
				break;
            }
			else if(v<V)
                l=mid+eps;
            else
                t=mid-eps;
        }
        printf("%.6lf\n",l);
    }
    return 0;
}

据说二分半径的话会卡精度，如果二分出半径r1再乘的话误差可能较大，没有尝试，因为我太菜了