HDU-2289 Cup

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

const double eps = 1e-9;
const double pi = acos(-1.0);
double r1, r2, H, V;
int t;

inline double volume(double h)
{
    double R = h / H * (r2 - r1) + r1;
    double v = (r1*R + r1*r1 + R*R) * h * pi / 3;
    return v;
}
int main()
{
    while(cin >> t)
    {
        while(t--)
        {
           cin >> r1 >> r2 >> H >> V;
           double l = 0, r = 100;
           while(r - l > eps)
           {
               double mid = (r + l) / 2;
               double vm = volume(mid);
               if(vm >= V) r = mid;
               else l = mid;
           }
           printf("%.6f\n",r);
        }
    }


    return 0;
}

题意：输入t（杭电输入的一种习惯），之后输入一个水杯（圆台（上大下小））的 下底半径r1，上底半径r2，杯高度H，体积V。求水的高度。

题解：二分查找。现在做起来满简单的。一开始又犯了一个低级数学错误（直接把上底r2当成了水面的高度，想想怎么这么简单..看来是敲太久晕了）。求水面上底的半径R = h / H * (r2 - r1) + r1（相似比）。其他就很简单了。（可以r = mid或者r = mid - eps 都过）