Description
在一个2维平面上有两条传送带,每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。FTD在AB上的移动速度为P,在CD上的移动速度为Q,在平面上的移动速度R。现在FTD想从A点走到D点,他想知道最少需要走多长时间
对于100%的数据
1<=Ax,Ay,Bx,By,Cx,Cy,Dx,Dy<=1000
1<=P,Q,R<=10
输出保留两位小数
Analysis
首先,我没什么头绪,想到将线段拆成很多个点,跑最短路。
后来发现根本不需要这样。
首先还是枚举线段AB上的一个出发点P,从A走到P
然后从P走到CD上一点Q,再从Q走到D。
可以发现CD上有一个最优点满足代价最小,因此可以三分CD上的点。
这样就能过了,但是时间太慢了。
再发现,线段AB上的那个点也是可以三分的。
于是变成了喜闻乐见的三分套三分。快得飛飛飛上天。
Code
#include<cstdio>
#include<cmath>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
typedef double db;
const db N=500000,eps=1e-5;
db P,Q,R,k,bb;
struct point
{
db x,y;
}a,b,c,d;
db sqr(db x){return x*x;}
db dis(point a,point b)
{
return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
}
db calc(point p,point p1)
{
return dis(a,p)/P+dis(p,p1)/R+dis(p1,d)/Q;
}
db get(point p)
{
point l=c,r=d,p1,p2;
while(abs(l.x-r.x)>1e-5 || abs(l.y-r.y)>1e-5)
{
p1.x=(l.x+r.x)/2,p1.y=(l.y+r.y)/2;
p2.x=(p1.x+r.x)/2,p2.y=(p1.y+r.y)/2;
if(calc(p,p1)<calc(p,p2)) r=p2;
else l=p1;
}
return calc(p,l);
}
int main()
{
scanf("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y,&d.x,&d.y,&P,&Q,&R);
db ans=N;
point l=a,r=b,p1,p2;
while(abs(l.x-r.x)>1e-5 || abs(l.y-r.y)>1e-5)
{
p1.x=(l.x+r.x)/2,p1.y=(l.y+r.y)/2;
p2.x=(p1.x+r.x)/2,p2.y=(p1.y+r.y)/2;
if(get(p1)<get(p2)) r=p2;
else l=p1;
}
printf("%.02lf",get(l));
return 0;
}