hdu-3400 有质量的三分题

题目链接

恩，比较有质量的一个三分的题，属于嵌套三分。

这个题想了很久，只写出几个独立的方程式。看看解题报告说是三分的题。然后，开始想着怎么个三分法。一直没能想出来(检讨下，没有独立的去分析)。

看看大牛的解题报告才知道怎么个三分法;

这个写的不错。http://hi.baidu.com/1093782566/blog/item/2bb8a017124a3735dd540126.html

这里主要解释一下，如何三分。

总耗时t=X/P+Y/R+Z/Q;
X=sqrt((Xx-Ax)^2+(Xy-Ay)^2)，单调递增函数
Y=sqrt((-Yx+Dx)^2+(Dy-Yy)^2)，单调递减函数
Z=sqrt((Yx-Xx)^2+(Yy-Xy)^2) 凸函数。

由于，Z由X,Y确定，Z先增后减。先三分X，再三分Y。X与x(属于【A，B】)有关，Y与y(属于【C,D】)有关。即可以在三分查找X的基础上嵌套三分查找Y。

代码：

#include<iostream>
#include <math.h>
#include <stdio.h>
using namespace std;
#define INF 0.00001
class point
{
public:
	double x;
	double y;
};
point A,B,C,D;
double P,R,Q;
double cal_dis(point a,point b)
{
	return sqrt((b.y-a.y)*(b.y-a.y)+(b.x-a.x)*(b.x-a.x));
}
double cal_y(point a)
{
	point first=C;point end=D;
	point mid,midmid;
	double t1,t2;
	do
	{
		mid.x=(first.x+end.x)/2.0;mid.y=(first.y+end.y)/2.0;
		midmid.x=(mid.x+end.x)/2.0;midmid.y=(mid.y+end.y)/2.0;
		t1=cal_dis(mid,D)/Q+cal_dis(a,mid)/R;
		t2=cal_dis(midmid,D)/Q+cal_dis(a,midmid)/R;
		if(t1>t2) first=mid;
		else end=midmid;
	}while(fabs(t1-t2)>INF);
	return t1;
}
double cal_t()
{
	point first=A;point end=B;
	point mid,midmid;
	double t1,t2;
	do
	{
		mid.x=(first.x+end.x)/2.0;mid.y=(first.y+end.y)/2.0;
		midmid.x=(mid.x+end.x)/2.0;midmid.y=(mid.y+end.y)/2.0;
		t1=cal_dis(mid,A)/P+cal_y(mid);
		t2=cal_dis(midmid,A)/P+cal_y(midmid);
		if(t1>t2) first=mid;
		else end=midmid;
	}while(fabs(t1-t2)>INF);
	return t1; 
}

int main()
{
	int t;
	cin>>t;
	while(t--)
	{
		cin>>A.x>>A.y>>B.x>>B.y>>C.x>>C.y>>D.x>>D.y>>P>>Q>>R;
		printf("%.2lf\n",cal_t());
	}
	return 0;
}

最近心情波动比较大，加油！！！吖飒~~~~相信自己可以的！！！！