【SCOI 2010】传送带

【题目】

传送门

题目描述：

在一个 $2$ 维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段 AB 和线段 CD。lxhgww 在 AB 上的移动速度为 P，在 CD 上的移动速度为 Q，在平面上的移动速度 R。现在 lxhgww 想从 A 点走到 D 点，他想知道最少需要走多长时间

输入格式：

输入数据第一行是 $4$ 个整数，表示 A 和 B 的坐标，分别为 Ax,Ay,Bx,By

第二行是 $4$ 个整数，表示 C 和 D 的坐标，分别为 Cx,Cy,Dx,Dy

第三行是 $3$ 个整数，分别是 P,Q,R

输出格式：

输出数据为一行，表示 lxhgww 从 A 点走到 D 点的最短时间，保留到小数点后 $2$ 位

样例数据：

输入
0 0 0 100
100 0 100 100
2 2 1

输出
136.60

提示：

【数据范围】
对于 $100\%$ 的数据，$1$ ≤ Ax,Ay,Bx,By,Cx,Cy,Dx,Dy ≤ $1000$，$1$ ≤ P,Q,R ≤ $10$

【分析】

其实题目意识就是在 AB 上选一点 E，在 CD 上选一点 F，使得 $\frac{AE}{P}+\frac{EF}{R}+\frac{DF}{Q}$ 最小。

如果我们固定 $E$ 点，相当于就是找一点 $F$ 使 $\frac{EF}{R}+\frac{DF}{Q}$ 最小。

然后很容易可以看出，后面的那堆是个单峰函数（可以用几何画板之类的画一下，不会证）。

那么就可以用三分来求解了。

同理，对于找 $E$，也可以发现它是单峰的，也可以用三分来找。

就直接用三分套三分啦。

实际操作的时候直接三分 $\frac{AE}{AB}$（$\frac{DF}{CD}$）即可。

【代码】

#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;
}