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));
}