【题目】
题目描述:
在一个 2 2 2 维平面上有两条传送带,每一条传送带可以看成是一条线段。两条传送带分别为线段 AB 和线段 CD。lxhgww 在 AB 上的移动速度为 P,在 CD 上的移动速度为 Q,在平面上的移动速度 R。现在 lxhgww 想从 A 点走到 D 点,他想知道最少需要走多长时间
输入格式:
输入数据第一行是 4 4 4 个整数,表示 A 和 B 的坐标,分别为 Ax,Ay,Bx,By
第二行是 4 4 4 个整数,表示 C 和 D 的坐标,分别为 Cx,Cy,Dx,Dy
第三行是 3 3 3 个整数,分别是 P,Q,R
输出格式:
输出数据为一行,表示 lxhgww 从 A 点走到 D 点的最短时间,保留到小数点后 2 2 2 位
样例数据:
输入
0 0 0 100
100 0 100 100
2 2 1
输出
136.60
提示:
【数据范围】
对于
100
%
100\%
100% 的数据,
1
1
1 ≤ Ax,Ay,Bx,By,Cx,Cy,Dx,Dy ≤
1000
1000
1000,
1
1
1 ≤ P,Q,R ≤
10
10
10
【分析】
其实题目意识就是在 AB 上选一点 E,在 CD 上选一点 F,使得 A E P + E F R + D F Q \frac{AE}{P}+\frac{EF}{R}+\frac{DF}{Q} PAE+REF+QDF 最小。
如果我们固定 E E E 点,相当于就是找一点 F F F 使 E F R + D F Q \frac{EF}{R}+\frac{DF}{Q} REF+QDF 最小。
然后很容易可以看出,后面的那堆是个单峰函数(可以用几何画板之类的画一下,不会证)。
那么就可以用三分来求解了。
同理,对于找 E E E,也可以发现它是单峰的,也可以用三分来找。
就直接用三分套三分啦。
实际操作的时候直接三分 A E A B \frac{AE}{AB} ABAE( D F C D \frac{DF}{CD} CDDF)即可。
【代码】
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define eps 1e-8
using namespace std;
double ax,ay,bx,by,cx,cy,dx,dy,p,q,r,len1,len2;
double dist(double x1,double y1,double x2,double y2)
{
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double check(double r1,double r2)
{
double x1=ax+(bx-ax)*r1;
double y1=ay+(by-ay)*r1;
double x2=dx-(dx-cx)*r2;
double y2=dy-(dy-cy)*r2;
return len1*r1/p+len2*r2/q+dist(x1,y1,x2,y2)/r;
}
double solve2(double r1)
{
double l=0,r=1;
while(r-l>eps)
{
double m1=l+(r-l)/3;
double m2=r-(r-l)/3;
if(check(r1,m1)<check(r1,m2)) r=m2;
else l=m1;
}
return check(r1,l);
}
double solve1()
{
double l=0,r=1;
while(r-l>eps)
{
double m1=l+(r-l)/3;
double m2=r-(r-l)/3;
if(solve2(m1)<solve2(m2)) r=m2;
else l=m1;
}
return solve2(l);
}
int main()
{
scanf("%lf%lf%lf%lf",&ax,&ay,&bx,&by);
scanf("%lf%lf%lf%lf",&cx,&cy,&dx,&dy);
scanf("%lf%lf%lf",&p,&q,&r);
len1=dist(ax,ay,bx,by);
len2=dist(cx,cy,dx,dy);
printf("%.2lf",solve1());
return 0;
}