题意
给出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;
}