Description
在一个2维平面上有两条传送带,每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。FTD在AB上的移动速度为P,在CD上的移动速度为Q,在平面上的移动速度R。现在FTD想从A点走到D点,他想知道最少需要走多长时间
Input
输入数据第一行是4个整数,表示A和B的坐标,分别为Ax,Ay,Bx,By 第二行是4个数,表示C和D的坐标,分别为Cx,Cy,Dx,Dy 第三行是3个整数,分别是P,Q,R
Output
输出数据为一行,表示lxhgww从A点走到D点的最短时间,保留到小数点后2位
Sample Input
0 0 0 100
100 0 100 100
2 2 1
Sample Output
136.60
Data Constraint
对于30%的数据
1<=Ax,Ay,Bx,By,Cx,Cy,Dx,Dy<=10
1<=P,Q,R<=5
对于100%的数据
1<=Ax,Ay,Bx,By,Cx,Cy,Dx,Dy<=1000
1<=P,Q,R<=10
The Solution
Analysis
很显然的,答案的路径就是现在 AB 上走,然后走到 CD 的一点,再走完剩下的
一段,于是就可以先三分在 AB 上的转折点,套个三分来三分 CD 上的转折点,
然后直接算就好了。
CODE
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstring>
#define fo(i,a,b) for (int i=a;i<=b;i++)
using namespace std;
struct note
{
double x,y;
}a,b,c,d;
double P,Q,R;
double Sqr(double x) {return x*x;}
double Dis(note a,note b)
{
return sqrt(Sqr(a.x-b.x)+Sqr(a.y-b.y));
}
double Ask(note l,note r)
{
return Dis(a,l)/P+Dis(l,r)/R+Dis(r,d)/Q;
}
double Calc(note q)
{
note l=c,r=d,m1,m2;
while (abs(l.x-r.x)>0.00001 || abs(l.y-r.y)>0.00001)
{
m1.x=(r.x-l.x)/3+l.x;
m1.y=(r.y-l.y)/3+l.y;
m2.x=(r.x-l.x)/3*2+l.x;
m2.y=(r.y-l.y)/3*2+l.y;
if (Ask(q,m1)<Ask(q,m2)) r=m2;
else l=m1;
}
return Ask(q,l);
}
int main()
{
scanf("%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y);
scanf("%lf%lf%lf%lf",&c.x,&c.y,&d.x,&d.y);
scanf("%lf%lf%lf",&P,&Q,&R);
note l=a,r=b,m1,m2;
while (abs(l.x-r.x)>0.00001 || abs(l.y-r.y)>0.00001)
{
m1.x=(r.x-l.x)/3+l.x;
m1.y=(r.y-l.y)/3+l.y;
m2.x=(r.x-l.x)/3*2+l.x;
m2.y=(r.y-l.y)/3*2+l.y;
if (Calc(m1)<Calc(m2)) r=m2;
else l=m1;
}
printf("%.02lf",Calc(l));
return 0;
}