hdu-3400 有质量的三分题

来源：互联网发布：福州网络安全教育平台编辑：程序博客网时间：2024/04/28 03:36

题目链接

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

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

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

这个写的不错。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.00001class 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;}

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