点集最短欧几里德距离

这个是经典的分治算法，大多数书上都有描述。

我却很少找到完整的代码，于是自己写了一个

 

整体是把归并排序的框架做了修改，复杂度是的O(n*logn)，经测试100万个点5秒左右能出结果

不过有个奇怪的地方一直不明白……我用java写了一样的代码，运行速度和c++居然差不多（貌似还略快），100万个点也是5秒左右……

感觉是这个c++的代码还不够完善……

 

#include<vector> #include<cstdio> #include<algorithm> #include<cmath> #define EPI 0.000001 using namespace std; const int maxn=100000; struct point{ double x,y; int ID; point(double _x,double _y,int _ID):x(_x),y(_y),ID(_ID){} point(){} }; typedef const point* PP; inline double Distance(const point &a,const point &b){ return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } inline double min(double a,double b){ return a<b?a:b; } inline bool within_strip(const point &a,double mid,double delta){ return a.x>mid-delta && a.x&l