POJ 2507 Crossed ladders 二分 几何

#include <cstdio>
#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
//感谢博友文章 
int main()
{
	double  x,y,c;
	while(1)
	{
		cin>>x>>y>>c;
		if(cin.eof())
		{
			break;
		}
		double l=0;
		double r=min(x,y);
		double w=(l+r)/2;
		while(1)
		{
			//cout<<l<<" "<<r<<" "<<w<<endl;
			if(abs((c/(sqrt(y*y-w*w)))-((sqrt(x*x-w*w))/(sqrt(y*y-w*w)+sqrt(x*x-w*w))))<0.0000001)//看博友程序加一位小数点 
			{
				break;
			}
			w=(l+r)/2;
			if(((c/(sqrt(y*y-w*w)))-((sqrt(x*x-w*w))/(sqrt(y*y-w*w)+sqrt(x*x-w*w))))>0)
			{
				r=w;
			}else 
			{
				l=w;
			}	
		}
		//cout<<abs((c/(sqrt(y*y-w*w)))-((sqrt(x*x-w*w))/(sqrt(y*y-w*w)+sqrt(x*x-w*w))))<<endl;
		//cout<<sqrt(x*x-w*w)<<" "<<sqrt(y*y-w*w)<<endl;
		//cout<<(c/(sqrt(y*y-w*w)))<<endl;
		//cout<<((sqrt(x*x-w*w))/(sqrt(y*y-w*w)+sqrt(x*x-w*w)))<<endl;
		cout<<setiosflags(ios::fixed)<<setprecision(3)<<w<<endl;
	}
}