应用示例:
假设一个金属片上钻n个孔,如果孔距离的太近就会发生断裂,需要判断最近距离就是最近对问题
double Closest(point S[], int low, int high); //实现求最近对距离
double Distance(point a, point b); //求两点之间距离
int Partition(point r[], int first, int end); //划分
void QuickSort(point r[], int first, int end); //快速排序-分治法
#include <iostream>
#include<iomanip>
#include<cmath>
using namespace std;
struct point{
int x,y;
};
int Partition(point r[], int first, int end){ //划分
int i = first, j = end; //初始化待划分区间
while (i < j) {
while (i < j && r[i].y <= r[j].y) j--; //右侧扫描
if (i < j) {
point temp = r[i]; r[i] = r[j]; r[j] = temp; //将较小记录交换到前面
i++;
}
while (i < j && r[i].y <= r[j].y) i++; //左侧扫描
if (i < j) {
point temp = r[i]; r[i] = r[j]; r[j] = temp; //将较大记录交换到后面
j--;
}
}
return i; // 返回轴值记录的位置
}
double Distance(point a,point b){
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
void QuickSort(point r[], int first, int end){ //快速排序
int pivot;
if (first < end) {
pivot = Partition(r, first, end); //划分,pivot是轴值在序列中的位置
QuickSort(r, first, pivot - 1); //求解子问题1,对左侧子序列进行快速排序
QuickSort(r, pivot + 1, end); //求解子问题2,对右侧子序列进行快速排序
}
}
double Closest(point S[],int low,int high)
{
double d1,d2,d3,d;
int mid,i,j,index;
point P[8];
if(high-low==1)
return Distance(S[low],S[high]);
if(high-low==2){
d1=Distance(S[low],S[low+11]);
d2=Distance(S[low+1],S[high]);
d3=Distance(S[low],S[high]);
if((d1<d2)&&(d2<d3))
return d1;
else if(d2<d3)
return d2;
else
return d3;
}
mid=(low+high)/2;
d1=Closest(S,low,mid);
d2=Closest(S,low+1,high);
if(d1<=d2)
d=d1;
else
d=d2;
index=0;
for(i=mid;(i>=low)&&(S[mid].x-S[i].x<d);i--)
{
P[index++]=S[i];
}
for(i=mid+1;(i<=high)&&(S[i].x-S[mid].x<d);i++)
{
P[index++]=S[i];
}
QuickSort(P,0,index-1);
for(i=0;i<index;i++)
{
for(j=i+1;j<index;j++)
{
if(P[j].y-P[i].y>=d)
break;
else
{
d3=Distance(P[i],P[j]);
if(d3<d)
d=d3;
}
}
}
return d;
}
void main(){
point S[8] = { {10,2},{15,2},{13,74},{42,22},{65,8},{70,12},{60,54},{62,2} };
double minDist = Closest(S, 0, 8-1);
//输出最近点对的距离
cout << "最近点对距离为:"<< minDist << endl;
}
运行结果:
本人原创,如需转载需告知;