关闭

hdu-3400 有质量的三分题

757人阅读 评论(0) 收藏 举报

题目链接

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

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

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

这个写的不错。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;
}
最近心情波动比较大,加油!!!吖飒~~~~相信自己可以的!!!!

2
1

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:106874次
    • 积分:2135
    • 等级:
    • 排名:第18460名
    • 原创:91篇
    • 转载:4篇
    • 译文:0篇
    • 评论:75条
    文章分类
    最新评论