BZOJ1857: [Scoi2010]传送带(三分套三分)

Time Limit: 1 SecMemory Limit: 64 MB

Submit: 2005Solved: 1091

[Submit][Status][Discuss]

你还在为没有学习平台而苦恼吗？你还在为没有学习资料而烦心吗？你还在为没人指导而担忧吗？可以私信小编 C++或者 加群 710520381 验证灵狐，为你提供学习的平台和资料。

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

Source

Day2

很显然最优的路线一定是从A走到AB上的一点走到CD上的一点再走到D

然后三分套三分就可以了

#include#include#include#defineeps 1e-3usingnamespace std;

inline int read() {

    charc = getchar();intx =0, f =1;

    while(c <'0'|| c >'9') {if(c =='-') f = -1; c = getchar();}

    while(c >='0'&& c <='9') x = x *10+ c -'0',c = getchar();

    returnx * f;

}int Ax, Ay, Bx, By, Cx, Cy, Dx, Dy, P, Q, R;doubledis(doublex,doubley,doublex2,double y2) {

    returnsqrt((x2 - x) * (x2 -x) + (y2 - y) * (y2 - y));

}doublecheck2(doublex,doubley,doublex2,double y2) {

    returndis(Ax, Ay, x, y) / P + dis(x, y, x2, y2) / R + dis(x2, y2, Dx, Dy) / Q;

}doublecheck(doublex,double y) {

    doublelx = Cx, ly = Cy, rx = Dx, ry = Dy, a, b;

    while(abs(rx - lx) > eps || abs(ry - ly) > eps) {

        doublewx1 = (lx *2+ rx) /3, wy1 = (ly *2+ ry) /3,

              wx2 = (lx + rx *2) /3, wy2 = (ly + ry *2) /3;

        a = check2(x, y, wx1, wy1);

        b = check2(x, y, wx2, wy2);

        if(a > b) lx = wx1, ly = wy1;

        elserx = wx2, ry = wy2;

    }

    return check2(x, y, lx, ly);

}int main() {

    #ifdef WIN32

    freopen("a.in","r",stdin);

    #endif    Ax = read(), Ay = read(), Bx = read(), By = read(), Cx = read(), Cy = read(), Dx = read(), Dy = read(), P = read(), Q = read(), R = read();

    doublelx = Ax, ly = Ay, rx = Bx, ry = By;

    while(abs(rx - lx) > eps || abs(ry - ly) > eps) {

        doublewx1 = (lx *2+ rx) /3, wy1 = (ly *2+ ry) /3,

              wx2 = (lx + rx *2) /3, wy2 = (ly + ry *2) /3;

        if(check(wx1, wy1) > check(wx2, wy2)) lx = wx1, ly = wy1;

        elserx = wx2, ry = wy2;

    }

    printf("%.2lf", check(lx, ly));

    return0;

}