1. 程式人生 > >【三分套三分】傳送帶

【三分套三分】傳送帶

題目

在一個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)); }