HDU1007：Quoit Design——题解

http://acm.hdu.edu.cn/showproblem.php?pid=1007

题目大意：给n个点，求点对最短距离/2.

——————————————————————

平面分治裸题。

暂时还不想讲为什么这么做。

所以原理暂割。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef double dl;
const dl INF=1e20;
const int N=;
struct node{
    dl x;
    dl y;
}p[N],a[N],b[N],c[N];
bool cmp(node A,node B){
    return A.x<B.x;
}
inline dl dis(int i,int j){
    return sqrt(pow(b[i].x-c[j].x,)+pow(b[i].y-c[j].y,));
}
dl solve(int l,int r){
    if(l>=r)return INF;
    int mid=(l+r)>>;
    dl x0=(p[mid].x+p[mid+].x)/2.0;
    dl d=min(solve(l,mid),solve(mid+,r));
    int l1=l,r1=mid+,bnum=,cnum=;
    for(int i=l;i<=r;i++){
    if(l1<=mid&&(r1>r||p[l1].y<p[r1].y)){
        a[i]=p[l1++];
        if(x0-d<a[i].x)b[++bnum]=a[i];
    }else{
        a[i]=p[r1++];
        if(a[i].x<x0+d)c[++cnum]=a[i];
    }
    }
    for(int i=l;i<=r;i++)p[i]=a[i];

    for(int i=,j=;i<=bnum||j<=cnum;){
    if(i<=bnum&&(j>cnum||b[i].y<c[j].y)){
        for(int k=j-;k>=;k--){
        if(b[i].y-d>=c[k].y)break;
        d=min(d,dis(i,k));
        }
        i++;
    }else{
        for(int k=i-;k>=;k--){
        if(c[j].y-d>=b[k].y)break;
        d=min(d,dis(k,j));
        }
        j++;
    }
    }
    return d;
}
int main(){
    int n;
    while(scanf("%d",&n)!=EOF&&n){
    for(int i=;i<=n;i++){
        scanf("%lf%lf",&p[i].x,&p[i].y);
    }
    sort(p+,p+n+,cmp);
    printf("%.2f\n",solve(,n)/2.0);
    }
    return ;
}