SCOI2010 传送带

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

感觉三分法很玄妙啊~单峰函数证不来全靠感性理解。。。（可能是做的题太少= =，反正这种题二分不过样例就三分o(╯□╰)o）

比较有启发的是枚举长度，用相似三角形计算坐标，避免了枚举x or y坐标的一系列问题（感兴趣的去枚举x或者y看看？）

#include<bits/stdc++.h>
using namespace std;
#define Inc(i,L,r) for(register int i=(L);i<=(r);++i)
const double eps = 1e-8;
double p,q,r,ax,ay,bx,by,cx,cy,dx,dy;
double dist(double x,double y){
	return sqrt(x*x+y*y);
}
double Calc(double x1,double y1,double lcd){
	double discd=dist(dx-cx,dy-cy);
	double x2=(discd>eps?cx+(dx-cx)*lcd/discd:cx);
	double y2=(discd>eps?cy+(dy-cy)*lcd/discd:cy);
	return dist(x2-x1,y2-y1)/r+(discd-lcd)/q;
}
double calc(double lab){
	double disab=dist(bx-ax,by-ay);
	double x1=(disab>eps?ax+(bx-ax)*lab/disab:ax);
	double y1=(disab>eps?ay+(by-ay)*lab/disab:ay);
	double L=0,r=dist(dx-cx,dy-cy);
	while(r-L>eps){
		double Mid1=((L*2)+r)/3;
		double Mid2=((L)+r*2)/3;
		if(Calc(x1,y1,Mid1)>Calc(x1,y1,Mid2))L=Mid1;
		else r=Mid2;
	}
	return lab/p+Calc(x1,y1,L);
}
void solv(){
	double L=0,r=dist(bx-ax,by-ay);
	while(r-L>eps){
		double Mid1=(L*2+r)/3;
		double Mid2=((L)+r*2)/3;
		if(calc(Mid1)>calc(Mid2))L=Mid1;
		else r=Mid2;
	}
	printf("%.2lf\n",calc(L));
}
void init(){
	cin>>ax>>ay>>bx>>by;
	cin>>cx>>cy>>dx>>dy;
	cin>>p>>q>>r;
}
int main(){
	init();
	solv();
	return 0;
}