HDU 3400 Line belt

题目大意

给出两条平行的线段AB， CD，然后一个人在线段AB的A点出发，走向D点，其中，人在线段AB上的速度为P， 在线段CD上的速度为Q，在其他地方的速度为R，求人从A点到D点的最短时间。

解题思路

如上图所示，红线部分为人走的路径。人所用的时间为 T = X / P + Y / Q + Z / R。

然后我们做一个变形，令人在线段AB上花的时间为：F(X) = X / P ，假设X是一个确定值的前提下，人走完Z和Y所花的时间为：G(Y) = Z / R + Y / Q。当X 和Y都确定的情况下，Z也是一个确定值。所以，所求的函数可以写成： T（X，Y) = F(X) + G(Y)。因为T是一个关于X和Y的函数，而G只是一个关于Y的函数，所以，可以用两层嵌套的三分法解这个方程。首先以X为变量，对T进行三分。在求G的值的时候，以Y为变量，对G进行三分，求出G在对应的X下的最小值。这样就能求出T的最小值了。

如何证明这个题能够用三分解？

首先，F(X)函数是一个单调递增的函数，而G(Y)很明显是一个先递减后递增的函数。两个函数叠加，所得的T（X,Y）函数应该也是一个先递减后递增的函数。故可用三分法解之。

/*
    problem:HDU 3400 Line belt
    author:DSL_HN_2002
*/
#include <iostream>
#include <cstdlib>
#include <cmath>
using namespace std;

const int Inf = 2147483647;
const double eps = 1e-4;
struct point {
    double x, y;
}a, b, c, d;
double p, q, R;

inline double dist(point a, point b)
{
    return sqrt(pow(a.x - b.x, 2) + pow(a.y - b.y, 2));
}

double calc_2(point a, point b, point to)
{
    point l = a, r = b;
    double res = (double)Inf, res1, res2;
    do {
      point mid1 = (point) {l.x + (r.x - l.x) / 3., l.y + (r.y - l.y) / 3.},
            mid2 = (point) {r.x - (r.x - l.x) / 3., r.y - (r.y - l.y) / 3.};
      res1 = dist(mid1, to) / R + dist(mid1, b) / q,
      res2 = dist(mid2, to) / R + dist(mid2, b) / q;
      if(res1 < res2) r = mid2, res = min(res, res1);
      else l = mid1, res = min(res, res2);
    }while(fabs(res1 - res2) > eps);
    return res;
}

double calc_1(point a, point b)
{
    point l = a, r = b;
    double res = (double)Inf, res1, res2;
    do {
      point mid1 = (point) {l.x + (r.x - l.x) / 3., l.y + (r.y - l.y) / 3.},
            mid2 = (point) {r.x - (r.x - l.x) / 3., r.y - (r.y - l.y) / 3.};
      res1 = calc_2(c, d, mid1) + dist(mid1, a) / p, 
      res2 = calc_2(c, d, mid2) + dist(mid2, a) / p;
      if(res1 < res2) r = mid2, res = min(res, res1);
      else l = mid1, res = min(res, res2);
    }while(fabs(res1 - res2) > eps);
    return res;
}

int main()
{
    freopen("hdu_3400.in", "r", stdin);
    freopen("hdu_3400.out","w", stdout);
    int T;
    double ans;
    cin >> T;
    while(T--) {
      cin >> a.x >> a.y >> b.x >> b.y;
      cin >> c.x >> c.y >> d.x >> d.y;
      cin >> p >> q >> R;
      ans = calc_1(a, b);
      printf("%.2lf\n", ans);
    }
    return 0;
}