BZOJ_P1857 [SCOI2010]传送带(三分套三分)

BZOJ传送门

Time Limit: 1 Sec Memory Limit: 64 MB
Submit: 975 Solved: 524
[Submit][Status][Discuss]
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

蒟蒻求助！！！！！！
为毛可以三分！！！！
求解为什么QuQ
求私信评论！

#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
#define eps 1e-3
struct Point{
    double x,y;
    void init(){scanf("%lf%lf",&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};}
};
Point A,B,C,D;double P,Q,R;
double Dis(Point &a,Point &b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
double GetTime(Point &a,Point &b){return Dis(A,a)/P+Dis(a,b)/R+Dis(b,D)/Q;}
double Calc(Point &a){
    Point l,r,l2,r2,mm;double t1,t2;
    l=C,r=D;
    while(Dis(l,r)>eps){
        mm=(r-l)/3.0;
        l2=l+mm,r2=l2+mm;
        t1=GetTime(a,l2),t2=GetTime(a,r2);
        if(t1<t2) r=r2;else l=l2;
    }
    return GetTime(a,l);
}
int main(){
    A.init(),B.init(),C.init(),D.init();
    scanf("%lf%lf%lf",&P,&Q,&R);
    Point l,r,l2,r2,mm;double t1,t2;
    l=A,r=B;
    while(Dis(l,r)>eps){
        mm=(r-l)/3.0;
        l2=l+mm,r2=l2+mm;
        t1=Calc(l2),t2=Calc(r2);
        if(t1<t2) r=r2;else l=l2;
    }
    printf("%.2lf\n",Calc(l));
}