恩,比较有质量的一个三分的题,属于嵌套三分。
这个题想了很久,只写出几个独立的方程式。看看解题报告说是三分的题。然后,开始想着怎么个三分法。一直没能想出来(检讨下,没有独立的去分析)。
看看大牛的解题报告才知道怎么个三分法;
这个写的不错。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;
}
最近心情波动比较大,加油!!!吖飒~~~~相信自己可以的!!!!