三分套三分
证明很恶心。。大概就是确定了在ab上走多长
就可以列出一个关于在cd上走的路程关于时间的函数式,然后这个函数是单峰的
相似的,在ab上走的路程关于时间的函数也是单峰的
所以就可以三分套三分
注意:a和b,c和d可能在一个点上
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
double ax,ay,bx,by,cx,cy,dx,dy,p,q,r,ans,ab,cd;
const double eps=1e-5;
inline double dis(double x1,double y1,double x2,double y2){
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
inline double calc(double x,double y){
double nx1,ny1,nx2,ny2;
if(ab==0) nx1=ax,ny1=ay;
else nx1=(bx-ax)*x/ab+ax,ny1=(by-ay)*x/ab+ay;
if(cd==0) nx2=dx,ny2=dy;
else nx2=(cx-dx)*y/cd+dx,ny2=(cy-dy)*y/cd+dy;
return x/p+y/q+dis(nx1,ny1,nx2,ny2)/r;
}
inline double solve(double a){
double l=0,r=cd,res=0;
do{
double midl=(l*2+r)/3,midr=(l+r*2)/3;
if(calc(a,midl)<=calc(a,midr)) r=midr,res=midl;
else l=midl,res=midr;
}while(r-l>eps);
return calc(a,res);
}
int main(){
scanf("%lf%lf%lf%lf",&ax,&ay,&bx,&by); scanf("%lf%lf%lf%lf",&cx,&cy,&dx,&dy);
scanf("%lf%lf%lf",&p,&q,&r);
ab=dis(ax,ay,bx,by),cd=dis(cx,cy,dx,dy);
double l=0,r=ab;
do{
double midl=(l*2+r)/3,midr=(l+r*2)/3;
if(solve(midl)<=solve(midr)) r=midr,ans=midl;
else l=midl,ans=midr;
}while(r-l>eps);
printf("%.2lf\n",solve(ans));
return 0;
}