Description

Input
Output
Sample Input
3 2 1 0.5 2 0.5 3 4 3 4
Sample Output
1.000
0.750
4.000
以前做这道题的时候,写出了式子,但是一直都算不出来。现在才知道这是要用三分法做的。
因为这个函数不是单增或单减的;
做法:
首先,先根据三角形相似写出式子,然后再(0,h)中进行查找。
#include<stdio.h>
#include<iostream>
using namespace std;
double H,h,D,L;
double equ(double x){
return x+(D*h-D*x)/(H-x);
}
double fd(double l,double r){
//double l=0,r=1000.0;
double mid=0,mmid=0;
while(r-l>1e-10){
mid=(l+r)/2.0;
mmid=(mid+r)/2.0;
if(equ(mid)>equ(mmid))
r=mmid;
else l=mid;
}
return equ(mid)>equ(mmid)? equ(mid):equ(mmid);
}
int main(){
int T;
double ans,t,p;
cin>>T;
while(T--){
cin>>H>>h>>D;
//L=D*h/H;
//因为我这里使用的是以竖面上的投影为x,所以x的最大值是h,所以要从0~h进行判断才行;
ans=fd(0,h);
printf("%.3lf\n",ans);
}
}
我一开始不知道三分的left点和right点,我还以为只要从0到最大值就好了;但是这样是不对的,因为如果每个点都在这个范围内的话,那么最大值肯定不会再那个范围内;
比如说x就会超过它的最大值h;
所以三分的范围就是(0,h);
最后三分之后,返回的值是两者中的最大值,因为它要求的是最长影子长度。
原文地址:http://blog.csdn.net/acmer_hades/article/details/43866713