求一个覆盖所有点的最小圆 最小圆覆盖

题目大意：平面上有n个点，求绘制一个半径最小的圆，覆盖所有的点  精度0.1  点的坐标最大为 100000

 

方法1：http://wenku.baidu.com/view/584b6d3e5727a5e9856a610d.html   O(n)

方法2：三分套三分暴力求解如下。O(1600 n)时间开销  1600~(log  2/3  0.1/100000) ^2

float disy(float y,float x)
{
    int n=pvec.size();
    float Rmin=0;
    
    for(int i=0;i<n;i++)
    {
        Rmin=max(Rmin,sqrt((pvec[i].y-y)*(pvec[i].y-y)+ (pvec[i].x-x)*(pvec[i].x-x)));
    }
    return Rmin;
}
float dis(float x)
{
    ymax=10000,ymin=0;
    while(ymax-ymin>accu)
    {
        if(disy(ymin+(ymax-ymin)*2.0/3,x)<disy(ymin+(ymax-ymin)*1.0/3,x))
        {
            ymin=ymin+(ymax-ymin)*1.0/3;
        }
        else
            ymax=ymin+(ymax-ymin)*2.0/3;
    }
    return disy((ymax+ymin)/2,x);
}
int findcenter(POINT &pstart)
{
    xmax=0,xmin=10000;
    int n=pvec.size();    
    for(int i=0;i<n;i++)
    {
        xmax=max(xmax,pvec[i].x);
        xmin=min(xmin,pvec[i].x);
    }
    while((xmax-xmin)/2>accu)
    {
        if(dis(xmin+(xmax-xmin)*1.0/3)<dis(xmin+(xmax-xmin)*2.0/3))
        {
            xmax=xmin+(xmax-xmin)*2.0/3;
        }
        else
            xmin=xmin+(xmax-xmin)*1.0/3;
    }
    pstart.x=(xmax+xmin)/2;
    pstart.y=(ymax+ymin)/2;
    int r=disy((ymax+ymin)/2,(xmax+xmin)/2);
    return r;
}

 

光顾着想用什么新方法比较好了，结果没想到最好的方法。。。