分治法求解最近点对问题

1.问题

分治法解最近对问题。

2.解析

1.当n<=3，通过蛮力算法求解；
2.当n>3：对所有的点按照x坐标进行排序，根据x进行点集划分，求中位线，将点集以中位线为界划分，在中轴线处做一条垂线，将点集分为n/2和n/2两个子集Pl和Pr，不断对左右两个部分分别求最近点对的距离d，然后比较选出最短距离dmin =min(d1,d2)；选出到中位线的距离均小于min的点，对集合Pl和Pr中的每个点 ，在y坐标区间[y,y+d]内最对取出8个候选点，计算与点p的最近距离d3；计算该部分点集中最小距离d3，dmin=min(dmin ,d3)。

3.设计

int m = (L + H)/2;//分治区间
double d1 = FF(p,L,m), d2 = FF(p,m+1,H), d = min(d1,d2);
while(p[l].x<p[m].x-d && l<=H){
    ++l;
}
while(p[r].x>p[m].x+d && l<=L){
    --r;
}

sort(p+l,p+r+1,CY); //筛选后的所有点以y坐标排序
for(int i=l;i<r;i++){
    for(int j=i+1;j<r;j++){
        if(p[j].y-p[i].y > d){
            break;
        }
        else{
            int d3 = Dis(p[i],p[j]);
            if(d3 < d){
                d = d3;
            } 
        }
    }
}

}

// 按x坐标再次排序
sort(p,p+n,CX);
double d = FF(p,1,n);
cout << d  ;
return 0;

}

double Dis(P p1,P p2){
//计算两点之间的距离函数
return sqrt((p1.x - p2.x)(p1.x - p2.x) + (p1.y - p2.y)(p1.y - p2.y));
}

4.分析

时间复杂度复杂度算上排序：O(NlogN)

5.源码

#include
#include
#include
#include
#include
#include
#define max 10

using namespace std;

struct P{
int x;
int y;
};//P点集，x升序排列
double Dis(P p1,P p2);
bool CY(P p1,P p2){
return p1.y < p2.y;
}
bool CX(P p1,P p2){
return p1.x < p2.x;
}
double FF(P p[],int L,int H)
{ int l=L,r=H;
double dmin;
if(H - L == 2){
int d1 = Dis(p[L],p[L+1]);
int d2 = Dis(p[L],p[H]);
int d3 = Dis(p[L+1],p[H]);
d1=min(d1,d2);
d3=min(d1,d3);
return d3;
}
if(H - L == 1){
return Dis(p[L],p[H]);//调用两点函数
}
int m = (L + H)/2;//分治区间
double d1 = FF(p,L,m), d2 = FF(p,m+1,H), d = min(d1,d2);
while(p[l].x<p[m].x-d && l<=H){
++l;
}
while(p[r].x>p[m].x+d && l<=L){
–r;
}

sort(p+l,p+r+1,CY); //筛选后的所有点以y坐标排序
for(int i=l;i<r;i++){
    for(int j=i+1;j<r;j++){
        if(p[j].y-p[i].y > d){
            break;
        }
        else{
            int d3 = Dis(p[i],p[j]);
            if(d3 < d){
                d = d3;
            } 
        }
    }
}

}
int main(){
P p[max];
int n;
cin >> n;
for(int i=1;i<=n;i++){
cin >> p[i].x >> p[i].y;
}
// 按x坐标再次排序
sort(p,p+n,CX);
double d = FF(p,1,n);
cout << d ;
return 0;
}

double Dis(P p1,P p2){
//计算两点之间的距离函数
return sqrt((p1.x - p2.x)(p1.x - p2.x) + (p1.y - p2.y)(p1.y - p2.y));
}