Quoit Design---hdu1007（最近点对问题 分治法）

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <bitset>
#include <iostream>
#include <time.h>
#include <vector>
#include <queue>

typedef long long LL;

using namespace std;

const int N = 1e6+1;
const double eps = 1e-10;
const int INF = 0x3f3f3f3f;
const int mod = 1000000007;
const double PI = 4*atan(1.0);

struct point
{
    double x, y;
}p[N], t[N];

double dist(point p1, point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y));
}

bool cmpX(point p1, point p2)
{
    return p1.x < p2.x;
}
bool cmpY(point p1, point p2)
{
    return p1.y < p2.y;
}

double Find(int L, int R)
{
    if(L+1 == R)///当只有两个点的时候；
        return dist(p[L], p[R]);
    if(L+2 == R)///三个点的时候；
        return min(min(dist(p[L], p[L+1]), dist(p[L], p[R])), dist(p[L+1], p[R]));

    int Mid = (L+R)/2;

    double Min_d  = min(Find(L, Mid), Find(Mid+1, R));///找到两边的最小值， 下面更新中间部分的；

    int cnt = 0;

    for(int i=L; i<=R; i++)
    {
        if(fabs(p[i].x-p[Mid].x) <= Min_d)
            t[cnt++] = p[i];///把可能是最近点对中的点加入t集合中去；
    }

    sort(t, t+cnt, cmpY);///排序，下面找到，t集合中最近的两点间的距离，跟新Min_d;

    for(int i=0; i<cnt; i++)
    {
        for(int j=i+1; j<cnt; j++)
        {
            if( t[j].y - t[i].y > Min_d )
                break;
            Min_d = min(Min_d, dist(t[i], t[j]));
        }
    }
    return Min_d;
}

int main()
{
    int n;

    while(scanf("%d", &n), n)
    {
        for(int i=0; i<n; i++)
            scanf("%lf %lf", &p[i].x, &p[i].y);

        sort(p, p+n, cmpX);///按横坐标x排序；

        double ans = Find(0, n-1);///递归 求解；

        printf("%.2f\n", ans/2);
    }
    return 0;
}