分治策略:
#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;
}