ZOJ 3203 Light Bulb

设 人离墙的距离是 x   则 F（x） =(H*x - D*h )/(x - D) + x;

分析易得 是   a/x + x 小于 0 的部分 ， 再加上一个常数。

但是重要的一点是 ： F（x） 的定义域并不是 0到 D 。  而是 0到 D*h /H  .

当超出这个区间时，在墙上部分  (H*x - D*h )/(x - D) 变成负值。无意思。

在三分求函数时一定要注意定义域。

#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <stack>
#include <cstring>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <assert.h>
#include <queue>
#define REP(i,n) for(int i=0;i<n;i++)
#define TR(i,x) for(typeof(x.begin()) i=x.begin();i!=x.end();i++)
#define ALLL(x) x.begin(),x.end()
#define SORT(x) sort(ALLL(x))
#define CLEAR(x) memset(x,0,sizeof(x))
#define FILLL(x,c) memset(x,c,sizeof(x))
using namespace std;
const double eps = 1e-5;
#define LL long long 
#define pb push_back
double  H , h ,D;

double f(double x){
//	 return D-x+H-(H-h)*D/x;
     return  (H*x - D*h )/(x - D) + x;
} 
void solve(){
     double left =0;;
     double right =D/H*h;
     int t = 100;
     while(t--){
     	double lmid = (left +right )/2;
     	double rmid = (lmid +right)/2;
     	 if(f(lmid)< f(rmid)){
     	 	  left = lmid ;
     	 }else{
     	 	right = rmid;
     	 }
     }
     double ans = f(right);
     printf("%.3f\n",ans);
}
int main(){
    int t;
    cin >>t ;
    while(t--){
    	scanf("%lf%lf%lf",&H,&h,&D);
    	solve();
    }
    return 0;
}