图论论文最小生成树算法城市高速公路问题中的应用

1、xxxx研究生堂下考试答卷2012-2013学年第 1 学期2012年 12 月 18 日最小生成树在城市高速公路问题中的应用摘 要： 城市高速公路问题就是以最短高速路程连接一组城市的问题，在城市规划和建设中应用广泛。本文以最小生成树在城市高速公路问题中的应用为例，利用最小生成树的三种算法的分析和研究，阐明了最小生成树在最优化方面的作用。关键词：城市高速公路问题 prim算法 kruskal算法 简易算法一 引言图论是数学的一个分支。它以图为研究对象。在图论的课程体系中，图结构是一种非常重要的非线性数据结构。带权图的最小生成树尤其被广泛应用在解决工程技术及科学管理等各个领域的最优化问题中。二

2、背景知识1 图和树：图论中的图是由若干给定的点及连接两点的线所构成的图形，这种图形通常用来描述某些事物之间的某种特定关系，用点代表事物，用连接两点的线表示相应两个事物间具有这种关系。树是五圈连通无向图，如果树t的节点数为n，那么树的边数为n-1。2 生成树：连通图 g 上的一个子图，该子图连通，无回路且包含图g 的所有节点，称为连通图的极小连通子图。一个连通图可以有多棵不同的生成树。3 最小生成树：对一个带权连通图，也有多可不同的生成树。由于该图是带权图，各边的权值不一定相等，因此这些生成树的各边权值之和也不一定相同，其中权值最小的生成树被称为该带权连通图的最小生成树。4 高速公路问题：假设有

3、 n 个城市，第 i 个城市的位置笛卡尔坐标为 (xi，yi)，每条公路可以连接两个城市。目前原有的公路有m 条，但是不能实现所有城市之间的连通，因此需要继续修建公路，在费用最低的原则下，实现 n 个城市的连通，还需要修建哪些条公路。由于修路的费用与公路的长短是成正比的，所以这个问题就可以转化成求修建哪几条公路能够实现所有城市的连通，同时满足所修公路总长最短。三 最小生成树的求解方法构造最小生成树可以有多种算法。大多数图论教材中介绍了其中的两种算法 prim 算法和 kruskal 算法，本文另介绍一种简易算法来实现最小生成树的构造。1 prim 算法思想：普里姆算法通过逐个往生成树上添加顶点

4、来构造连通网的最小生成树。算法具体步骤：（1） 开始：选取连通网中的任意一个顶点添加到最小生成树中。（2） 重复执行以下操作：1） 连通网的顶点集合分成两个部分：已经添加到最小生成树中的顶点集合和尚未添加到最小生成树中的顶点集合；2） 找出所有连通这两个集合中顶点的边；3） 从中选取一条权值最小的边添加到生成树中，同时将与这条边相连的顶点也添加到生成树中。（3）结束：所有的顶点都被添加到最小生成树中。2 kruskal 算法思想：通过逐个往生成树上添加边来构造连通网的最小生成树。算法具体步骤：（1） 将连通网中的所有顶点添加到最小生成树中，构造一个森林；（2） 将各边按照权值从小到大排序；（3

5、） 按照排好的顺序向生成树中添加不使森林中产生回路的边 (若构成回路则不添加，继续考察下一条边)。直至该森林变成一棵树为止。3 简易算法思想：通过逐个从连通网中删除边来构造最小生成树。算法具体步骤：（1） 将连通网中各边按照权值从大到小排序；（2） 按照排好的顺序从连通网中删除权值最大的边，条件是使删除该边后的子图仍然保持连通(若删除后子图不连通则改边保留，继续删除下一条边)。直至子图中任何一条边都不能删除 （即删除任意一条边都会造成该子图不连通）为止。4 三种算法的比较（1） 普里姆算法：主要是对顶点进行操作；采用邻接矩阵作为存储结构，在行过程中对连通网中的每一个顶点都考察到了，因此普里姆算

6、法的时间复杂度为（n 为连通网中顶点的个数）。普里姆算法适用于求边稠密的连通网的最小生成树。（2） 克鲁斯卡尔算法：主要是对边进行操作，时间复杂度主要取决于对边按照权值进行排序的时间，边的个数为e，排序的时间复杂度可以做到o (eloge)，因此算法的时间复杂度为 o(eloge)。克鲁斯卡尔算法适用于求边稀疏的连通网的最小生成树。（3） 简易算法：主要是对边进行操作，时间复杂度主要取决于对边按照权值进行排序的时间，边的个数为e，排序的时间复杂度可以做到 o (eloge)，因此算法的时间复杂度为 o(eloge)。该算法适用于求边稀疏的连通网的最小生成树。四 应用利用最小生成树来解决高速公路

