Bzoj 1857: [Scoi2010]傳送帶(三分套三分)

1857: [Scoi2010]傳送帶
Time Limit: 1 Sec Memory Limit: 64 MB
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

/* 三分套三分. 比較顯然的凸形函數. 凸性函數套凸性函數應該仍是凸性函數吧.... 不會證實. 別忘了討論斜率爲0的狀況.... */
#include<iostream>
#include<cstdio>
#include<cmath>
#define MAXN 10001
#define eps 1e-9
using namespace std;
double ax,ay,bx,by,cx,cy,dx,dy,P,Q,R,ans;
double Cola(double x1,double y1,double x2,double y2)
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double sanfenfuben2(double x,double y)//CD無斜率.
{
    double mid,midmid,k,t1,t2,l,r;
    l=min(cy,dy),r=max(cy,dy);
    while(l+eps<r)
    {
        mid=(2*l+r)/3,midmid=(l+2*r)/3;
        t1=Cola(cx,mid,dx,dy)/Q+Cola(cx,mid,x,y)/R;
        t2=Cola(cx,midmid,dx,dy)/Q+Cola(cx,midmid,x,y)/R;
        if(t1<t2) r=midmid;
        else l=mid;
    }
    t1=Cola(cx,l,dx,dy)/Q+Cola(cx,l,x,y)/R;
    return t1;
}
double sanfen2(double x,double y)//CD有斜率.
{
    if(cx==dx) return sanfenfuben2(x,y);
    double mid,midmid,k,t1,t2,l,r,b;
    l=min(cx,dx),r=max(cx,dx);
    k=(cy-dy)/(cx-dx);
    b=cy-k*cx;
    while(l+eps<r)
    {
        mid=(2*l+r)/3,midmid=(l+2*r)/3;
        t1=Cola(mid,mid*k+b,dx,dy)/Q+Cola(mid,mid*k+b,x,y)/R;
        t2=Cola(midmid,midmid*k+b,dx,dy)/Q+Cola(midmid,midmid*k+b,x,y)/R;
        if(t1<t2) r=midmid;
        else l=mid;
    }
    t1=Cola(l,l*k+b,dx,dy)/Q+Cola(l,l*k+b,x,y)/R;
    return t1;
}
void sanfen1()//AB有斜率.
{
    double mid,midmid,k,t1,t2,l,r,b;
    l=min(ax,bx),r=max(ax,bx);
    k=(ay-by)/(ax-bx);
    b=ay-k*ax;
    while(l+eps<r)
    {
        mid=(2*l+r)/3,midmid=(l+2*r)/3;
        t1=sanfen2(mid,k*mid+b)+Cola(mid,mid*k+b,ax,ay)/P;//+Cola(mid,mid*k,dx,dy)/R;
        t2=sanfen2(midmid,midmid*k+b)+Cola(midmid,midmid*k+b,ax,ay)/P;//+Cola(midmid,midmid*k,dx,dy)/R;
        if(t1<t2) r=midmid;
        else l=mid;
    }
    t1=sanfen2(l,l*k+b)+Cola(l,l*k+b,ax,ay)/P;
    ans=min(ans,t1);
}
void sanfenfuben1()//AB無斜率.
{
    double mid,midmid,k,t1,t2,l,r;
    l=min(ay,by),r=max(ay,by);
    while(l+eps<r)
    {
        mid=(2*l+r)/3,midmid=(l+2*r)/3;
        t1=sanfen2(ax,mid)+Cola(ax,mid,ax,ay)/P;//+Cola(mid,mid*k,dx,dy)/R;
        t2=sanfen2(ax,midmid)+Cola(ax,midmid,ax,ay)/P;//+Cola(midmid,midmid*k,dx,dy)/R;
        if(t1<t2) r=midmid;
        else l=mid;
    }
    t1=sanfen2(ax,l)+Cola(ax,l,ax,ay)/P;
    ans=min(ans,t1);
}
int main()
{
    freopen("walk.in","r",stdin);
    freopen("walk.out","w",stdout);
    cin>>ax>>ay>>bx>>by>>cx>>cy>>dx>>dy;
    cin>>P>>Q>>R;
    ans=Cola(ax,ay,dx,dy)/R;
    if(ax!=bx) sanfen1();
    else sanfenfuben1();
    printf("%.2lf",ans);
    return 0;
}