HDU 2289 几何+圆台

题意;给定一个圆台杯子的R，r，H和里面水的体积V

求h

二分。。。。。。。。。。。。。。。。。

因为解不出来h。。。。。。。。。。。。

/*
几何
圆台体积
V=1/3*pi*h*(r1*r1+r2*r2+r1*r2)

*/
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<queue>
//#include<map>
#include<math.h>
using namespace std;
typedef long long ll;
//typedef __int64 int64;
const int maxn = 105;
const int inf = 0x7fffffff;
const double pi=3.141592653;
double r,R,H,v;

double solve( double h ){
    double r2=h*(R-r)/H+r;
    return pi*h*( r2*r2+r*r+r2*r )/3.0;
}
const double eps = pow(10.0,-8);
int main(){
    int ca;
    scanf("%d",&ca);
    while( ca-- ){
        scanf("%lf%lf%lf%lf",&r,&R,&H,&v);
        double lo,hi;
        lo=0;
        hi=H;
        double mid;
        while( hi-lo>eps ){
            mid=( lo+hi )/2.0;
            if( solve( mid )==v )
                break;
            if( solve( mid )>v )
                hi=mid-eps;
            else
                lo=mid+eps;
        }
        printf("%.6lf\n",mid);
    }
    return 0;
}