hdu 4717 The Moving Points 三分法(2013 ACM/ICPC Asia Regional Online —— Warmup2 1002)

#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
using namespace std;
const double eps=1e-14;
const int maxn=303;
int n;
struct node{
    double x,y,vx,vy;
}e[maxn*maxn];
double cross(node a,node b,double t)
{
    return (a.x+a.vx*t-b.x-b.vx*t)*(a.x+a.vx*t-b.x-b.vx*t)+(a.y+a.vy*t-b.y-b.vy*t)*(a.y+a.vy*t-b.y-b.vy*t);
}
double find(double t)
{
    double ans=0;
    int i,j;
    for(i=0;i<n;i++)
    for(j=i+1;j<n;j++)
        ans=max(ans,cross(e[i],e[j],t));
    return ans;
}
int main()
{
    int T,tt=0;
    scanf("%d",&T);
    while(T--)
    {
        int i,j,k;
        scanf("%d",&n);
        for(i=0;i<n;i++)
        scanf("%lf%lf%lf%lf",&e[i].x,&e[i].y,&e[i].vx,&e[i].vy);
        double l,r,m1,m2;
        l=0,r=1e8;
        for(i=0;i<100;i++)
        {
            m1=l+(r-l)/3;
            m2=r-(r-l)/3;
            if(find(m1)<find(m2)) r=m2;
            else l=m1;
        }
        printf("Case #%d: %.2lf %.2lf\n",++tt,l,sqrt(find(l)));
    }
    return 0;
}
/*
    每组点和速度，根据变量t，横纵坐标就是关于t的以为函数。
    求两点距离的平方就是下凸函数。fi(x)=ai*t*t+bi*t+ci;
    而求n*(n-1)/2个距离的最大值，s(x)=max{fi(x)},也是下凸函数。
    因而可以用三分法求得结果。
*/