Two Smallest Difference of Ad

刚开始做这道题的时候，只知道是去找二维坐标下的两个点最近的距离，然后除2即为所求。

但是并不知道该怎么求，discuss里面有关于这道题的讨论，就是把所有的点分治。。。。。好吧，其实我不懂

我的思路很简单，就是按照每个点y的坐标排序，从下到上以此这样，如果y坐标相等，则按x坐标从小到大排序，然后二重循环找任意两个点的距离，然后记录最小值

然而。。。在计算最小值的时候总是超时，这里面有一个小技巧，就是算距离的时候是按照勾股定理然后开根号所得的，那么现有的三角形只要有一条直角边大于现有的距离，那么这个三角形的斜边必然大于最小距离。。。so，这算是一个剪枝吧。。。

下面的代码是手写的快排和归并排序，没错，是的，这是一个老师布置的作业

#include <iostream>#include <stdio h="">#include <fstream>#include <algorithm>#include <math h="">using namespace std;typedef struct node{ double x,y;}point;const int inf=(1<<31)-1;const int maxn=100005;point num[maxn],ln[maxn],rn[maxn];int n;double d=inf;bool cmp_y(point a,point b){ if(a.y!=b.y) return a.y<b.y; else return a.x<b.x;} double dis(point a,point b){ return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));} int pivot(point *num,int left,int right){ point a=num[left]; int p=left; for(int j=left+1;j<=right;j++) { if(a.y>num[j].y) { p++; swap(num[p],num[j]); } } swap(num[p],num[left]); return p;} void Merge(point *num,int left,int mid,int right){ int n1=mid-left+1,n2=right-mid,i,j,k; for(i=0;i<n1;i++) ln[i]=num[i+left]; ln[n1].x=ln[i].y=inf; for(i=0;i<n2;i++) rn[i]=num[mid+i+1]; rn[n2].x=rn[n2].y=inf; i=j=0; for(k=left;k<=right;k++) { if(ln[i].y<rn[j].y){ num[k]=ln[i++]; } else num[k]=rn[j++]; } return ;} void Sort(point *num,int left,int right){ if(left<right) {// int mid=pivot(num,left,right); int mid=(left+right)/2;// Sort(num,left,mid-1); Sort(num,left,mid); Sort(num,mid+1,right); Merge(num,left,mid,right); } return ;} int main(){ while(~scanf("%d",&n)) { if(n==0) break; for(int i=0; i<n; i++) scanf("%lf %lf",&num[i].x,&num[i].y);// sort(num,num+n,cmp_y); Sort(num,0,n-1);// for(int i=0;i<n;i++)// printf("%.2f %.2f\n",num[i].x,num[i].y);// int flag=0; d=inf; double x=inf; for(int i=0;i<n;i++){ for(int j=i+1;j<n j="" if="" num="" j="" y-num="" i="" y="">d) break; d=min(d,dis(num[i],num[j])); } } d/=2.0; printf("%.2f\n",d); } return 0;}</n></math></algorithm></fstream></stdio></iostream>