BZOJ 1857 SCOI 2010 传送带 三分套三分

分别三分离开AB的点和到达CD的点

传说中的三分套三分。

可以考虑写函数指针节省代码

注意到如果已知AB点到目标点的最小距离，那么这些最小距离组成的是凸函数，对于CD亦是，我们可以采用三分取得极值。

#include <cstdio>
#include <cmath>
const double eps = 1e-4;

struct Point {
	double x, y;
	
	Point operator + (const Point& a) { return (Point) {x+a.x, y+a.y}; }
	Point operator - (const Point& a) { return (Point) {x-a.x, y-a.y}; }
	Point operator * (const double& a) { return (Point) {x*a, y*a}; }
	Point operator / (const double& a) { return (Point) {x/a, y/a}; }
	
	void read() { scanf("%lf%lf", &x, &y); }
	
} A, B, C, D;

double p, q, r;

inline double dis(Point a, Point b) {
	return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

inline double dis2(Point a, Point b) {
	return dis(A, a) / p + dis(a, b) / r + dis(b, D) / q;
}

double calc(Point p0) {
	Point l = C, r = D;
	while (dis(l, r) > eps) {
		Point delta = (r - l) / 3;
		Point l2 = l + delta, r2 = l2 + delta;
		double ret1 = dis2(p0, l2), ret2 = dis2(p0, r2);
		if (ret1 < ret2) r = r2; else l = l2;
	}
	return dis2(p0, l);
}

int main() {
	A.read(); B.read(); C.read(); D.read();
	scanf("%lf%lf%lf", &p, &q, &r);
	
	Point l = A, r = B;
	while (dis(l, r) > eps) {
		Point delta = (r - l) / 3;
		Point l2 = l + delta, r2 = l2 + delta;
		double ret1 = calc(l2), ret2 = calc(r2);
		if (ret1 < ret2) r = r2; else l = l2;
	}
	printf("%.2lf\n", calc(l));
	
	return 0;
}

Description

在一个2维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。lxhgww在AB上的移动速度为P，在CD上的移动速度为Q，在平面上的移动速度R。现在lxhgww想从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

HINT

对于100%的数据，1<= Ax，Ay，Bx，By，Cx，Cy，Dx，Dy<=1000
1<=P，Q，R<=10