【三分套三分】傳送帶
阿新 • • 發佈:2019-01-30
題目
在一個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));
}