7、问题，将高速公路问题中的城市看做图中的顶点，城市之间修建的道路看做图中顶点之间的边，城市之间所修修建的公路的长度看做是图中个边上的权值。这样我们就把高速公路问题转换成了求一个有向连通网的最小生成树问题。此时假设城市个数为6，分别为a，b，c，d，e，f。并设其对应城市之间的公路距离权值及初始状态的连通无向图如下所示：边(a,b）(a,c）(a,d）(b,c)(b,e) (c,d)(c,e)(c,f)(d,f)(e,f)权值61553564261 简易算法来求解最小生成树（1） 实现步骤：从权值最大的边开始进行删除 (e,f),(c,e) ,(a,b),(a,d)被删除后都没有破坏连通性，所以这

8、些边可以从图中删除，得下图：当删除第五条边(b,c)时，造成了图的连通性的破坏，所以该条边不能被删除必须保留。下图为最终构造好的最小生成树：（2） 算法实现：1）连通网的存储结构表示城市高速公路网络的无向连通网采用邻接矩阵的存储方式进行存储。由于需要区分已经存在的公路和需要修建的公路，所以在每条边上增加一个标志位，同时为了给所有的边排序，因此单独建立一个表示边信息的结构体数组结构。具体实现如下：typedef char towntype; typedef float roadtype;typedef struct int n; /*图的顶点个数 */int m; /* 图的边个数 */town

9、type towns maxvex ; /* 顶点信息 */roadtype roadsmaxvexmaxvex ;/*边信息 */ graphmatrix;typedef struct int start_town, stop_town; /* 边的起点和终点 */roadtype weight; /* 边的权 */enum exist,unexist ex;/*区别已经建好的公路和未修建的公路 */ edge;edge mst 50 ;2） 判断图的连通性函数判断无向图的的连通性有很多方法，这里采用的是通过对图进行深度优先搜索，统计遍历过的顶点个数，如果顶点个数比图中顶点个数少，说明该图不

10、连通，相反说明该图是连通图。具体实现如下：void dfsmatrix (graphmatrix gm, int i, int n , int &visited ) int j;printf “(%d”,i) ;visited i =1;for (int j=0;jn;j+)if ( gm i j ! =0 & gm i j ! =maxvalue& ! visited j)dfsmatrix (gm,j,n) ;遍历结束后可以通过统计 visited 数组中的值为 1 的元素的个数 v 来确定访问过的结点个数，如果 vn 说明该无向图不连通，反之说明该无向图连通。graphthrough (

11、graphmatrix gm)for (i=0;in;i+) visited i =0dfsmatrix (gm, 0, n , visited) ;for (i=0;in;i+)if ( visited i =1) num+;if (num=n) return 1;else return 0;3） 对权值进行由大到小排序由于整个算法的时间复杂度主要取决于排序算法的效率高低，因此在这里我们采用的是快速排序算法qsort (edge&mst, int i, int j )，排序的时间复杂度可以做到o (eloge)。具体算法实现就不再赘述。4） 简易算法的主体部分int easy (graphm

12、atrix &gm , edge mst )int i, j, num = 0, th,weight;for (i=0;igm.m;i+)if (mst i .ex=unexist)weight=gm. roads start_town stop_town ;gm. roads start_town stop_town =0;th=graphthrough (gm) ;if (th=1) continue;else gm. roads start_town stop_town =weight;2 prim算法求解最小生成树（1）实现步骤：从节点a开始，依次添加节点c，f，d，b，e，并依次添加

13、对应的边，各个步骤如下图所示：（2） 算法实现：#include#include#include#define infinity 1000#define max_name 50#define max_vertex_num 50typedef char vertexmax_name;/顶点名字串typedef int adjmatrixmax_vertex_nummax_vertex_num;/邻接距阵typedef struct vertex adjvex; /邻接矩阵 int lowcost; /权值closemax_vertex_num;/定义一个结构以便在后面closedge 使用typ

14、edef struct/定义图 vertex vexsmax_vertex_num; /顶点集 adjmatrix arcs; /边 int vexnum,arcnum;/点个数，边个数mgraph;int locatevex(mgraph g,vertex u)/若g中存在顶点u,则返回该点在图中位置;否则返回其他信息; int i; for(i=0;ig.vexnum;+i) if(strcmp(u,g.vexsi)=0) return i; return 1;void creategraph(mgraph &g) int i,j,k,w; vertex v1,v2; printf(输入无

15、向图顶点数和边数: n); scanf(%d %d,&g.vexnum,&g.arcnum); printf(输入各顶点的值：n, g.vexnum); for(i=0;ig.vexnum;+i) /构造顶点集 scanf(%s,&g.vexsi); for(i=0;ig.vexnum;+i) /初始化邻接方阵 for(j=0;jg.vexnum;+j) g.arcsij=infinity; printf(输入一条边依附的顶点及权值： n,g.arcnum);/输入一条边依附的顶点及权值 for(k=0;kg.arcnum;+k) scanf(%s%s%d,v1,v2,&w); i=locat

16、evex(g,v1);/v1在图中位置 j=locatevex(g,v2);/v2在图中位置 g.arcsij=g.arcsji=w; /置于对称弧 int minimum(close c,mgraph g)/求出下一个节点 第k个顶点 int i=0,j,k,min; min=infinity; /初始化 k=-1; for(j=0;j=g.vexnum;j+)/求最小 if(cj.lowcost0) min=cj.lowcost; k=j; return k;void prim(mgraph g,vertex u) int i,j,k=0; close closedge;/一个结构 boo

17、l isbreak=false; k=locatevex(g,u);/u在图中位置 返回g.vexsi中下标 for(j=0;j=g.vexnum;+j) /辅助数组初始化 closedge从o 开始 if(j!=k)/没有自己到自己的 closedgek.lowcost=0; strcpy(closedgej.adjvex,u); closedgej.lowcost=g.arcskj;/列 int flag1000; flag0=0; int count=1; for(i=1;ig.vexnum;+i) k=minimum(closedge,g); if(k=-1)isbreak=true;

18、break;printf(%s-%s%dn,closedgek.adjvex,g.vexsk,g.arcsklocatevex(g,closedgek.adjvex); /输出生成树的边 closedgek.lowcost=0; / 第k个顶点并入u集flagcount=k;count+; for(j=0;jg.vexnum;+j) if(g.arcskjclosedgej.lowcost)/新顶点并入u后重新选择最小边 strcpy(closedgej.adjvex,g.vexsk); closedgej.lowcost=g.arcskj; if(isbreak) printf(此图不连通，

19、无最小支撑树 ！n访问过的点为：n);for(i=0;icount;i+) printf(%s ,g.vexsflagi);printf(n); void main() mgraph g; creategraph(g); printf(最小生成树的各条边为: n); prim(g,g.vexs0); 3 kruskal算法求解最小生成树（1）实现步骤：依次添加边（a,c）, （d,f）, （b,e）, （c,f）, （b,c），并依次添加对应节点，各个步骤结果如下图：（2） 算法实现：#include#include#define m 20#define max 20typedef struc

20、t /构造边 int begin; int end; int weight; /权值edge;typedef struct int adj; int weight;adjmatrixmaxmax;typedef struct adjmatrix arc; int vexnum, arcnum; /顶点数和边数mgraph;void creatgraph(mgraph *); /函数申明 构造图void sort(edge* ,mgraph *); /函数申明 对边的排序void minispantree(mgraph *); /最小生成树int find(int *, int ); void

21、swapn(edge *, int, int); /交换两条边的权值和它们的起点和终点 void creatgraph(mgraph *g) /构件图g int i, j,n, m; printf(请输入图的顶点数和边数:); scanf(%d %d,&g-vexnum,&g-arcnum); for (i = 1; i vexnum; i+) /初始化图 for ( j = 1; j vexnum; j+) g-arcij.adj = g-arcji.adj = 0; for ( i = 1; i arcnum; i+) /输入边和权值 printf(n请输入边的起始点和终点:); scan

22、f(%d %d,&n,&m); while(n g-vexnum | m g-vexnum) printf( n);printf( n); printf(输入的数字不符合要求 请重新输入:); scanf(%d%d,&n,&m); g-arcnm.adj = g-arcmn.adj = 1; getchar(); printf(n请输入%d与%d之间的权值: , n, m); scanf(%d,&g-arcnm.weight); /输入权值 void sort(edge edges,mgraph *g) /对权值进行排序 int i, j; for ( i = 1; i arcnum; i+)

23、 for ( j = i + 1; j arcnum; j+) if (edgesi.weight edgesj.weight) swapn(edges, i, j); printf( n); printf( n); printf( n); printf( 权 从 小 到 大 排 序 之 后 为:n); printf(n); printf( 起点 终点 权n); for (i = 1; i arcnum; i+) printf( %dn, edgesi.begin, edgesi.end, edgesi.weight); void swapn(edge *edges,int i, int j)

24、 /交换权值 以及头和尾 int temp; temp = edgesi.begin; edgesi.begin = edgesj.begin; edgesj.begin = temp; temp = edgesi.end; edgesi.end = edgesj.end; edgesj.end = temp; temp = edgesi.weight; edgesi.weight = edgesj.weight; edgesj.weight = temp; void minispantree(mgraph *g)/生成最小生成树 int i, j, n, m; int k = 1; int parentm; edge edgesm; for ( i = 1; i vexnum; i+) for (j = i + 1; j vexnum; j+) if (g-arcij.adj = 1) edgesk.begin = i; edgesk.end = j; edgesk.weight = g-arcij.weight; k+; sort(edges, g); for (i = 1; i arcnum; i+) parenti = 0; printf(n);printf( n);printf( n);printf( n); printf( 最 小 生 成