题目大意是找一个最小的圆覆盖给你一堆点中的两个点,输出改圆的半径
基础最近点对模板题
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <map>
#define LL __int64
using namespace std;
const int maxn = 100010;
struct point
{
double x,y;
};
point p1[maxn],p2[maxn];
int n;
bool cmpx(point a,point b)//按x升序
{
return a.x < b.x;
}
bool cmpy(point a,point b)//按y降序
{
return a.y < b.y;
}
double dis(point x,point y)
{
double v = (x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y);
return sqrt(v);
}
double solve(int l,int r)
{
if(l+1==r) return dis(p1[l],p1[r]);//这里有点像搜索,如果到了最近几个数,我就直接比较输出
if(l+2==r) return min(dis(p1[l],p1[l+1]),min(dis(p1[l+1],p1[l+2]),dis(p1[l],p1[r])));
int mid = (l+r)>>1;//二分思想
double ans = min(solve(l,mid),solve(mid+1,r));
int cnt = 0;
for(int i=l;i<=r;i++){ //在这个区域找点
if(fabs(p1[i].x-p1[mid].x)<=ans)
p2[cnt++] = p1[i];
}
sort(p2,p2+cnt,cmpy);
for(int i=0;i<cnt;i++){
for(int j=i+1;j<cnt;j++){
if(p2[j].y - p2[i].y >= ans) break;//优化判断
ans = min(ans,dis(p2[i],p2[j]));
}
}
return ans;
}
int main()
{
while(scanf("%d",&n)&&n){
for(int i=0;i<n;i++) scanf("%lf%lf",&p1[i].x,&p1[i].y);
sort(p1,p1+n,cmpx);
double ans = solve(0,n-1);//找出直径
printf("%.2f\n",ans/2.0);//输出半径
}
return 0;
}