传送带（洛谷-P2571）

题目描述

在一个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

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

输出格式：

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

输入输出样例

输入样例#1：

0 0 0 100
100 0 100 100
2 2 1

输出样例#1：

136.60

思路：

假设从 AB 上一点 E 开始走向 CD ，到达 CD 的 F 点，那么 AE+EF+FD 就是走的总路程

假设 AB 上的 E 点固定，那么对于另一个点 F 来说，显然该点存在最小值，是一个单峰函数，此时可利用三分来求出这个点

同理，将 F 点固定后，对另一个点 E，也存在最小值，因此不能单独的三分 E 点或 F 点，要同时三分

源代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<utility>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<bitset>
#define PI acos(-1.0)
#define INF 0x3f3f3f3f
#define LL long long
#define Pair pair<int,int>
const double EPS = 1E-10;
const int MOD = 1E9+7;
const int N = 10000+5;
const int dx[] = {-1,1,0,0,-1,-1,1,1};
const int dy[] = {0,0,-1,1,-1,1,-1,1};
using namespace std;
struct Node{
    double x,y;
    Node(){}
    Node(double x,double y):x(x),y(y){}
}a,b,c,d;
double p,q,r;
double getDis(double x1,double y1,double x2,double y2){
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double judge(Node e){
    Node left(c.x,c.y);
    Node right(d.x,d.y);
    int times=100;//控制三分次数
    while(times--){
        Node lmid,rmid;
        lmid.x=left.x+(right.x-left.x)/3;
        lmid.y=left.y+(right.y-left.y)/3;
        rmid.x=right.x-(right.x-left.x)/3;
        rmid.y=right.y-(right.y-left.y)/3;
        double dis1=getDis(lmid.x,lmid.y,d.x,d.y)/q+getDis(lmid.x,lmid.y,e.x,e.y)/r;
        double dis2=getDis(rmid.x,rmid.y,d.x,d.y)/q+getDis(rmid.x,rmid.y,e.x,e.y)/r;

        if(dis1<dis2){
            right.x=rmid.x;
            right.y=rmid.y;
        }
        else{
            left.x=lmid.x;
            left.y=lmid.y;
        }
    }
    return getDis(left.x,left.y,d.x,d.y)/q+getDis(left.x,left.y,e.x,e.y)/r+getDis(e.x,e.y,a.x,a.y)/p;

}
int main(){
    scanf("%lf%lf",&a.x,&a.y);
    scanf("%lf%lf",&b.x,&b.y);
    scanf("%lf%lf",&c.x,&c.y);
    scanf("%lf%lf",&d.x,&d.y);
    scanf("%lf%lf%lf",&p,&q,&r);

    Node left(a.x,a.y);
    Node right(b.x,b.y);
    int times=100;//控制三分次数
    while(times--){
        Node lmid,rmid;
        lmid.x=left.x+(right.x-left.x)/3;
        lmid.y=left.y+(right.y-left.y)/3;
        rmid.x=right.x-(right.x-left.x)/3;
        rmid.y=right.y-(right.y-left.y)/3;

        if(judge(lmid)<judge(rmid)){
            right.x=rmid.x;
            right.y=rmid.y;
        }
        else{
            left.x=lmid.x;
            left.y=lmid.y;
        }
    }
    printf("%.2lf",judge(left));
    return 0;
}