HDU4717

这道题是看了题解之后才做出来的，之前二分一直错

思路：

1.枚举 某一个秒数的时候 的两两之间的最大距离

2.用三分求出对应的时间

至于为什么用三分，因为二分求出的应该是一个具有单调性的结果，对于这种移动中的两点的距离 呈一个有极限的圆锥曲线

我的理解是

如右边的图，所求的极限点应该为x 如果我们用二分在第一步中求得m 而这个m恰好小于于r点对应的y值，那么下一步二分的区间就变成了[l,m]了哦，永远到不了正确答案x的地方了..所以要三分法

/*Code By Aquariuslt*/
/*HDU 4717 The Moving Points*/
#include<iostream>
#include<stdio.h>
#include<string>
#include<algorithm>
#include<string.h>
#include<set>
#include<vector>
#include<queue>
#include<stack>
#include<limits.h>
#include<math.h>
#define FOR(i,a,b) for(int i=(a);i<(b);i++)
#define FORD(i,a,b) for(int i=(a);i<=(b);i++)
#define REP(i,b) FOR(i,0,b)
#define CLR(a) memset(a,0,sizeof(a))
using namespace std;
#define maxn 10001
#define eps 1e-6
int n;
double x[maxn],y[maxn],vx[maxn],vy[maxn];
double dist(int i,int j,double t) {
    double ret=0.0;
    double x1 = x[i] + vx[i]*t;
    double y1 = y[i] + vy[i]*t;
    double x2 = x[j] + vx[j]*t;
    double y2 = y[j] + vy[j]*t;
    ret = sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
    return ret;
}
double maxdist(double t) {
    double maxtemp=0.0;
    REP(i,n) {
        REP(j,n) {
            if(i!=j) {
                maxtemp=max(maxtemp,dist(i,j,t));
            }
        }
    }
    //printf("maxtemp:%.2lf\n",maxtemp);
    return maxtemp;
}
int main() {
    int T;
    scanf("%d",&T);
    REP(ti,T) {
        scanf("%d",&n);
        REP(i,n) {
            scanf("%lf %lf %lf %lf",&x[i],&y[i],&vx[i],&vy[i]);
        }
        double l=0.0,r=1e9;
        while(r-l>eps) {
            double l1=l+(r-l)/3.0;
            double r1=r-(r-l)/3.0;
            if(maxdist(r1)>maxdist(l1))r=r1;
            else l=l1;
            //printf("l:%.2lf r:%.2lf\n",l,r);
            //printf("dl:%.2lf dr:%.2lf\n",maxdist(l),maxdist(r));
        }
        printf("Case #%d: %.2lf %.2lf\n",ti+1,l,maxdist(l));
    }
    return 0;
}