【三分套三分】传送带

题目

在一个2维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。FTD在AB上的移动速度为P，在CD上的移动速度为Q，在平面上的移动速度R。现在FTD想从A点走到D点，他想知道最少需要走多长时间在一个2维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。FTD在AB上的移动速度为P，在CD上的移动速度为Q，在平面上的移动速度R。现在FTD想从A点走到D点，他想知道最少需要走多长时间

分析

一开始想用小船过河的方法做
发现只能得到特殊的答案
发现对于每条线段每个点的情况是上凸的函数
于是就上三分套三分了

压行代码

#include
#include
#include
#define fr first
#define se second
#define mp(p,q) make_pair(p,q)
using namespace std;
typedef pair pr;
const double eps=1e-6;
double mid,k[2],b[2],u,v,w;
pr a[4],l,r,p,q;
double Dis(pr x,pr y){
    return sqrt((x.fr-y.fr)*(x.fr-y.fr)+(x.se-y.se)*(x.se-y.se));
}
double Cal(pr x,pr y){
    return Dis(a[0],x)/u+Dis(x,y)/w+Dis(y,a[3])/v;
}
double Jud(pr x){
    pr l,r,p,q;
    for(l=a[2],r=a[3],p=mp((2*l.fr+r.fr)/3,(2*l.se+r.se)/3),q=mp((l.fr+2*r.fr)/3,(l.se+2*r.se)/3);Dis(l,r)>eps;p=mp((2*l.fr+r.fr)/3,(2*l.se+r.se)/3),q=mp((l.fr+2*r.fr)/3,(l.se+2*r.se)/3))Cal(x,p)>Cal(x,q)?l=p:r=q;
    return Cal(x,l);
}
int main(){
    freopen("data.txt","r",stdin);
    for(int i=0;i<4;i++)scanf("%lf%lf",&a[i].fr,&a[i].se);
    scanf("%lf%lf%lf",&u,&v,&w);
    for(l=a[0],r=a[1],p=mp((2*l.fr+r.fr)/3,(2*l.se+r.se)/3),q=mp((l.fr+2*r.fr)/3,(l.se+2*r.se)/3);Dis(l,r)>eps;p=mp((2*l.fr+r.fr)/3,(2*l.se+r.se)/3),q=mp((l.fr+2*r.fr)/3,(l.se+2*r.se)/3))Jud(p)>Jud(q)?l=p:r=q;
    printf("%.2f\n",Jud(l));
}