1006 Line belt

1006  Line belt 

给出两条传送带的起点到末端的坐标，其中ab为p的速度，cd为q的速度 其他地方为r的速度，求a到d点的最短时间.

三分法，先三分第一条线段，找到一个点，然后根据这个点再三分第二条线段.

#include<iostream>  
#include<cstdio>  
#include<cmath>  
#define eps 1e-15  
using namespace std;  
struct point  
{  
    double x,y;  
};  
point A,B,C,D;  
double P,Q,R;  
int sgn(double a)  
{  
    return (a>eps)-(a<-eps);  
}  
double dist(point a,point b)  
{  
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));  
}  
double calcu_lin(point a,point b,point c,point d)  
{  
    return dist(a,b)/P+dist(b,c)/R+dist(c,d)/Q ;  
}  
point MID(point a,point b)  
{  
    point c;  
    c.x=(a.x+b.x)*0.5;  
    c.y=(a.y+b.y)*0.5;  
    return c;  
}  
double sanfen1(point a,point b,point c,point d,point &mid)  
{  
   point mid1,midmid;  
    point l,rr;l=c;rr=d;  
    double t1,t2;  
   do  
   {  
     mid1=MID(l,rr);  
     midmid=MID(mid1,rr);  
     t1=calcu_lin(a,b,mid1,d);  
     t2=calcu_lin(a,b,midmid,d);  
     if(t1<t2)  
     rr=midmid;  
     else  
     l=mid1;  
   }  
   while(sgn(t1-t2)!=0);  
   mid=mid1;  
   return  t1;  
}  
double sanfen2(point d,point c,point b,point a,point &mid)  
{  
   point mid1,midmid;  
    point l,rr;l=a;rr=b;  
    double t1,t2;  
   do  
   {  
     mid1=MID(l,rr);  
     midmid=MID(mid1,rr);  
     t1=calcu_lin(a,mid1,c,d);  
     t2=calcu_lin(a,midmid,c,d);  
     if(t1<t2)  
     rr=midmid;  
     else  
     l=mid1;  
   }  
   while(sgn(t1-t2)!=0);  
   mid=mid1;  
   return  t1;  
}  
void work()  
{  
    double t1,t2;point mid1,mid2;  
    mid1.x=B.x;mid1.y=B.y;  
    do  
    {  
        //printf("okok\n");  
      t1=sanfen1(A,mid1,C,D,mid2);  
      t2=sanfen2(D,mid2,B,A,mid1);  
    }  
    while(sgn(t1-t2)!=0);  
    printf("%.2lf\n",t1);  
}  
int main()  
{  
    int t;  
    double ff;  
    scanf("%d",&t);  
     while(t--)  
    {  
        scanf("%lf%lf%lf%lf",&A.x,&A.y,&B.x,&B.y);  
        scanf("%lf%lf%lf%lf",&C.x,&C.y,&D.x,&D.y);  
        scanf("%lf%lf%lf",&P,&Q,&R);  
          work();  
  
  
    }  
  
    return 0;  
}