贪心算法设计实验报告.docx

1、贪心算法设计实验报告贪心算法设计实验报告 算法设计技巧与分析课程实验报告 五、实验步骤 1.用Kruskal算法实现最小生成树 算法描述：假设 WN=(V,E) 是一个含有 n 个顶点的连通网，则按照克鲁斯卡尔算法构造最小生成树的过程为：先构造一个只含 n 个顶点，而边集为空的子图，若将该子图中各个顶点看成是各棵树上的根结点，则它是一个含有 n 棵树的一个森林。之后，从网的边集 E 中选取一条权值最小的边，若该条边的两个顶点分属不同的树，则将其加入子图，也就是说，将这两个顶点分别所在的两棵树合成一棵树；反之，若该条边的两个顶点已落在同一棵树上，则不可取，而应该取下一条权值最小的边再试之。依次类

2、推，直至森林中只有一棵树，也即子图中含有 n-1条边为止。 下面给出c语言代码实现及说明 本程序对树的存储主要是以边为存储对象，即边的结构体里面有这样几个参数：1，边的权值。2，边的一个顶点。3，边的另一个顶点。4，边是否属于生成树的一条边（即最小生成树边标志）。由该程序的存储结构决定了该算法比较适用于边稀疏的情形，至于边稠密的情况会在下面的Prim算法中给出。 在Kruskal算法中有两个比较重要的部分1，对边按权重排序。2，对一条边加入子树后是否会产生回路的判断即判断边的两个节点是否在同一个树中（集合里） 对于问题1：可以有各种排序算法，读者可以自行选择自己喜欢的排序算法来替换代码中的排序

3、算法。（本处使用选择排序算法（效率较低），读者可以自己修改为快速排序或者是对排序） 下面主要讲解解决问题2，解决这个问题一般借用不相交集的思想，即在本程序中每次以同集合中所有点的编号最小的数来标识本集合。（对于根节点即标号最小的节点则标记为0）例如对于上图实例，下面给出详细的不相交集的维护过程（给出前几步的详细说明） a,初始状态（*分别在数组的第1，2，3，4，5，6，7）0号位空置（均为0） b,选择第一条边A-D(将D的标记改为1)（因为AD最短） 对表进行维护（维护后仍同上表，因为还没有两个集合合并） C，选择第二条边C-E（修改上表） 对上表进行维护（任同上表，因为还没有两个集合合并

4、） D，选择第三条边（D-F）(根据条件DF两点不再同一集合，改边可选) 然后就合并DF两点所在的集合D的前去是1，即A标记为0，E的标记也为0，合并因为61所以表修改如下 以后几步均如上判断两点是否在一个集合从而判断改边是否可取，并维护上表 下面附上源代码 /* Kruskal算法的实现 09网一 殷赛 *-*13 输入：图G(用结构体数组来存储每条边，包含每条边的节点） 输出：图G的最小生成树树 */ #includestdio.h #includestdlib.h typedef struct Edge char dot_1; char dot_2; int weight; int le

5、ap; Edge; Edge* selectionsort(Edge *array,int n)/选择排序（对边按权重由高到低排序） int i,j,min,temp; for(i=0;ii+) min=i; for(j=i+1;jj+) if(arraymin.weightarrayj.weight) min=j; if(min!=i) temp=arrayi.weight; Edge *Kruskal(Edge *Graph,int num_e,int *V,int num_v)/克鲁斯卡尔算法实现 if(Vm!=Vnm!=Vnn!=Vm)/如果边的两个顶点不再一个集合则边是生成树的边(注

6、意首节点的标记和集合里非首节点的标记不同) arraymin.weight=temp; temp=arrayi.dot_1; arrayi.dot_1=arraymin.dot_1; arraymin.dot_1=temp; temp=arrayi.dot_2; arrayi.dot_2=arraymin.dot_2; arraymin.dot_2=temp; return array; int m,n,test; int i,j,t,k; for(i=0;inum_e;i+) for(j=1;jnum_v+1;j+) if(Graphi.dot_1=Vj) m=j; if(Graphi.do

7、t_2=Vj) n=j; Graphi.leap=1; if(Vn=0) else /维护不相交集 k=1;/对每个节点都检查是否为标记合格节点 while(knum_v+1) if(Vk!=0)/只要标记不为0都进行整理,只不过有些节点的标记整理前后是 if(Vm=0) else if(VmVn) else VVn=Vm; VVm=Vn; if(mVn) else Vm=Vn; VVn=m; if(nVm) else Vn=Vm; VVm=n; 一样的(即标记符合标准的节点) if(Vm=0Vn=0)/如果边的两个顶点是两个集合的首节点则可以合并 Graphi.leap=1; if(mn)

