Sicily 1-D closet pair

写在前面：题意应该是求closest pair，老师大概觉得这题太naïve了

题目：

Description
Given a set S = {p1, p2, …, pn} of n points in the 1 dimensional space, find the two points of S whose distance is the smallest.

Input
There are multiple cases.
For each case, the first line is the number of points n (2<=n<=100000). The second line contains n real numbers.

Output
For each case, output the smallest distance, the answer is rounded to 6 digits after the decimal point.

Sample Input
4
1 3 6 10

Sample Output
2.000000

Hint
This problem is simple. We require that you use the divide-and-conquer algorithm to solve it. Please do not use any sort functions.

//
题意说明：给你一堆一维坐标上的点，让你求最短的两点间距离

思路：
用分治法且不用sort函数

1.只有一个点，返回无穷
2.只有两个点，返回两点间距离
3.其他情况将问题分割：
分为左右两个区间，分别找出左右区间最短的两点间距离
再找到分立于两区间的最短的两点间距离
三者取最小值

简单看来很简单，其实很多细节要注意
比如，如何分割？
如何找分立于两个区间的点的最小距离？

这里分割直接模仿了快速排序的partition，并且用了简单的版本
其他具体内容见代码

tips：要分清物理上（内存里）一排点，和抽象的数轴上的一排点
做这道题很容易（可能只有我）把内存上相邻的点想成数轴上相邻，而忽略了分割需要调整位置的步骤

#include <cstdio>
#include <cstdlib>
#include <algorithm>
using namespace std;

#define INF 1<<30 
//一个大数，后面会用到

double points[100001];//用于存储输入的点的信息

int partition(int begin, int end) {//对一个给定的区间，找一个主元，把数组分为
    int i = begin - 1;                   //一边小于主元，一边大于，就是快速排序里的分割方法
    for (int j = begin; j < end; j++) {//不过这里直接选最后一个为主元，会引起一些问题
        if (points[j] <= points[end]) {//但是后面会处理
            i++;
            swap(points[i], points[j]);
        }
    }
    swap(points[i + 1], points[end]);
    return i+1;
}

double get_closest(int begin,int end) {
    if (begin-end+1 == 1)//当前区间只有一个点，返回无穷，用一个大数代替
        return INF;
    else if (begin - end + 1 == 2)//两个点，直接返回这两个点距离
        return abs(points[begin]-points[end]);
    else {
        int mid = partition(begin, end);
        //分割区间并得到分割的主元，注意这个主元还不确定是左区间的最后一个点
        int left_end = mid == begin ? mid : mid - 1;
        //由于可能刚好选的主元是原来区间的最小最大值，如果把mid当成左区间的末尾或者
        //右区间的起始，都会有一种情况无法分割造成递归不终止
        //所以刚好是最小值时，就把mid当做左边区间的末尾
        //否则就是右边区间的起始，这样至少会分割一个点出去
        int right_begin = left_end + 1;
        //以上三行就把当前区间分成了左右两部分
        //左：[begin,left_end]，右：[right_begin,end]

        double left_max = points[begin], right_min = points[right_begin];
        for (int i = begin + 1; i <= left_end; i++)
            if (points[i] > left_max)
                left_max = points[i];
        for (int i = right_begin + 1; i <= end; i++)
            if (points[i] < right_min)
                right_min = points[i];
        //两个循环找到左区间的最大值和右区间的最小值
        double ans_mid = abs(left_max - right_min);
        //得到两个区间间最近的两个点的距离

        double ans_left = get_closest(begin, left_end);
        double ans_right = get_closest(right_begin, end);
        //递归求得左右区间最短距离

        double ans = ans_left < ans_right ? ans_left : ans_right;
        ans = ans_mid < ans ? ans_mid : ans;
        return ans;
    }
}

int main() {
    int n;
    while (scanf("%d",&n)!=EOF) {
        for (int i = 0; i < n; i++) 
            scanf("%lf", &points[i]);//因为输入还蛮多的，用scanf算了
        printf("%.6lf\n", get_closest(0,n-1));
    }
}