用分治算法解平面最接近点对问题.docVIP

一. 用分治算法解平面最接近点对问题 1.题目 关于最接近点对问题: 给定平面上n个点，找出其中一对点，使得在n个点所构成的所有点对中，该点对的距离最小。 2. 程序详细介绍(各模块的功能等) 本程序主要包括两个类:类Point和类Ppoint.其中类Point为处理一些的基本数据传递等.类Ppoint为该程序的主要实现模块,该类中有输入点对的函数shuru,对所输入的点对按X轴排序的函数sort,求各点对的距离的函数xiao等. 假设S中的点为平面上的点，它们都有2个坐标值x和y。为了将平面上点集S线性分割为大小大致相等的2个子集S1和S2，我们选取一垂直线lx=m）来作为分割直线。其中m为S中各点x坐标的中位数。由此将S分割为S1={p∈S|px≤m}和S2={p∈S|pxm}。从而使S1和S2分别位于直线l的左侧和右侧，且S=S1∪S2 。由于m是S中各点x坐标值的中位数，因此S1和S2中的点数大致相等。递归地在S1和S2上解最接近点对问题，我们分别得到S1和S2中的最小距离δ1和δ2 4. 调试与测试：调试方法，测试结果（包括输入数据和输出结果）的分析与讨论 运行该程序时,屏幕上会出现一个界面,首先该界面会提示输入要处理的点对个数,输入点对个数后从键盘输入数字0即可显示出处理后的各个结果,会出现如下结果: 5.程序代码(源程序) #includeiostream #includecmath #includefstream using namespace std; int i,j,k,d,m,n; double p2,q,s,r,t; class Point //创建一个点类// { public: double x; double y; double getx() { return x; } double gety() { return y; } friend class Ppoint; }; class Ppoint { int sum; double juli[10][10]; double min[11]; //min[10]用来存放每组中最短的距离// double mini[11]; //mini[10]用来存放每组中距离最短的点对中的第一个点// double minj[11]; //minj[10]用来存放每组中距离最短的点对中的第二个点// Point p[100]; Point p1; public: void shuru() { cout请输入要处理的点的个数endl; cinsum; for(i=0;isum;i++) { cout请输入点对 endl; cinp[i].x; cinp[i].y; coutp[i].x,p[i].yendl; } } void sort() { cout以下是按x轴上由小到大排序后的点对endl; for(i=0;isum-1;i++) { for(j=i+1;jsum;j++) { if(p[i].xp[j].x) { p1=p[i]; p[i]=p[j]; p[j]=p1; } } } for(i=0;isum;i++) { coutp[i].x,p[i].yendl; } } void xiao() { cout以下是对每个模块中的点求距离endl; for(k=0;ksum/10;k++) { cout按任意键继续endl; cini; cout以下是第k个模块中的距离endl; for(i=1;i10;i++) { for(j=0;ji;j++) { r=abs(p[k*10+i].x-p[k*10+j].x); t=abs(p[k*10+i].y-p[k*10+j].y); juli[i][j]=sqrt(r*r+t*t); coutjuli[i][j]\t; } coutendl; } min[k]=juli[k][0],mini[k]=10*k+1,minj[k]=10*k; for(i=1;i10;i++)