实验1 递归与分治算法.doc

淮海工学院计算机工程学院 实验报告书 课程名： 《算法分析与设计》 题 目： 实验1 递归与分治算法 班 级： 软件102班 学 号： 111003215 姓 名： 鹿 迅 实验1 递归与分治算法 实验目的和要求 （1）进一步掌握递归算法的设计思想以及递归程序的调试技术； （2）理解这样一个观点：分治与递归经常同时应用在算法设计之中。 （3）分别用蛮力法和分治法求解最近对问题； （4）分析算法的时间性能，设计实验程序验证分析结论。 实验内容 设p1=(x1, y1), p2=(x2, y2), …, pn=(xn, yn)是平面上n个点构成的集合S，设计算法找出集合S中距离最近的点对。 实验环境 Turbo C 或VC++ 实验学时 2学时，必做实验 数据结构与算法 #includeiostream #includeCmath #includealgorithm #define N 100 using namespace std; struct point { int x,y; }; bool cmpx(point a,point b) { return a.xb.x; } bool cmpy(point a,point b) { return a.yb.y; } int Sqrt(point a,point b)//两点间的距离的平方 { int k; k=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y); return k; } double min(double d1,double d2){ if(d1=d2) return d1; else return d2; } bool different(point p[],int n,int start,int end) { for(int i=start;iend;i++) { if(p[i].x==p[end].xp[i].y==p[end].y) return true; } return false; } int ClosestPoints1(point p[],int n,point one,point two)//求出两最近点的相关信息 { int d=9999; for(int i=0;in;i++) for(int j=i+1;jn;j++) { int temp=Sqrt(p[i],p[j]); if(tempd) { one=p[i]; two=p[j]; d=temp; } } return d; } int ClosestPoints2(point S[],int n) { if(n2) return 10000; if(n==2) { int d=Sqrt(S[0],S[1]); return d; } sort(S,S+n,cmpx);//按照X大小排序 int m=S[n/2].x; int t=n/2;int i=0; point S1[5000],S2[5000]; for(i=0;it;i++) { S1[i]=S[i]; } for(i=t;in;i++) { S2[i-t]=S[i]; } int d1=ClosestPoints2(S1,t); int d2=ClosestPoints2(S2,n-t); int d=min(d1,d2); point p1[N],p2[N]; int j=0; int p1l=0,p2l=0; //对x坐标差值在2d之间的点进行归类，找到这些点 for(i=0;it;i++) { if(abs(S1[i].x-m)d) { p1[p1l]=S1[i]; p1l++; } } for(i=0;in-t;i++) { if(abs(S2[i].x-m)d) { p2[p2l]=S2[i]; p2l++; } } //对两个区间内的点沿Y坐标轴进行排序 sort(p1,p1+p1l,cmpy); sort(p2,p2+p2l,cmpy); int md=9999; for(i=0;ip1l;i++) {