hdu3400 Line belt (双三分)
题意:给出两条平行的线段AB, CD,然后一个人在线段AB的A点出发,走向D点,其中,人在线段AB上的速度为P, 在线段CD上的速度为Q,在其他地方的速度为R,求人从A点到D点的最短时间。
分析:居然要用三分,哎,怎么也没想到,先在AB上三分,固定一个点,接着在CD上三分,就是嵌套的三分搜索,思路清晰的话,代码也不是很难敲
注意精度还有边界的处理
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
struct point
{
double x,y;
};
point A,B,C,D;
int p,q,r;
double length(point X,point Y)
{
return sqrt((X.x-Y.x)*(X.x-Y.x)+(X.y-Y.y)*(X.y-Y.y));
}
double time(double a,double b)
{
point X,Y;
X.x=a*(B.x-A.x)+A.x;
X.y=a*(B.y-A.y)+A.y;
Y.x=b*(C.x-D.x)+D.x;
Y.y=b*(C.y-D.y)+D.y;
return length(A,X)/p+length(D,Y)/q+length(X,Y)/r;
}
double ThiDiv(double lena)
{
double left=0.0,right=1.0,lm,rm;
while(right-left>1e-6)
{
lm=(left*2.0+right)/3.0;
rm=(left+right*2.0)/3.0;
if(time(lena,lm)<=time(lena,rm))
right=rm;
else left=lm;
}
return time(lena,left);
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%lf %lf %lf %lf",&A.x,&A.y,&B.x,&B.y);
scanf("%lf %lf %lf %lf",&C.x,&C.y,&D.x,&D.y);
scanf("%d %d %d",&p,&q,&r);
double left=0.0,right=1.0,lm,rm;
while(right-left>1e-6)//这个大于号居然写出了小于号,测试例子居然也过了,WA了N次
{
//这里三分的是在AB上走的距离占AB比例
lm=(left*2+right)/3.0;
rm=(left+2*right)/3.0;
if(ThiDiv(lm)<=ThiDiv(rm))
right=rm;
else left=lm;
}
printf("%.2f\n",ThiDiv(left));
}
return 0;
}