Description
在一个2维平面上有两条传送带,每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。lxhgww在AB上的移动速度为P,在CD上的移动速度为Q,在平面上的移动速度R。现在lxhgww想从A点走到D点,他想知道最少需要走多长时间
Input
输入数据第一行是4个整数,表示A和B的坐标,分别为Ax,Ay,Bx,By
第二行是4个整数,表示C和D的坐标,分别为Cx,Cy,Dx,Dy
第三行是3个整数,分别是P,Q,R
Output
输出数据为一行,表示lxhgww从A点走到D点的最短时间,保留到小数点后2位
Sample Input
0 0 0 100
100 0 100 100
2 2 1
Sample Output
136.60
Hint
【数据范围】
对于100%的数据,1<= Ax,Ay,Bx,By,Cx,Cy,Dx,Dy<=1000,1<=P,Q,R<=10
感觉三分法很玄妙啊~单峰函数证不来全靠感性理解。。。(可能是做的题太少= =,反正这种题二分不过样例就三分o(╯□╰)o)
比较有启发的是枚举长度,用相似三角形计算坐标,避免了枚举x or y坐标的一系列问题(感兴趣的去枚举x或者y看看?)
#include<bits/stdc++.h>
using namespace std;
#define Inc(i,L,r) for(register int i=(L);i<=(r);++i)
const double eps = 1e-8;
double p,q,r,ax,ay,bx,by,cx,cy,dx,dy;
double dist(double x,double y){
return sqrt(x*x+y*y);
}
double Calc(double x1,double y1,double lcd){
double discd=dist(dx-cx,dy-cy);
double x2=(discd>eps?cx+(dx-cx)*lcd/discd:cx);
double y2=(discd>eps?cy+(dy-cy)*lcd/discd:cy);
return dist(x2-x1,y2-y1)/r+(discd-lcd)/q;
}
double calc(double lab){
double disab=dist(bx-ax,by-ay);
double x1=(disab>eps?ax+(bx-ax)*lab/disab:ax);
double y1=(disab>eps?ay+(by-ay)*lab/disab:ay);
double L=0,r=dist(dx-cx,dy-cy);
while(r-L>eps){
double Mid1=((L*2)+r)/3;
double Mid2=((L)+r*2)/3;
if(Calc(x1,y1,Mid1)>Calc(x1,y1,Mid2))L=Mid1;
else r=Mid2;
}
return lab/p+Calc(x1,y1,L);
}
void solv(){
double L=0,r=dist(bx-ax,by-ay);
while(r-L>eps){
double Mid1=(L*2+r)/3;
double Mid2=((L)+r*2)/3;
if(calc(Mid1)>calc(Mid2))L=Mid1;
else r=Mid2;
}
printf("%.2lf\n",calc(L));
}
void init(){
cin>>ax>>ay>>bx>>by;
cin>>cx>>cy>>dx>>dy;
cin>>p>>q>>r;
}
int main(){
init();
solv();
return 0;
}