8、else Vn=m; Vm=n; k+; t=Vk; while(Vt!=0) Vk=t; t=Vt; /维护不相交集 k=1; while(knum_v+1) if(Vk!=0) t=Vk; while(Vt!=0) Vk=t; t=Vt; k+; /* printf(不相交集的情况：n); for(test=1;testnum_v+1;test+) printf(%-4c,Vtest); printf(n); for(test=1;testnum_v+1;test+) printf(%-4d,Vtest); printf(n);*/ return Graph; void main() in

9、t i,j,num_v,num_e,cost=0; Edge *Graph=NULL; int *V=NULL; printf(请输入土中有多少个顶点!n); scanf(%d,num_v); V=(int*)malloc(sizeof(int*)*2); for(i=0;ii+) Vi=(int*)malloc(sizeof(int)*(num_v+1); for(i=0;ii+) for(j=0;jnum_v+1;j+) Vij=0; for(i=1;inum_v+1;i+) printf(请输入第%d个顶点:,i); scanf( %c,Vi); printf(请输入图中有多少条边!n)

10、; scanf(%d,num_e); Graph=(Edge*)malloc(sizeof(Edge)*num_e); for(i=0;inum_e;i+) Graph=selectionsort(Graph,num_e); printf(请输入第%d条边的权值和两个顶点！n,i+1); scanf(%d %c %c,Graphi.weight,Graphi.dot_1,Graphi.dot_2); Graphi.leap=0; /以上部分是存储图 /- Graph=Kruskal(Graph,num_e,V,num_v); printf(构成最小生成树的边和顶点分别是：n); printf(

11、顶点1-顶点2-边n); for(i=0;inum_e;i+) if(Graphi.leap=1) %c printf( %c - - %dn,Graphi.dot_1,Graphi.dot_2,Graphi.weight); cost=cost+Graphi.weight; printf(最小生成树的权值是：%dn,cost); 2.用Prim算法实现最小生成树 算法描述：假设V是图中顶点的集合,E是图中边的集合,TE为最小生成树中的边的集合,则prim算法通过以下步骤可以得到最小生成树: 1:初始化:U=u 0,TE=f。此步骤设立一个只有结点u 0的结点集U和一个空的边集TE作为最小生成

12、树的初始形态,在随后的算法执行中,这个形态会不断的发生变化,直到得到最小生成树为止。 2:在所有uU,vVU的边(u,v)E中,找一条权最小的边(u 0,v 0)，将此边加进集合TE中，并将此边的非U中顶点加入U中。此步骤的功能是在边集E中找一条边,要求这条边满足以下条件:首先边的两个顶点要分别在顶点集合U和VU中,其次边的权要最小。找到这条边以后,把这条边放到边集TE中,并把这条边上不在U中的那个顶点加入到U中。这一步骤在算法中应执行多次,每执行一次,集合TE和U都将发生变化,分别增加一条边和一个顶点,因此,TE和U是两个动态的集合,这一点在理解算法时要密切注意。 3:如果U=V,则算法结束

13、;否则重复步骤2。可以把本步骤看成循环终止条件。我们可以算出当U=V时,步骤2共执行了n1次(设n为图中顶点的数目),TE中也增加了n1条边,这n1条边就是需要求出的最小生成树的边 （图任如上图） 下面给出具体c语言代码和详解 该代码的存储结构是邻接矩阵上三角用来存储最小生成树的边，树边用-1标记，下三角用来存储原图，方便输出最小生成树上图的最后邻接矩阵如下：（具体实现和解释见代码）(INF代表无穷，即不相邻) /* 09网一 殷赛 *-*13 Prim算法的实现 输入：图G 输出：图G的最小生成树 */ #includestdio.h #includestdlib.h #define INF

14、 * /* Prim算法最重要的部分在于从可以挑选的边中间挑选出最小值 并且对加入后会产生回路的边标记为不可选（此代码中将该边权值标记为INF即无穷大） 每选中一条边就会加入一个点（就必须对与该点连接的边进行维护，即多产生回路边标记为无穷） 找最小边则是从上三角中所有可选的行和列中选择 选中的树边则标记为-1（其权值从下三角读出） */ void Prim(int *Graph,int num_v)/转化为对上三角的维护 int i,j,leap=0,temp; int m,n,min; int *a=(int *)malloc(sizeof(int)*num_v); for(i=0;inum

15、_v;i+) ai=0; a=1; while(leap!=num_v-1) min=INF; for(i=0;inum_v;i+)/搜索上三角中的最小值 if(ai=1) for(j=i+1;jnum_v;j+)/行中搜索 if(minGraphijGraphij0) min=Graphij; m=j; n=i; temp=0; for(j=0;jj+)/列中搜索 if(minGraphjiGraphji0) min=Graphji; /* for(i=0;inum_v;i+)/检查矩阵中在那些行（列）可以选择最小值（即那些列和行是候选行和列） printf(%-4d,ai); m=j; n

16、=i; temp=1; if(temp=0)/为了区别出是在行还是在列中搜索到的元素 else Graphmn=-1; Graphnm=-1; for(i=0;ii+)/去掉行中回路边 for(i=m+1;inum_v;i+)/去掉列中回路边 am=1; if(ai=1Graphmi0) Graphmi=INF; if(ai=1Graphim0) Graphim=INF; printf(n);*/ /* for(i=0;inum_v;i+)/检验对上三角的维护，和下三角是否修改了（下三角保存了原树和用于输出最小生成树） void main() Graph=(int*)malloc(sizeof

17、(int*)*num_v);/动态生成二维数组用来存储图（邻接矩阵） for(i=0;inum_v;i+) for(i=0;inum_v;i+)/初始化矩阵 for(j=0;jnum_v;j+) */ leap+; for(j=0;jnum_v;j+) printf(%-4dt,Graphij); printf(n); int i,j; int num_v,num_e; int *Graph=NULL; char *V=NULL; char ch_1,ch_2; int weight; int m,n; printf(请输入图的顶点数：); scanf(%d,num_v); V=(char*)

18、malloc(sizeof(char)*num_v); Graphi=(int*)malloc(sizeof(int)*num_v); Graphij=INF; for(i=0;inum_v;i+) Graphii=0; for(i=0;inum_v;i+) printf(请输入第%d个顶点:,i+1); scanf( %c,Vi); printf(请输入图的边数：); scanf(%d,num_e); for(i=0;inum_e;i+) printf(请输入第%d条边的顶点和权值:,i+1); scanf( %c %c%d,ch_1,ch_2,weight); for(j=0;jnum_v

19、;j+) if(Vj=ch_1) m=j; if(Vj=ch_2) n=j; Graphmn=weight; Graphnm=weight; /以上是对图用邻接矩阵存储 /- Prim(Graph,num_v); printf(最小生成树如下：n); printf(顶点-顶点-权值n); weight=0; for(i=0;inum_v;i+) for(j=i+1;jnum_v;j+) if(Graphij=-1) printf( %-2c - %-2c - %-2dn,Vi,Vj,Graphji); weight=weight+Graphji; printf(最小生成树的权重是：%dn,weight); 六、实验结果与分析 Kruskal算法适用于边稀疏的情形，而Prim算法适用于边稠密的情形