UVA 10245(分治)

题目链接：https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1186

 

题目大意：求出最近两点间距离

 

题目思路：使用分治算法，先根据x排序，然后取中间点把一个平面分成左右两部分，左右都是n/2个点左右，最短距离要么是左边部分的点之间的最短距离，要么是右边部分的点之间的最短距离，要么是跨区域之间的最短距离。假设我们知道左边部分的点之间的最短距离和右边部分的点之间的最短距离，那么我们只需要求出跨区域的最短距离即可。跨区域只用比x到分界线距离<左右部分最短距离的点就行。求左边部分和右边部分最短距离就直接继续分两部分，只有一个点的时候就返回一个大值，两个点就求出两个点距离，以此类推。

 

以下是代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<math.h>
using namespace std;
#define inf 0x3f3f3f3f
#define MAXN 10005
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)
#define LL long long
#define MOD 998244353
struct node{
    double x,y;
}a[MAXN];
bool cmp(node a,node b){
    return a.x<b.x;
}
double dis(node a,node b){
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double solve(int l,int r){
    if(l==r)return inf;
    int mid=(l+r)>>1;
    int ll=mid,rr=mid;
    double d=min(solve(l,mid),solve(mid+1,r));
    while(ll>=l&&a[mid].x-a[ll].x<d)ll--;
    while(rr<=r&&a[rr].x-a[mid].x<d)rr++;
    rep(i,ll+1,mid){
        rep(j,mid+1,rr-1){
            double temp=dis(a[i],a[j]);
            d=min(d,temp);
        }
    }
    return d;
}
int main()
{
    int n;
    while(~scanf("%d",&n)&&n){
        rep(i,1,n){
            scanf("%lf%lf",&a[i].x,&a[i].y);
        }
        sort(a+1,a+n+1,cmp);
        double ans=solve(1,n);
        if(ans>=10000){
            printf("INFINITY\n");
        }
        else{
            printf("%.4lf\n",ans);
        }
    }
    return 0;
}