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>