hdu 4717(三分) The Moving Points

Posted 鲨鱼辣椒

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了hdu 4717(三分) The Moving Points相关的知识,希望对你有一定的参考价值。

链接:http://acm.hdu.edu.cn/showproblem.php?pid=4717

n个点,给出初始坐标和移动的速度和移动的方向,求在哪一时刻任意两点之间的距离的最大值的最小。

对于最大值的最小问题首先就会想到二分,然而由于两点之间的距离是呈抛物线状,所以三分的方法又浮上脑海,当然二分也可以做,不过思维上更复杂一些

 1 #include<cmath>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<iostream>
 5 #include<algorithm>
 6 using namespace std;
 7 
 8 const double eps = 0.000001;
 9 const int M = 310;
10 int n;
11 
12 struct point{
13     double x,y,vx,vy;
14     void read(){scanf("%lf%lf%lf%lf",&x,&y,&vx,&vy);}
15 }po[M];
16 
17 double max(double x,double y){return x>y?x:y;}
18 
19 double dis(double x1,double y1,double x2,double y2)
20 {
21     return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
22 }
23 
24 double cal(double t)
25 {
26     double ans=0;
27     for (int i=1 ; i<=n ; i++){
28         for (int j=i+1 ; j<=n ; j++){
29             double num=dis(po[i].x+t*po[i].vx,po[i].y+t*po[i].vy,po[j].x+t*po[j].vx,po[j].y+t*po[j].vy);
30             ans=max(num,ans);
31         }
32     }
33     return ans;
34 }
35 
36 double solve(double l,double r)
37 {
38     while (r-l>=eps)
39     {
40         double mid1=(l*2+r)/3;
41         double mid2=(l+r*2)/3;
42         if (cal(mid2)-cal(mid1)>eps)
43             r=mid2;
44         else
45             l=mid1;
46     }
47     return l;
48 }
49 
50 int main()
51 {
52     int t,cas=0;
53     scanf("%d",&t);
54     while (t--)
55     {
56         scanf("%d",&n);
57         for (int i=1 ; i<=n ; i++)
58             po[i].read();
59         double ans=solve(0,1e8);
60         printf("Case #%d: ", ++cas);
61         if (n==1) printf("0.00 0.00\n");
62         else printf("%.2lf %.2lf\n",ans,cal(ans));
63     }
64     return 0;
65 }

 

以上是关于hdu 4717(三分) The Moving Points的主要内容,如果未能解决你的问题,请参考以下文章

HDU 4355 Party All the Time 三分算法

hdu 2438Turn the corner 三分

三分法hdu2438 Turn the corner

HDU4355-Party All the Time-三分

[hdu4355]Party All the Time(三分)

HDU 4355 Party All the Time(三分|二分)