最近对问题-递归与分治算法.docVIP

最近对问题-递归与分治算法

《算法分析与设计》实验报告 -PAGE4-

实验1递归与分治算法

一，实验目的和要求

（1）进一步掌握递归算法的设计思想以及递归程序的调试技术；

（2）理解这样一个观点：分治与递归经常同时应用在算法设计之中。

（3）分别用蛮力法和分治法求解最近对问题；

（4）分析算法的时间性能，设计实验程序验证分析结论。

二，实验内容

设p1=(x1,y1),p2=(x2,y2),…,pn=(xn,yn)是平面上n个点构成的集合S，设计算法找出集合S中距离最近的点对。

三，实验环境

TurboC或VC++

四，实验学时

2学时，必做实验

五，数据结构与算法

#includeiostream.h

#includecmath

#defineTRUE1

#defineFALSE0

typedefstructNode

{

doublex;

doubley;

}Node;//坐标

typedefstructList

{

Node*data;//点

intcount;//点的个数

}List;

typedefstructCloseNode

{

Nodea;

Nodeb;//计算距离的两个点

doublespace;//距离平方

}CloseNode;

intn;//点的数目

}

//左右各距中线d的区域的最近对算法

voidmiddle(constListL,CloseNodecnode,intmid,doublemidX)

{

inti,j;//分别表示中线左边，右边的点

doubled=sqrt(cnode.space);

i=mid;

while(i=0L.data[i].x=(midX-d))//在左边的d区域内

{

j=mid;

while(L.data[++j].x=(midX+d)j=L.count)//在右边的d区域内

{

if(L.data[j].y(L.data[i].y-d)||L.data[j].y(L.data[i].y+d))//判断纵坐标是否在左边某固定点的2d区域内

continue;

doublespace=square(L.data[i],L.data[j]);

if(cnode.spacespace)//在满足条件的区域内依次判断

{

cnode.a=L.data[i];

cnode.b=L.data[j];

cnode.space=space;

}

}

--i;

}

}

//分治法求最近对

voidDivideConquer(constListL,CloseNodeclosenode,intbegin,intend)

{

if(begin!=end)

{

intmid=(begin+end)/2;//排列后的中间的那个点

doublemidX=L.data[mid].x;

DivideConquer(L,closenode,begin,mid);//继续在左半边用分治法求最近对

DivideConquer(L,closenode,mid+1,end);//继续在右半边用分治法求最近对

middle(L,closenode,mid,midX);//判断左右各距中线d的区域，是否有最近对

}

}

voidmain()

{

//初始化

Listlist;

CloseNodeclosenode;

closenode.space=10000;//最近点的距离

create(list);//输入各点到NList中

cout各点坐标为:endl;

for(inti=0;ilist.count;++i)

coutX=list.data[i].xY=list.data[i].y\n;

BruteForce(list,closenode,0,list.count-1);

cout用蛮力法求