ZOJ 3203 Light Bulb （三分算法）

/*人左右走动，求影子L的最长长度。
根据图，很容易发现当灯，人的头部和墙角成一条直线时(假设此时人站在A点)，此时的长度是影子全在地上的最长长度。当人再向右走时，影子开始投影到墙上，当人贴着墙，影子长度即为人的高度。所以当人从A点走到墙，函数是先递增再递减，为凸性函数，所以我们可以用三分法来求解。
*/

 

#include<functional>
#include<algorithm>
#include<iostream>
#include<fstream>
#include<sstream>
#include<iomanip>
#include<numeric>
#include<cstring>
#include<cassert>
#include<cstdio>
#include<string>
#include<vector>
#include<bitset>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<list>
#include<set>
#include<map>


using namespace std;


const double eps=1e-9;

int T;

double D,H,h;


double cal(double x)
{
    return D-x+H-(H-h)*D/x;
}


int main()
{
    int T;
    double ans;
    scanf("%d",&T);
    while(T--){
      scanf("%lf%lf%lf",&H,&h,&D);
      double lp=(H-h)*D/H,rp=D;
      while(lp+eps<rp){
        double mid1=(lp*2.0+rp)/3.0;
        double mid2=(lp+rp*2.0)/3.0;
        double val1=cal(mid1);
        double val2=cal(mid2);
        if(val1>=val2) rp=mid2;
        else lp=mid1;
        ans=val1;
      }
      printf("%.3lf\n",ans);
    }
    return 0;
}