[ZOJ 1450] Minimal Circle（最小圆覆盖问题）

题目描述：

给出一组点的坐标，求出能够覆盖掉他们的最小圆的坐标及半径.

算法分析：

解决这个问题有一种算法，随机化增量算法，在随机数据下可以O（N）的解决这个问题.

具体实现方法：

1.先将所有的点随机化处理

2.按顺序把点一个一个的加入（一步一步的求前i个点的最小覆盖圆），每加入一个点就进入步骤3

3.判断当前点是否在当前的最下覆盖圆内，如果不在进入4，在进入2

4.将圆心设为当前点 i 半径为 0，然后把前 i-1个点加入进这个最小覆盖圆中，每加入一个点进入5

5.如果发现当前j号点在当前的最小圆的外面，那么说明点j也一定在前j个点（包括i）的最小覆盖圆边界上，我们转到6来再进一步确定这个圆，否则前j个点（包括i）的最小覆盖圆与前i-1个点（包括i）的最小覆盖圆是一样的，则不需要更新，返回

6.此时已经确认点i，j一定在前j个点（包括i）的最小覆盖圆的边界上了，那么我们可以把当前圆的圆心设为第i个点与第j的点连线的中点，半径为到这两个点的距离（就是找一个覆盖这两个点的最小圆），然后重新把前j-1个点加入这个圆中（还是类似上面的步骤，只不过这次我们提前确定了两个点在圆上，目的是求出包含点i，j的前k个点的最小覆盖圆），每加入一个点就进入7

7.如果发现当前k号点在当前的最小圆的外面，那么说明点k也一定在前k个点（包括i，j）的最小覆盖圆边界上，我们不需要再进一步确定这个圆了（因为三个点能确定一个圆！），直接求出这三点共圆，否则前k个点（包括i，j）的最小覆盖圆与前k-1个点（包括i，j）的最小覆盖圆是一样的，则不需要更新。

题目链接：

传送门

代码如下：

#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <cstdlib>
const int maxm=1005;
const double eps=1e-9;
struct Point{
    double x,y;
    int tim;
};
int n;
Point a[maxm],o;
double r;
double distance(Point p1,Point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
void makecir(Point p1,Point p2,Point p3)//三点确定一个圆
{
    double x1=2*(p2.x-p1.x),y1=2*(p2.y-p1.y);
    double c1=p2.x*p2.x+p2.y*p2.y-p1.x*p1.x-p1.y*p1.y;
    double x2=2*(p3.x-p1.x),y2=2*(p3.y-p1.y);
    double c2=p3.x*p3.x+p3.y*p3.y-p1.x*p1.x-p1.y*p1.y;
    o.x=(y1*c2-y2*c1)/(y1*x2-y2*x1);
    o.y=(x2*c1-x1*c2)/(y1*x2-y2*x1);
    r=distance(o,p1);
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        if(!n) break;
        for(int i=1;i<=n;i++)
         scanf("%lf%lf",&a[i].x,&a[i].y);
        std::random_shuffle(a+1,a+n+1);
        o=a[1],r=0;
        for(int i=2;i<=n;i++)
        if(distance(o,a[i])>r+eps)
        {
            o=a[i],r=0;
            for(int j=1;j<=i-1;j++)
            if(distance(o,a[j])>r+eps)
            {
                o.x=(a[i].x+a[j].x)/2;
                o.y=(a[i].y+a[j].y)/2;
                r=distance(o,a[j]);
                for(int k=1;k<=j-1;k++)
                 if(distance(o,a[k])>r+eps)
                  makecir(a[i],a[j],a[k]);
            }
        }   
        printf("%.2lf %.2lf %.2lf\n",o.x,o.y,r);
    }
    return 0;
}