1、8.1 图的基本概念,第8章 图,8.2 图的存储结构,8.3 图的遍历,8.4 生成树和最小生成树,8.5 最短路径,8.6 拓扑排序,本章小结,8.7 AOE网与关键路径,8.1 图的基本概念 8.1.1 图的定义 图(Graph)G由两个集合V(Vertex)和E(Edge)组成,记为G=(V,E),其中V是顶点的有限集合,记为V(G),E是连接V中两个不同顶点(顶点对)的边的有限集合,记为E(G)。,在图G中,如果代表边的顶点对是无序的,则称G为无向图,无向图中代表边的无序顶点对通常用圆括号括起来,用以表示一条无向边。 如：无向边 ( V1 , V2) 如果表示边的顶点对是有序的,则称
2、G为有向图,在有向图中代表边的顶点对通常用尖括号括起来 。 如：有向边 ,8.1.2 图的基本术语 端点和邻接点在一个无向图中,若存在一条边(vi,vj),则称vi和vj为此边的两个端点,并称它们互为邻接点。在一个有向图中,若存在一条边,则称此边是顶点vi的一条出边,同时也是顶点vj的一条入边；称vi和vj分别为此边的起始端点(简称为起点)和终止端点(简称终点)；称vi邻接到vj 。,2. 顶点的度、入度和出度在无向图中,顶点所具有的边的数目称为该顶点的度。在有向图中,以顶点vi为终点的入边的数目,称为该顶点的入度。以顶点vi为始点的出边的数目,称为该顶点的出度。一个顶点的入度与出度的和为该顶
3、点的度。若一个图中有n个顶点和e条边,每个顶点的度为di(1in),则有：,3. 完全图若无向图中的每两个顶点之间都存在着一条边,有向图中的每两个顶点之间都存在着方向相反的两条边,则称此图为完全图。显然,完全无向图包含有n(n-1)/2条边,完全有向图包含有n(n-1)条边。例如,图(a)所示的图是一个具有4个顶点的完全无向图,共有6条边。图(b)所示的图是一个具有4个顶点的完全有向图,共有12条边。,4. 稠密图、稀疏图当一个图接近完全图时,则称为稠密图。相反,当一个图含有较少的边数(即当en(n-1)时,则称为稀疏图。 5. 子图设有两个图G=(V,E)和G=(V,E),若V是V的子集,即
4、VV,且E是E的子集,即EE,则称G是G的子图。例如图(b)是图(a)的子图,而图(c)不是图(a)的子图。,6. 路径和路径长度在一个图G=(V,E)中,从顶点vi到顶点vj的一条路径是一个顶点序列(vi,vi1,vi2,vim,vj),若此图G是无向图,则边(vi,vi1),(vi1,vi2),(vim-1,vim),(vim,vj)属于E(G)；若此图是有向图,则,属于E(G)。路径长度是指一条路径上经过的边的数目。若一条路径上除开始点和结束点可以相同外,其余顶点均不相同,则称此路径为简单路径。例如,有图中,(v0,v2,v1)就是一条简单路径,其长度为2。,7. 回路或环若一条路径上的
5、开始点与结束点为同一个顶点,则此路径被称为回路或环。开始点与结束点相同的简单路径被称为简单回路或简单环。例如,右图中,(v0,v2,v1,v0)就是一条简单回路,其长度为3。,9. 连通、连通图和连通分量在无向图G中,若从顶点vi到顶点vj有路径,则称vi和vj是连通的。 若图G中任意两个顶点都连通,则称G为连通图,否则称为非连通图。无向图G中的极大连通子图称为G的连通分量。显然,任何连通图的连通分量只有一个,即本身,而非连通图有多个连通分量。,9. 强连通图和强连通分量在有向图G中,若从顶点vi到顶点vj有路径,则称从vi到vj是连通的。若图G中的任意两个顶点vi和vj都连通,即从vi到vj
6、和从vj到vi都存在路径,则称图G是强连通图。例如,右边两个图都是强连通图。有向图G中的极大强连通子图称为G的强连通分量。显然,强连通图只有一个强连通分量,即本身,非强连通图有多个强连通分量。,10. 权和网图中每一条边都可以附有一个对应的数值,这种与边相关的数值称为权。权可以表示从一个顶点到另一个顶点的距离或花费的代价。边上带有权的图称为带权图,也称作网。,例 有n个顶点的有向强连通图最多需要多少条边？最少需要多 少条边？,答：有n个顶点的有向强连通图最多有n(n-1)条边(构成一个有向完全图的情况)；最少有n条边(n个顶点依次首尾相接构成一个环的情况)。,8.2 图的存储结构,8.2.1
7、邻接矩阵存储方法邻接矩阵是表示顶点之间相邻关系的矩阵。设G=(V,E)是具有n(n0)个顶点的图,顶点的顺序依次为(v0,v1,vn-1),则G的邻接矩阵A是n阶方阵,其定义如下：(1)如果G是无向图,则：(2)如果G是有向图,则：,(3)如果G是带权无向图,则：(4)如果G是带权有向图,则：,邻接矩阵的特点如下：(1)图的邻接矩阵表示是惟一的。(2)无向图的邻接矩阵一定是一个对称矩阵。因此,按照压缩存储的思想,在具体存放邻接矩阵时只需存放上(或下)三角形阵的元素即可。(3)不带权的有向图的邻接矩阵一般来说是一个稀疏矩阵,因此,当图的顶点较多时,可以采用三元组表的方法存储邻接矩阵。(4)对于无
8、向图,邻接矩阵的第i行(或第i列)非零元素(或非元素)的个数正好是第i个顶点vi的度。,(5)对于有向图,邻接矩阵的第i行(或第i列)非零元素(或非元素)的个数正好是第i个顶点vi的出度(或入度)。(6)用邻接矩阵方法存储图,很容易确定图中任意两个顶点之间是否有边相连。但是,要确定图中有多少条边,则必须按行、按列对每个元素进行检测,所花费的时间代价很大。这是用邻接矩阵存储图的局限性。,邻接矩阵的数据类型定义如下： #define MAXV typedef struct int no; /*顶点编号*/InfoType info; /*顶点其他信息*/ VertexType; /*顶点类型*/
9、typedef struct /*图的定义*/ int edgesMAXVMAXV; /*邻接矩阵*/int n, e; /*顶点数,弧数*/VertexType vexsMAXV; /*存放顶点信息*/ MGraph;,MGraph g;g:,edgesNN是邻接矩阵,vexsN是顶点的数组,8.2.2 邻接表存储方法图的邻接表存储方法是一种顺序分配与链式分配相结合的存储方法。在邻接表中,对图中每个顶点建立一个单链表,第i个单链表中的结点表示依附于顶点vi的边(对有向图是以顶点vi为尾的弧)。每个单链表上附设一个表头结点。表结点和表头结点的结构如下：表结点（边结点） 表头结点,邻接表的特点如
10、下：(1)邻接表表示不惟一。这是因为在每个顶点对应的单链表中,各边结点的链接次序可以是任意的,取决于建立邻接表的算法以及边的输入次序。(2)对于有n个顶点和e条边的无向图,其邻接表有n个顶点结点和2e个边结点。显然,在总的边数小于n(n-1)/2的情况下,邻接表比邻接矩阵要节省空间。(3)对于无向图,邻接表的顶点vi对应的第i个链表的边结点数目正好是顶点vi的度。(4)对于有向图,邻接表的顶点vi对应的第i个链表的边结点数目仅仅是vi的出度。其入度为邻接表中所有adjvex域值为i的边结点数目。,邻接表存储结构的定义如下： typedef struct ANode /*弧的结点结构类型*/ i
11、nt adjvex; /*该弧的终点位置*/struct ANode *nextarc; /*指向下一条弧的指针*/InfoType info; /*该弧的相关信息*/ ArcNode; typedef struct Vnode /*邻接表头结点的类型*/ Vertex data; /*顶点信息*/ArcNode *firstarc; /*指向第一条弧*/ VNode; typedef VNode AdjListMAXV; /*AdjList是邻接表类型*/ typedef struct AdjList adjlist; /*邻接表*/int n,e; /*图中顶点数n和边数e*/ ALGra
12、ph; /*图的类型*/,例 给定一个具有n个结点的无向图的邻接矩阵和邻接表。(1) 设计一个将邻接矩阵转换为邻接表的算法；(2) 设计一个将邻接表转换为邻接矩阵的算法；(3) 分析上述两个算法的时间复杂度。解：(1)在邻接矩阵上查找值不为0的元素,找到这样的元素后创建一个表结点并在邻接表对应的单链表中采用前插法插入该结点。算法如下：,void MatToList( MGraph g, ALGraph * ,(2)在邻接表上查找相邻结点,找到后修改相应邻接矩阵元素的值。算法如下： void ListToMat(ALGraph *G, MGraph ,(3)上述两个算法的时间复杂度均为O(n2)
13、。对于(2)的算法,若不计算给aij赋初值0的双重for循环语句,其时间复杂度为O(n+e),其中e为图的边数。,8.3 图的遍历 8.3.1 图的遍历的概念从给定图中任意指定的顶点(称为初始点)出发,按照某种搜索方法沿着图的边访问图中的所有顶点,使每个顶点仅被访问一次,这个过程称为图的遍历。如果给定图是连通的无向图或者是强连通的有向图,则遍历过程一次就能完成,并可按访问的先后顺序得到由该图所有顶点组成的一个序列。根据搜索方法的不同,图的遍历方法有两种：一种叫做深度优先搜索法(DFS)；另一种叫做广度优先搜索法(BFS)。,8.3.2 深度优先搜索遍历深度优先搜索遍历的过程是：从图中某个初始顶
14、点v出发,首先访问初始顶点v,然后选择一个与顶点v相邻且没被访问过的顶点w为初始顶点,再从w出发进行深度优先搜索,直到图中与当前顶点v邻接的所有顶点都被访问过为止。显然,这个遍历过程是个递归过程。以邻接表为存储结构的深度优先搜索遍历算法如下(其中,v是初始顶点编号,visited 是一个全局数组,初始时所有元素均为0表示所有顶点尚未访问过)：,void DFS(ALGraph *G,int v) ArcNode *p;visitedv=1; /*置已访问标记*/printf(“%d “,v); /*输出被访问顶点的编号*/p=G-adjlistv.firstarc; /*p指向顶点v的第一条弧
15、的弧头结点*/while (p!=NULL) if (visitedp-adjvex=0) DFS(G,p-adjvex); /*若p-adjvex顶点未访问,递归访问它*/p=p-nextarc; /*p指向顶点v的下一条弧的弧头结点*/ ,例如,以上图的邻接表为例调用DFS( )函数,假设初始顶点编号v=2,给出调用DFS( )的执行过程。,8.3.3 广度优先搜索遍历广度优先搜索遍历的过程是：首先访问初始点vi,接着访问vi的所有未被访问过的邻接点vi1,vi2,vit,然后再按照vi1,vi2,vit的次序,访问每一个顶点的所有未被访问过的邻接点,依次类推,直到图中所有和初始点vi有路
16、径相通的顶点都被访问过为止。以邻接表为存储结构,用广度优先搜索遍历图时,需要使用一个队列,以类似于按层次遍历二叉树遍历图。对应的算法如下(其中,v是初始顶点编号)：,void BFS(ALGraph *G, int v) ArcNode *p; int w,i;int queueMAXV,front=0,rear=0; /*定义循环队列*/int visitedMAXV; /*定义存放结点的访问标志的数组*/for (i=0;in;i+) visitedi=0; /*访问标志数组初始化*/printf(“%2d“,v); /*输出被访问顶点的编号*/visitedv=1; /*置已访问标记*/
17、rear=(rear+1)%MAXV;queuerear=v; /*v进队*/,while (front!=rear) /*若队列不空时循环*/ front=(front+1)%MAXV;w=queuefront; /*出队并赋给w*/p=G-adjlistw.firstarc; /*找w的第一个的邻接点*/while (p!=NULL) if (visitedp-adjvex=0) printf(“%2d”,p-adjvex); /*访问之*/visitedp-adjvex=1; rear=(rear+1)%MAXV;/*该顶点进队*/queuerear=p-adjvex;p=p-nexta
18、rc; /*找下一个邻接顶点*/printf(“n“); ,例如,以上图的邻接表为例调用BFS( )函数,假设初始顶点编号v=2,给出调用BFS( )的执行过程。,8.3.4 非连通图的遍历对于无向图来说,若无向图是连通图,则一次遍历能够访问到图中的所有顶点；但若无向图是非连通图,则只能访问到初始点所在连通分量中的所有顶点,其他连通分量中的顶点是不可能访问到的。为此需要从其他每个连通分量中选择初始点,分别进行遍历,才能够访问到图中的所有顶点；,对于有向图来说,若从初始点到图中的每个顶点都有路径,则能够访问到图中的所有顶点；否则不能访问到所有顶点,为此同样需要再选初始点,继续进行遍历,直到图中的
19、所有顶点都被访问过为止。,采用深度优先搜索遍历非连通无向图的算法如下： DFS1(ALGraph *G) int i;for (i=0;in;i+)if (visitedi=0) DFS(G,i); ,采用广度优先搜索遍历非连通无向图的算法如下： BFS1(ALGraph *G) int i;for (i=0;in;i+)if (visitedi=0) BFS(G,i); ,8.3.5 图遍历算法的应用 例8.3 假设图G采用邻接表存储，设计一个算法，判断无向图G是否连通。若连通则返回1；否则返回0. int Connect(ALGraph *G) int i, flag=1;for (i=0
20、; in; i+) visitedi=0;DFS(G, 0);for (i=0; in; i+)if (visitedi=0) flag=0;break; return flag ; ,例8.4 假设图G采用邻接表存储，设计一个算法，输出图G中从顶点u到顶点v的长度为l 的所有简单路径。 void PathAll(ALGraph *G, int u, int v, int l, int path , int d) / path0d是到当前为止已经走过的路径长度。调用时d初值为-1 int m, i;ArcNode *p;visitedu=1; d+;pathd=u; / 将当前顶点添加到路径中
21、if (u=v ,p=G-adjlistu.firstarc; while ( p!=NULL ) m=p-adjvex;if (visitedm=0)PathAll(G, m, v, l , path, d );p-p-nextarc; visitedu=0; / 恢复环境，使得该顶点可重新使用 ,主函数： viod main( ) int pathMAXV, u, v, l, i, j;MGraph g;ALGraph *G;int AMAXVMAXV= 0, 1, 0, 1, 0 ; 1, 0, 1, 0, 0 ; 0, 1, 0, 1, 1 ; 1, 0, 1, 0, 1 ; 0, 0
22、, 1, 1, 0 ; ;g.n=5; g.e=6;,for (i=0; ig.n; i+)for (j=0; jg.n; j+)g.edgesij=Aij; MatToList(g, G); for (i=0; ig.n; i+) visitedi=0; printf(“图G:”); DispAdj(G); / 输出邻接表 u=1; v=4; l=3; PathAll(G, u, v, l, path, -1); / 输出所有从u到v长度为l 的路径 ,8.4 生成树和最小生成树 8.4.1 生成树的概念一个连通图的生成树是一个极小连通子图,它含有图中全部顶点,但只有构成一棵树的(n-1)条
23、边。,如果在一棵生成树上添加一条边,必定构成一个环：因为这条边使得它依附的那两个顶点之间有了第二条路径。一棵有n个顶点的生成树(连通无回路图)有且仅有(n-1)条边,如果一个图有n个顶点和小于(n-1)条边,则是非连通图。如果它多于(n-1)条边,则一定有回路。但是,有(n-1)条边的图不一定都是生成树。,对于一个带权(假定每条边上的权均为大于零的实数)连通无向图G中的不同生成树,其每棵树的所有边上的权值之和也可能不同；图的所有生成树中具有边上的权值之和最小的树称为图的最小生成树。按照生成树的定义,n个顶点的连通图的生成树有n个顶点、n-1条边。因此,构造最小生成树的准则有三条：(1) 必须只
24、使用该图中的边来构造最小生成树；(2) 必须使用且仅使用n-1条边来连接图中的n个顶点；(3) 不能使用产生回路的边。,8.4.2 无向图的连通分量和生成树在对无向图进行遍历时,对于连通图,仅需调用遍历过程(DFS或BFS)一次,从图中任一顶点出发,便可以遍历图中的各个顶点。对非连通图,则需多次调用遍历过程,每次调用得到的顶点集连同相关的边就构成图的一个连通分量。设G=(V,E)为连通图,则从图中任一顶点出发遍历图时,必定将E(G)分成两个集合T和B,其中T是遍历图过程中走过的边的集合,B是剩余的边的集合：TB=,TB=E(G)。显然,G=(V,T)是G的极小连通子图,即G是G的一棵生成树。,
25、由深度优先遍历得到的生成树称为深度优先生成树；由广度优先遍历得到的生成树称为广度优先生成树。这样的生成树是由遍历时访问过的n个顶点和遍历时经历的n-1条边组成。对于非连通图,每个连通分量中的顶点集和遍历时走过的边一起构成一棵生成树,各个连通分量的生成树组成非连通图的生成森林。,8.4.3 普里姆算法普里姆(Prim)算法是求最小生成树的一种构造性算法。假设G=(V,E)是一个具有n个顶点的带权连通无向图,T=(U,TE)是G的最小生成树,其中U是T的顶点集,TE是T的边集,则由G构造最小生成树T的步骤如下：(1)初始化=v0。v0到其他顶点的所有边为候选边；(2)重复以下步骤n-1次,使得其他
26、n-1个顶点被加入到U中： 从候选边中挑选权值最小的边输出,设该边在V-U中的顶点是v,将v加入U中； 考察当前V-U中的所有顶点vi,修改候选边：若(v,vi)的权值小于原来和vi关联的候选边,则用(v,vi)取代后者作为候选边。,普里姆算法求解最小生成树的过程,将带权连通图用带权邻接矩阵costnn存储：,数组closest0n-1： closesti=j 不是生成树的顶点i到生成树U中最接近的顶点是j:,数组lowcost0n-1存放不在生成树的顶点到生成树U中最接近的顶点的边的权值： Lowcosti=0 若顶点i已经是生成树的顶点 0Lowcosti 若顶点i不是生成树的顶点，此时，
27、lowcosti=(i, closesti )边的权值lowcost =0, 3, 0, 5, 6, 4 如: lowcost5=4 点5到树的最近值=4 边(5,2)=4,如: closest3=2 不在生成树中的顶点3到树中最近的点是2 closest = _ , 0 , _, 2, 2, 2 ,普里姆(Prim)算法如下： #define INF 32767 /*INF表示*/ void Prim(int costMAXV, int n, int v) int lowcostMAXV, min;int closestMAXV,i,j,k;for (i=0;in;i+) /*给lowcos
28、t和closest置初值*/ lowcosti=costvi;closesti=v; ,for (i=1;in;i+) /*找出n-1个顶点*/ min=INF;for (j=0;jn;j+) /*在(V-U)中找出离U最近的顶点k*/if (lowcostj!=0 ,Prim()算法中有两重for循环,所以时间复杂度为O(n2)。,8.4.4 克鲁斯卡尔算法克鲁斯卡尔(Kruskal)算法是一种按权值的递增次序选择合适的边来构造最小生成树的方法。假设G=(V,E)是一个具有n个顶点的带权连通无向图,T=(U,TE)是G的最小生成树,则构造最小生成树的步骤如下：(1)置U的初值等于V(即包含有
29、G中的全部顶点),TE的初值为空集(即图T中每一个顶点都构成一个分量)。(2)将图G中的边按权值从小到大的顺序依次选取：若选取的边未使生成树T形成回路,则加入TE；否则舍弃,直到TE中包含(n-1)条边为止。,克鲁斯卡尔算法求解最小生成树的过程,为了简便,在实现克鲁斯卡尔算法Kruskal()时,参数E存放图G中的所有边,假设它们是按权值从小到大的顺序排列的。n为图G的顶点个数,e为图G的边数。typedef struct int u; /*边的起始顶点*/int v; /*边的终止顶点*/int w; /*边的权值*/ Edge; Edge Ee= 0,2,1, 3,5,2, 1,4,3,
30、2,5,4, 0,3,5, 2,3,5,1,2,5, 0,1,6, 2,4,6, 4,5,6 ;Kruskal()算法如下：,void Kruskal(Edge E ,int n,int e) / E已经从小到大排序 int i,j,m1,m2,sn1,sn2,k; int vsetMAXV;for (i=0;in;i+) vseti=i; /*初始化辅助数组*/k=1; /*k表示当前构造最小生成树的第几条边,初值为1*/j=0; /*E中边的下标,初值为0*/while (kn) /*生成的边数小于n时循环*/ m1=Ej.u; m2=Ej.v; /*取一条边的头尾顶点*/sn1=vset
31、m1;sn2=vsetm2; /*分别得到两个顶点所属的集合编号*/,if (sn1!=sn2) /*两顶点属于不同的集合,该边是最小生成树的一条边*/ printf(“(%d,%d):%dn“,m1,m2,Ej.w);k+; /*生成边数增1*/for (i=0;in;i+) /*两个集合统一编号*/if (vseti=sn2) /*集合编号为sn2的改为sn1*/vseti=sn1;j+; /*扫描下一条边*/,完整的克鲁斯卡尔算法应包括对边按权值递增排序，上述算法假设边已排序的情况下，时间复杂度为O(n2)。如果给定的带权连通无向图G有e条边,n个顶点，采用堆排序(在第11章中介绍)对边
32、按权值递增排序，那么用克鲁斯卡尔算法构造最小生成树的时间复杂度降为O(elog2e)。 由于它与n无关，只与e有关，所以说克鲁斯卡尔算法适合于稀疏图。,8.5 最短路径 8.5.1 路径的概念在一个无权的图中,若从一顶点到另一顶点存在着一条路径,则称该路径长度为该路径上所经过的边的数目,它等于该路径上的顶点数减1。由于从一顶点到另一顶点可能存在着多条路径,每条路径上所经过的边数可能不同,即路径长度不同,我们把路径长度最短(即经过的边数最少)的那条路径叫做最短路径,其路径长度叫做最短路径长度或最短距离。,对于带权的图,考虑路径上各边上的权值,则通常把一条路径上所经边的权值之和定义为该路径的路径长
33、度或称带权路径长度。从源点到终点可能不止一条路径,把带权路径长度最短的那条路径称为最短路径,其路径长度(权值之和)称为最短路径长度或者最短距离。,8.5.2 从一个顶点到其余各顶点的最短路径问题：给定一个带权有向图G与源点v,求从v到G中其他顶点的最短路径,并限定各边上的权值大于或等于0。,采用狄克斯特拉(Dijkstra)算法求解基本思想是：设G=(V,E)是一个带权有向图, 把图中顶点集合V分成两组：第一组为已求出最短路径的顶点集合(用S表示,初始时S中只有一个源点,以后每求得一条最短路径v,vk,就将vk加入到集合S中,直到全部顶点都加入到S中,算法就结束了)第二组为其余未确定最短路径的
34、顶点集合(用U表示)。按最短路径长度的递增次序依次把第二组的顶点加入S中。在加入的过程中,总保持从源点v到S中各顶点的最短路径长度不大于从源点v到U中任何顶点的最短路径长度。此外,每个顶点对应一个距离,S中的顶点的距离就是从v到此顶点的最短路径长度,U中的顶点的距离从v到此顶点只包括S中的顶点为中间顶点的当前最短路径长度。,狄克斯特拉算法的具体步骤如下：(1)初始时,S只包含源点,即S=v,v的距离为0。U包含除v外的其他顶点,U中顶点u距离为边上的权(若v与u有边)或(若u不是v的出边邻接点)。(2)从U中选取一个距离v最小的顶点k,把k加入S中(该选定的距离就是v到k的最短路径长度)。(3
35、)以k为新考虑的中间点,修改U中各顶点u的距离：若从源点v到顶点u(uU)的距离(经过顶点k)比原来距离(不经过顶点k)短,则修改顶点u的距离值,修改后的距离值的顶点k的距离加上边上的权。(4)重复步骤(2)和(3)直到所有顶点都包含在S中。,S U v0到06各顶点的距离dist (中间顶点只能是S中的点) 0 1,2,3,4,5,6 0, 4, 6, 6, , , 0,1 2,3,4,5,6 0, 4, 5, 6, 11, , 0,1,2 3,4,5,6 0, 4, 5, 6, 11, 9, 0,1,2,3 4,5,6 0, 4, 5, 6, 11, 9, 0,1,2,3,5 4,6 0,
36、 4, 5, 6, 10, 9, 17 0,1,2,3,5,4 6 0, 4, 5, 6, 10, 9, 16 0,1,2,3,5,4,6 0, 4, 5, 6, 10, 9, 16 则v0到v1v6各顶点的最短距离分别为4、5、6、10、9和16。,已求最短距离的顶点集,未知最短距离的顶点集,数组dist =从源点v0到每个点的最短路径长度(中间顶点只能是S中的点)disti=x 从源点v0到点vi 的最短路径长度x(中间顶点只能是S中的点) 数组Path =从源点v0到每个点i的最短路径中点i的直接前驱点 Path i=k 从源点v0到点i的最短路径中点i的直接前驱点k,例: 源点v0=0
37、dist6=1601-2-5-4-6path6=4, path4=5path5=2 path2=1path1=0 path0=-1,狄克斯特拉算法(n为图G的顶点数,v0为源点编号)： void Dijkstra(int cost MAXV,int n,int v0) int distMAXV,pathMAXV; int sMAXV;int mindis,i,j,u;for (i=0;in;i+) disti=costv0i; /*距离初始化*/si=0; /*s置空*/if (costv0iINF) /*路径初始化*/pathi=v0;elsepathi=-1;sv0=1;pathv0=0;
38、 /*源点编号v0放入s中*/,for (i=0;in;i+) /*循环直到所有顶点的最短路径都求出*/ mindis=INF;u=-1;for (j=0;jn;j+) /*选取不在s中且具有最小距离的顶点u*/if (sj=0 /*输出最短路径*/,void Ppath(int path,int i,int v0) /*输出从v0到i号结点路径的中间经过的点*/ int k;k=pathi;if (k=v0) return; /*找到了起点则返回*/Ppath(path,k,v0); /*输出从v0到k号结点路径的中间经过的点*/printf(“%d,“,k); /*输出k顶点*/ ,voi
39、d Dispath(int dist,int path,int s,int n,int v0) int i;for (i=0;in;i+)if (si=1) printf(“从%d到%d的最短路径长度为:%dt路径为:“,v0,i,disti);printf(“%d,“,v0); /*输出路径上的起点*/Ppath(path,i,v0); /*输出路径上的中间点*/printf(“%dn“,i); /*输出路径上的终点*/else printf(“从%d到%d不存在路径n“,v0,i); ,8.5.3 每对顶点之间的最短路径问题：对于一个各边权值均大于零的有向图,对每一对顶点vivj,求出vi
40、与vj之间的最短路径和最短路径长度。可以通过以每个顶点作为源点循环求出每对顶点之间的最短路径。除此之外,弗洛伊德(Floyd)算法也可用于求两顶点之间最短路径。,假设有向图G=(V,E)采用邻接矩阵cost存储,另外设置一个二维数组A用于存放当前顶点之间的最短路径长度,分量Aij表示当前顶点vi到顶点vj的最短路径长度。弗洛伊德算法的基本思想是递推产生一个矩阵序列A-1,A0,Ak,An-1,其中Akij表示从顶点vi到顶点vj的路径上所经过的中间顶点的编号不大于k的最短路径长度。,初始时,有A-1ij=costij。当求从顶点vi到顶点vj的路径上所经过的顶点编号不大于k+1的最短路径长度时
41、,要分两种情况考虑：一种情况是该路径不经过顶点编号为k+1的顶点,此时该路径长度与从顶点vi到顶点vj的路径上所经过的顶点编号不大于k的最短路径长度相同；另一种情况是从顶点vi到顶点vj的最短路径上经过编号为k+1的顶点。那么,该路径可分为两段,一段是从顶点vi到顶点vk+1的最短路径,另一段是从顶点vk+1到顶点vj的最短路径,此时最短路径长度等于这两段路径长度之和。这两种情况中的较小值,就是所要求的从顶点vi到顶点vj的路径上所经过的顶点编号不大于k+1的最短路径。,弗洛伊德思想可用如下的表达式来描述：A-1ij=costijAk+1ij=minAkij,Akik+1+Akk+1j (0k
42、 n-1)Pathk+1ij=从点i到点j的最短路径中点i的直接后继点,采用弗洛伊德算法求解过程,考虑顶点v0,A0ij表示由vi到vj,经由顶点v0的最短路径。只有从v2到v1经过v0的路径和从v2到v3经过v0的路径,不影响v2到v1和v2到v3的路径长度,因此,有：,考虑顶点v1,A1ij表示由vi到vj,经由顶点v1的最短路径。存在路径v0-v1-v2,路径长度为9,将A02改为9,path02改为1,因此,有：,考虑顶点v2,A2ij表示由vi到vj,经由顶点v2的最短路径。存在路径v3-v2-v0,长度为4,将A30改为4；存在路径v3-v2-v1,长度为4,将A31改为4。存在路
43、径v1-v2-v0,长度为7,将A10改为7。将path30、path31和path10均改为2。因此,有：,考虑顶点v3,A3ij表示由vi到vj,经由顶点v3的最短路径。存在路径v0-v3-v2,长度为8比原长度短,将A02改为8；存在路径v1-v3-v2-v0,长度为6(A13+A30=2+4=6)比原长度短,将A10改为6；存在路径v1-v3-v2,长度为3,比原长度短,将A12改为3；将path02、path10后path12均改为3。因此,有：,因此,最后求得的各顶点最短路径长度A和Path矩阵为：,从0到0路径为：0,0 路径长度为：0 从0到1路径为：0,1 路径长度为：5 从
44、0到2路径为：0,3,2 路径长度为：8 从0到3路径为：0,3 路径长度为：7 从1到0路径为：1,3,2,0 路径长度为：6 从1到1路径为：1,1 路径长度为：0 从1到2路径为：1,3,2 路径长度为：3 从1到3路径为：1,3 路径长度为：2,从2到0路径为：2,0 路径长度为：3 从2到1路径为：2,1 路径长度为：3 从2到2路径为：2,2 路径长度为：0 从2到3路径为：2,3 路径长度为：2 从3到0路径为：3,2,0 路径长度为：4 从3到1路径为：3,2,1 路径长度为：4 从3到2路径为：3,2 路径长度为：1 从3到3路径为：3,3 路径长度为：0,弗洛伊德算法如下：
45、 void Floyd(int costMAXV,int n) int AMAXVMAXV,pathMAXVMAXV;int i,j,k;for (i=0;i(Aik+Akj) Aij=Aik+Akj; pathij=k; Dispath(A, path, n); /*输出最短路径*/,8.6 拓扑排序设G=(V,E)是一个具有n个顶点的有向图,V中顶点序列v1,v2,vn称为一个拓扑序列,当且仅当该顶点序列满足下列条件：若是图中的边(即从顶点vi到vj有一条路径),则在序列中顶点vi必须排在顶点vj之前。在一个有向图中找一个拓扑序列的过程称为拓扑排序。例如,计算机专业的学生必须完成一系列规定的基础课和专业课才能毕业,假设这些课程的名称与相应代号有如下关系：,课程之间的先后关系可用有向图表示 ：,对这个有向图进行拓扑排序可得到一个拓扑序列：C1-C3-C2-C4-C7-C6-C5。也可得到另一个拓扑序列：C2-C7-C1-C3-C4-C5-C6,还可以得到其他的拓扑序列。学生按照任何一个拓扑序列都可以顺序地进行课程学习。,拓扑排序方法如下：(1)从有向图中选择一个没有前驱(即入度为0)的顶点并且输出它。(2)从有向图中删去该顶点,并且删去从该顶点发出的全部有向边。(3)重复上述两步,直到剩余的有向图中不再存在没有前驱(入度为0)的顶点为止。,
