最小圆覆盖

强悍的随机化线性算法。

这类算法首先考虑暴力应该如何达到暴力最优复杂度（一般来说是有部分操作最快可以o（1）进行），然后用一定递归思想分割问题，并考虑在随机情况下最坏情况出现概率，从而找出最坏情况出现并不影响整体。

（注意先随机）对于最小圆覆盖来说，假设ci为前i个点的解，若i不在ci-1中，则i一定在ci的圆周上，递归至前i-1，已确定i在圆周上求解，同样可确定第二个点，再递归至第三层则可直接确定（3点确定一个圆，当然也有两点情况）

这种思想可拓至线性规划，但是暂时没有研究，单单是二维就觉得维护最高点很麻烦（或许可从暴力找灵感，因为其本质与暴力很像）。

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
struct point
{
  double x,y,z;
}a[200],o;
int n,m,st[4],t;
double r;
const double eps=1e-10;
double dist(point e,point u) {double x=e.x-u.x,y=e.y-u.y;return sqrt(x*x+y*y);}
void cross(point e,point u,point &v){v.x=e.y*u.z-e.z*u.y,v.y=e.z*u.x-e.x*u.z,v.z=e.x*u.y-e.y*u.x;}
void make3()
{
  point e=a[st[1]],u=a[st[2]],v=a[st[3]],p,q,l1,l2;
  cross(e,u,l1),cross(e,v,l2);
  p.x=(e.x+u.x)/2,p.y=(e.y+u.y)/2,p.z=1;
  q.x=p.x+l1.x,q.y=p.y+l1.y,q.z=1;
  cross(p,q,l1);
  p.x=(e.x+v.x)/2,p.y=(e.y+v.y)/2,p.z=1;
  q.x=p.x+l2.x,q.y=p.y+l2.y,q.z=1;  
  cross(p,q,l2);
  cross(l1,l2,o);
  o.x/=o.z,o.y/=o.z,o.z=1,r=dist(o,e);
}
void make2()
{
  point e=a[st[1]],u=a[st[2]];
  o.x=(e.x+u.x)/2,o.y=(e.y+u.y)/2,r=dist(e,u)/2;
}
int check(int i)
{
  double k=dist(a[i],o);
  if (r-k<-eps) return 1;else return 0;
}
void doit3(int k)
{
  for (int i=1;i<=k;i++) if (check(i)) st[t=3]=i,make3();
}
void doit2(int k)
{
  st[t=2]=1,make2();
  for (int i=2;i<=k;i++) if (check(i)) st[t=2]=i,make2(),doit3(i-1);
}
void init()
{
  int i,l,p,q;
  point e;
  memset(a,0,sizeof(a)),memset(st,0,sizeof(st));
  for (i=1;i<=n;i++) scanf("%lf%lf\n",&a[i].x,&a[i].y),a[i].z=1;
  for (i=1,l=n*10;i<=l;i++) p=1+rand()%n,q=1+rand()%n,e=a[p],a[p]=a[q],a[q]=e;
  st[1]=1,st[2]=2,t=2,make2();
  for (i=3;i<=n;i++) if (check(i)) st[t=1]=i,doit2(i-1);
  printf("%.3lf\n",r);
}
int main()
{
  srand((int)time(NULL));
    scanf("%d\n",&n);
    init();
  return 0;
}