算法上机（二）C语言用分治算法解决最近点对问题

最近点对问题

二维空间上有很多个点，每个点的坐标为（x，y），求距离最近的两个点的坐标和距离

解决方法

考虑两种方法。第一种方法，是暴力解法，简单易懂，但时间复杂度高，当点的数量多时，运行困难。第二种方法，利用分治的思想递归地解决最近点对问题。

暴力法

float baoli(Point points[], int length, Point &a, Point &b){
	float min = 99999;
	for(int i=0;i<length;i++){
		for(int j=0;j<i;j++){
			double temp = sqrt(pow(points[i].x-points[j].x,2)+pow(points[i].y-points[j].y,  2));
			if(temp<min){
				min = temp;
				a = points[i];
				b = points[j];
			} 
		}
	}
	return min;
}

暴力法的思想很简单，这里就不多作解释了

分治法

分治法解决最近点对问题的思路如下：
先将点对按横坐标的大小排序，分为左半区和右半区，递归调用函数求出左右两个半区的最近点对和最小距离d1，d2，再求出中间点的min（d1，d2）邻域的最近点对和最小距离，跟min（d1，d2）比较，取最小者。
代码如下。

#include <stdio.h>
#include <math.h> 

//Point 是一个点的结构体，x，y表示横纵坐标 
typedef struct Point{
	float x;	
	float y;
}Point;
//求两点距离的函数 
float distance(Point a,Point b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

//快排 
int partition(Point a[],int begin,int end){
	Point key = a[end];
	int i = begin-1;
	for(int j=begin;j<end;j++){
		if(a[j].x<=key.x){
			i++;
			Point temp = a[i];
			a[i] = a[j];
			a[j] = temp;
		}
	}
	i++;
	a[end] = a[i];
	a[i] = key;
	return i;
}
void quiksort(Point a[],int begin,int end){
	if(begin<end){
		int middle = partition(a,begin,end);
		quiksort(a,begin,middle-1);
		quiksort(a,middle+1,end);
	}
}
//分治法求解距离 
float closest (Point* points,int n,Point &a, Point &b){
	//递归函数出口,两个点直接返回距离，一个点返回一个最大值 
	if(n==2){
		a = points[0];
		b = points[1];
		return distance(a,b);
	} 
	else if(n<2){
		return 99999;
	}
	//分治 
	else{
		int n1 = n/2;
		int n2 = n-n1;
		float d1,d2,dis,mid;
		Point points1 [n1];
		Point points2 [n2];
		Point a1,a2,b1,b2;
		//按x坐标排点对 
		quiksort(points,0,n-1);
		mid = points[n1].x;
		for(int i=0;i<n1;i++){
			points1[i] = points[i];
		}
		for(int i=0;i<n2;i++){
			points2[i] = points[i+n1];
		}
		//左右两个分区 
		d1 = closest(points1,n1,a1,b1);
		d2 = closest(points2,n2,a2,b2);
		if(d1<=d2){
			a = a1;
			b = b1;
			dis = distance(a,b);
		}
		else{
			a = a2;
			b = b2;
			dis = distance(a,b);
		}
		//中间区域
		Point points3[n];
		int k = -1;
		for(int i=0;i<n;i++){
			if(abs(points[i].x-mid)<=dis){
				points3[++k] = points[i];
			}
		}
		for(int i=0;i<n;i++){
			for(int j=i+1;j<i+7&&j<k;j++){
				if(distance(points3[i],points3[j])<dis){
					a = points3[i];
					b = points3[j];
					dis = distance(a,b);
				}
			}
		} 
		return dis;
		
	}
}
int main(void){
	Point points[7] = {{1.0,2.0},{1.1,6.2},{5.0,7.0},{15.0,12.5},{2.6,2.5},{5.9,4.6},{9.1,15.8}};
	for(int i=0;i<7;i++){
		printf("(%f,%f)\n",points[i].x,points[i].y);
	}

	Point a,b;
	float dis = closest(points,7,a,b);
	printf("距离%f，点a（%f，%f），点b（%f，%f）",dis,a.x,a.y,b.x,b.y);
}