《二分》hdoj acm 4.1.5

#include<stdio.h>
#include<math.h>
double pi=acos(-1.0);
double x,y,l,w;
double cal(double a)
{
	 double s,h;
	 s=l*cos(a)+w*sin(a)-x;
	 h=s*tan(a)+w*cos(a);
	 return h;
}
int main()
{
	double left,right,mid,midmid;
	while(~scanf("%lf%lf%lf%lf",&x,&y,&l,&w))
	{
		left=0.0;
		right=pi/2;
		while(fabs(right-left)>1e-8)
		{
			mid=(left+right)/2;
			midmid=(right+mid)/2;
			if(cal(mid)>=cal(midmid))
				right=midmid;
			else
				left=mid;
		}
		if(cal(mid)<=y)
			printf("yes\n");
		else
			printf("no\n");
	}
	return 0;
}