最小圆覆盖,随机增量法.

#include<iostream>
#include<cstdio>
#include<cmath>
/*
    algorithm:
　　A、令Ci表示为前i个点的最小覆盖圆。当加入新点pi时如果pi不在Ci-1里那么pi必定在Ci的边界上。
　　B、再从新考虑这样一个问题，Ci为前i个点最小覆盖圆且p在Ci的的边界上！同理加入新点pi时如果p
　　i不在Ci-1里那么pi必定在Ci的边界上。这时我们就包含了两个点在这个最小圆的边界上。
　　C、再从新考虑这样一个问题，Ci为前i个点最小覆盖圆且有两个确定点再边界上！此时先让
　　O(N)的方法能够判定出最小圆。
 　　------------------------------------------------------------------------------------
　　analysis:
　　现在来分析为什么是线性的。
　　C是线性的这是显然的。
　　B<-C的过程中。考虑pi 他在园内的概率为 (i-1)/i 。在圆外的概率为 1/i 所以加入pi的期望复杂度为：
    (1-i)/i*O(1) +(1/i)*O(i) {前者在园内那么不进入C，只用了O(1)。
    后者进入C用了O(i)的时间}这样分析出来，复杂度实际上仍旧
　　是线性的。
　　A<-B的过程中。考虑方法相同，这样A<-B仍旧是线性。于是难以置信的最小圆覆盖的复杂度变成了线性的。
*/ 
using namespace std;
struct node
{
    double x,y;
};

int n;
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)  //给出两条直线ax+by+c=0,dx+ey+f=0 求交点
{                                  //注意到三角形里两条中垂线不可能平行，所以不会产生除0错误
      O.y=(c*d-f*a)/(b*d-e*a);
     O.x=(c*e-f*b)/(a*e-b*d);
}

int main()
{
     freopen("circle.in","r",stdin);
     freopen("circle.out","w",stdout);
     int m;
    scanf("%d",&m);
    scanf("%d",&n);
    int l,ll;
    for (int i=1;i<=n;++i) 
    {
      scanf("%lf%lf",&p[i].x,&p[i].y);
    }
     O=p[1];r=0;                  //初始C1
     
     for(int i=2;i<=n;++i)            //A
     if(dist(O,p[i])>r+1e-6) 
     {
         O=p[i];r=0;
         for (int j=1;j<=i-1;++j)        //B
         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)        //C
             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]);
             } 
         }
     }
     printf("%.2lf\n",r);
     //while(1);
     return 0;
 }