bzoj 1336/bzoj 1337: [Balkan2002]Alien最小圆覆盖 随机增量法

题意

给出N个点，让你画一个最小的包含所有点的圆。输出圆心和半径。
n<=100000。

分析

随机增量法裸题。

假设已经求出了前i-1个点构成的最小圆，且第i个点不在该圆内，那么显然我们取前i-1个点中的若干个与点i构成一个圆，则点i一定在边界上。
那么我们枚举1到i-1中的点j，假设已经求出了1到j-1中的点和i构成的最小圆，且点j不在该圆内，那么点j一定在1到j中的点和i构成的最小圆的边界上。
那么现在已经确定了两个边界上的点，考虑确定第三个点。枚举1到j-1中的点k，若点k不在1到k-1中的点与i,j两点构成的最小圆内，那么点k一定在1到k中的点与i,j两点构成的最小圆内。这样就确定下了边界上的三个点i,j,k，从而可以得到覆盖1到k中的点加上i,j两点的最小圆。

复杂度看似是n^3的，但其实随机情况下的期望复杂度是O(n)的。
证明：(from http://blog.csdn.NET/lthyxy/article/details/6661250)

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;

const int N=100005;

int n;
double r;
struct data{double x,y;}a[N],o;

double dis(data a,data b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

data get_o(data p1,data p2,data p3)
{
    double a=p1.x-p2.x,b=p1.y-p2.y,c=-((p1.x+p2.x)/2*a+(p1.y+p2.y)/2*b);
    double d=p1.x-p3.x,e=p1.y-p3.y,f=-((p1.x+p3.x)/2*d+(p1.y+p3.y)/2*e);
    data o;
    o.x=(b*f-c*e)/(a*e-b*d);o.y=(f*a-c*d)/(b*d-e*a);
    return o;
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%lf%lf",&a[i].x,&a[i].y);
    for (int i=1;i<=n;i++) swap(a[rand()%n+1],a[rand()%n+1]);
    o=a[1];
    for (int i=2;i<=n;i++)
    {
        if (dis(a[i],o)<=r) continue;
        o.x=(a[i].x+a[1].x)/2;o.y=(a[i].y+a[1].y)/2;r=dis(a[i],a[1])/2;
        for (int j=2;j<i;j++)
        {
            if (dis(a[j],o)<=r) continue;
            o.x=(a[i].x+a[j].x)/2;o.y=(a[i].y+a[j].y)/2;r=dis(a[i],a[j])/2;
            for (int k=1;k<j;k++)
            {
                if (dis(a[k],o)<=r) continue;
                o=get_o(a[i],a[j],a[k]);
                r=dis(a[i],o);
            }
        }
    }
    printf("%.10lf\n%.10lf %.10lf",r,o.x,o.y);
    return 0;
}