【三分套三分】传送带

题目

在一个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<cstdio>#include<cmath>#include<algorithm>#define fr first#define se second#define mp(p,q) make_pair(p,q)using namespace std;typedef pair<double,double> 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));}