【算法】最近对问题（分治法）

分治策略：

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<cmath>
using namespace std;

struct Point{
	int x,y;
};

double Distance(Point a,Point b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int Partition(Point a[],int first,int end) // 划分 
{
	int i=first,j=end; // 初始化待划分区间 
	while(i<j){
		while(i<j && a[i].y<=a[j].y) j--; // 右侧扫描 
		if(i<j){
			Point temp = a[i];
			a[i]=a[j];
			a[j]=temp;
			i++;
		}
		while(i<j && a[i].y<=a[j].y) i++; // 左侧扫描 
		if(i<j){
			int tx = a[i].x,ty = a[i].y;
			a[i].x=a[j].x; a[i].y=a[j].y;
			a[j].x=tx; a[j].y=ty;
			j--;
		}		
	}	
	return i;
} 

void QuickSort(Point a[],int first,int end)
{
	int pivot; // 枢轴在序列中的位置 
	if(first<end){
		pivot = Partition(a,first,end); // 划分
		QuickSort(a,first,pivot-1); // 对左侧子序列进行快速排序
		QuickSort(a,pivot+1,end);  // 对右侧子序列进行快速排序
	} 
}

double Closet(Point S[],int low,int high)
{
	double d1,d2,d3,d;
	Point P[100]; // 存放点的集合 
	if(high-low==1){ // 只有两个点，返回两个点的距离 
		return Distance(S[low],S[high]);
	}
	if(high-low==2){ // 只有三个点，求最近对距离 
		d1 = Distance(S[low],S[low+1]);
		d2 = Distance(S[low+1],S[high]);
		d3 = Distance(S[low],S[high]);
		if(d1<d2 && d2<d3) 
			return d1;
		else if(d2<d3)
			return d2;
		else 
			return d3;
	}
	// 多个点，分治
	int mid = (low+high)/2; // 计算中间点
	d1 = Closet(S,low,mid); // 递归求解子问题,左边 
	d2 = Closet(S,mid+1,high);  // 递归求解子问题，右边 
	if(d1<=d2) // 取两边的最短距离 为 d 
		d = d1; 
	else 
		d = d2;
	// 求跨两个部分的最短距离d3
	int index = 0;
	for(int i=mid;(i>=low)&&(S[mid].x-S[i].x<d);i--) // 建立点集合p1 
		P[index++] = S[i]; 
	for(int i=mid+1;(i<=high)&&(S[i].x-S[mid].x<d);i++) //建立点集合p2 
		P[index++] = S[i]; 
	QuickSort(P,0,index-1); // 对集合p1,p2按y坐标升序排列（快速排序）
	for(int i=0;i<index;i++){
		for(int j=i+1;j<index;j++){
			if(P[j].y-P[i].y>=d) // 超出y坐标范围，点P[i]处理完毕 
				break;
			else{
				d3 = Distance(P[i],P[j]);
				if(d3<d) d = d3;
			}
		}
	} 
	return d;
}

int main()
{
	int n;
	cin>>n;
	Point s[100];
	for(int i=0;i<n;i++){
		cin>>s[i].x;
		cin>>s[i].y;
	}
	double result = Closet(s,0,n-1);
    cout<<result<<endl;
	return 0;
}