【计算几何】【二分】【随机增量法】hdu6167 Missile Interception

n个半径为R的圆是否有公共部分，等价于询问是否存在一个半径小于R的圆，能覆盖所有n个圆的圆心。

对这n个点求最小圆覆盖即可。从网上扒了个随机增量法的代码。

这样算上二分，复杂度就是nlogn了。

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const double eps=0.000000001;
int n;
double V,x[103],y[103],dx[103],dy[103],v[103],cx[103],cy[103];
struct node  
{  
    double x,y;  
};   
node p[1000001];  
double r;  
node O;  
  
double dist(node a,node b)  
{
    return sqrt( (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y) );  
}
void calc(double a,double b,double c,double d,double e,double f)
{
    O.y=(c*d-f*a)/(b*d-e*a);  
    O.x=(c*e-f*b)/(a*e-b*d);  
}
double get()
{
    for (int i=1;i<=n;++i)   
    {  
        p[i].x=cx[i];
        p[i].y=cy[i];
    }  
     O=p[1];r=0;                
       
     for(int i=2;i<=n;++i)        
     if(dist(O,p[i])>r+1e-6)   
     {  
         O=p[i];r=0;  
         for (int j=1;j<=i-1;++j)      
         if (dist(O,p[j])>r+1e-6)   
         {  
             O.x=(p[i].x+p[j].x)/2;  
             O.y=(p[i].y+p[j].y)/2;  
             r=dist(O,p[j]);  
             for (int k=1;k<=j-1;++k)     
             if (dist(O,p[k])>r+1e-6)  
             {  
                calc(p[j].x-p[i].x,p[j].y-p[i].y,(p[j].x*p[j].x+p[j].y*p[j].y-p[i].x*p[i].x-p[i].y*p[i].y)/2,  
                     p[k].x-p[i].x,p[k].y-p[i].y,(p[k].x*p[k].x+p[k].y*p[k].y-p[i].x*p[i].x-p[i].y*p[i].y)/2);   
               r=dist(O,p[k]);  
             }   
         }  
     }  
     return r;
}
bool check(double t)
{
    for(int i=1;i<=n;++i)
    {
        cx[i]=x[i]+dx[i]*t*v[i];
        cy[i]=y[i]+dy[i]*t*v[i];
    }
    return get()<t*V+eps;
}
int main()
{
    while(scanf("%d%lf",&n,&V)!=EOF)
    {
        for(int i=1;i<=n;++i)
        {
            scanf("%lf%lf%lf%lf%lf",&x[i],&y[i],&dx[i],&dy[i],&v[i]);
            double l=sqrt(dx[i]*dx[i]+dy[i]*dy[i]);
            dx[i]/=l;
            dy[i]/=l;
        }
        double l=0,r=9999999;
        while(l<r-eps)
        {
            double mid=(l+r)*0.5;
            if(check(mid)) r=mid;
            else l=mid;
        }
        printf("%.4lf\n",r);
    }
    return 0;
}