ACcoders Problem 2049 题解

题意

在一个二维平面上有两个传送带 $AB$ 和 $CD$，在传送带 $AB$ 上的速度为 $p$，再传送带 $CD$ 上的速度为 $q$，不在传送带上的速度为 $r$，现已给出 $A,B,C,D$ 的坐标以及各自的速度，问从 $A$ 到 $D$ 的最短时间为多少。

思路

我们如果要从 $A$ 跑到 $D$，就绝对会有一部分的路程不在传送带上（如果他们两个重合，当我没说），那么绝对会有一点 $N$ 从传送带 $AB$ 上跑下来，还会有一点 $M$ 从不在传送带的的区域跑上来，所以我们便可以设从 $A$ 到 $D$ 的最短时间为
$$\frac{AN}{p}+\frac{NM}{r}+\frac{MB}{q}$$

两点具体距离可以用欧氏距离来求。

这时，这个式子包含两个未知数，找最小值的时候会变的非常麻烦，那该怎么办呢？

我们考虑将 $N$ 设为一个参数，那么这时候要找的就是
$$f(N)=\frac{NM}{r}+\frac{MB}{q}$$

的最小值。我们可以用软件画一下这个函数（我不会告诉你我其实不会证），很容易发现它是一个单峰函数。

考虑三分 $f(x)$ 的最小值，从而再返回得出整个问题的答案。

现在我们得出来了 $f(N)$，那么这时答案即为 $\frac{AN}{p}+f(N)$，但是我们并不知道 $N$，那该怎么办呢？

我们再用软件画一下此时的函数图像，发现又是一个单峰函数（还是不会证/kk），我们再三分当前函数的最小值，从而得出整个问题的答案（保留 $2$ 为小数）。

代码

#include<bits/stdc++.h>
using namespace std;
double ax,ay,bx,by;
double cx,cy,dx,dy;
double p,q,r;
double dist(double x1,double y1,double x2,double y2){
	return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double f(double x1,double y1,double x2,double y2){
	return dist(x1,y1,x2,y2)/r+dist(x2,y2,dx,dy)/q;
}
double work(double x,double y)
{
	double lx=cx,ly=cy,rx=dx,ry=dy;
	while(dist(lx,ly,rx,ry)>1e-6)
	{
		double midx=(rx-lx)/3.0,midy=(ry-ly)/3.0;
		double lmidx=lx+midx,lmidy=ly+midy;
		double rmidx=rx-midx,rmidy=ry-midy;
		double tmp1=f(x,y,lmidx,lmidy),tmp2=f(x,y,rmidx,rmidy);
		if(tmp1>tmp2) lx=lmidx,ly=lmidy;
		else rx=rmidx,ry=rmidy;
	}
	return f(x,y,lx,ly);
}
double calc()
{
	double lx=ax,ly=ay,rx=bx,ry=by;
	while(dist(lx,ly,rx,ry)>1e-6)
	{
		double midx=(rx-lx)/3.0,midy=(ry-ly)/3.0;
		double lmidx=lx+midx,lmidy=ly+midy;
		double rmidx=lmidx+midx,rmidy=lmidy+midy;
		double tmp1=work(lmidx,lmidy)+dist(ax,ay,lmidx,lmidy)/p;
		double tmp2=work(rmidx,rmidy)+dist(ax,ay,rmidx,rmidy)/p;
		if(tmp1>tmp2) lx=lmidx,ly=lmidy;
		else rx=rmidx,ry=rmidy;
	}
	return work(lx,ly)+dist(ax,ay,lx,ly)/p;
}
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);
	printf("%.2lf",calc());
	return 0;
}