HDU4717

来源：互联网发布：淘宝客招募贴范文编辑：程序博客网时间：2024/05/29 02:09

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

思路：

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-6int 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;}