数据结构课程设计Word格式文档下载.docx
《数据结构课程设计Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《数据结构课程设计Word格式文档下载.docx(20页珍藏版)》请在冰豆网上搜索。
。
G的邻接矩阵
当邻接矩阵的行表头、列表头顺序一定时,一个图的邻接矩阵表示是唯一的。
图的邻接矩阵表示,除了需用一个二维数组存储顶点之间的相邻关系的邻接矩阵外,通常还需要使用一个具有n个元素的一维数组来存储顶点信息,其中下标为i的元素存储顶点i的信息。
因此,图的邻接矩阵的存储结构定义如下:
#definfMVNum50//最大顶点数
typedefstruct
{
VertexTypevexs[MVNum];
//顶点数组,类型假定为char型
Adjmatrixarcs[MVNum][MVNum];
//邻接矩阵,假定为int型
}MGraph;
1.1.2单源最短路径
最短路径的提法很多。
在这里先讨论单源最短路径问题:
即已知有向图(带权),我们希望找出从某个源点S
V到G中其余各顶点的最短路径。
为了叙述方便,我们把路径上的开始点称为源点,路径的最后一个顶点为终点。
那么,如何求得给定有向图的单源最短路径呢?
迪杰斯特拉(Dijkstra)提出按路径长度递增产生诸点的最短路径算法,称之为迪杰斯特拉算法。
迪杰斯特拉算法求最短路径的实现思想是:
设G=(V,E)是一个有向图,结点集为,
,cost是表示G的邻接矩阵,cost[i][j]表示有向边<
i,j>
的权。
若不存在有向边<
,则cost[i][j]的权为无穷大(这里取值为32767)。
设S是一个集合,其中的每个元素表示一个顶点,从源点到这些顶点的最短距离已经求出。
设顶点v1为源点,集合S的初态只包含一个元素,即顶点v1。
数组dist记录从源点到其他顶点当前的最短距离,其初值为dist[i]=cost[v1][i],i=1,2,……,n。
从S之外的顶点集合V-S中选出一个顶点w,使dist[w]的值最小。
于是从源点到达w只通过S中顶点,把w加入集合S中,调整dist中记录的从源点到V-S中每个顶点v的距离:
从原来的dist[v]和dist[w]+cost[w][v]中选择较小的值作为新的dist[v]。
重复上述过程,直到V-S为空。
最终结果是:
S记录了从源点到该顶点存在最短路径的顶点集合,数组dist记录了源点到V中其余各顶点之间的最短路径,path是最短路径的路径数组,其中path[i]表示从源点到顶点i之间的最短路径的前驱顶点。
因此,迪杰斯特拉算法可用自然语言描述如下:
初始化S和D,置空最短路径终点集,置初始的最短路径值;
S[v1]=TRUE;
D[v1]=0;
//S集初始时只有源点,源点到源点的距离为0;
While(S集中顶点数<
n)
{
开始循环,每次求得v1到某个v顶点的最短路径,并加v到S集中;
S[v]=TRUE;
更新当前最短路径及距离;
}
1.1.3任意一对顶点间最短路径
任意一对顶点间最短路径问题,是对于给定的有向网络图G=(V,E),要对G中任意一对顶点有序对“v,w(v
w)”,找出v到w的最短路径。
要解决这个问题,我们可以依次把有向网络图中每个顶点作为源点,重复执行前面讨论的迪杰斯特拉算法n次,即可以求得每对顶点之间的最短路径。
这里还可以用另外一种方法,称作费洛伊德(Floyd)算法。
费洛伊德(Floyd)算法算法的基本思想是:
假设求从顶点vi到vj的最短路径。
如果从vi到vj存在一条长度为arcs[i][j]的路径,该路径不一定是最短路径,还需要进行n次试探。
首先考虑路径<
vi,v1>
和<
v1,vj>
是否存在。
如果存在,则比较<
vi,vj>
vi,v1,vj>
的路径长度,取长度较短者为当前所求得的最短路径。
该路径是中间顶点序号不大于1的最短路径。
其次,考虑从vi到vj是否包含有顶点v2为中间顶点的路径<
vi,…,v2,…,vj>
,若没有,则说明从vi到vj的当前最短路径就是前一步求出的;
若有,那么<
可分解为<
vi,…v2>
v2,…,vj>
而这两条路径是前一次找到的中间顶点序号不大于1的最短路径,将这两条路径长度相加就得到路径<
的长度。
将该长度与前一次中求出的从vi到vj的中间顶点序号不大于1的最短路径比较,取其长度较短者作为当前求得的从vi到vj的中间顶点序号不大于2的最短路径。
依此类推,直到顶点vn加入当前从vi到vj的最短路径后,选出从vi到vj的中间顶点序号不大于n的最短路径为止。
由于图G中顶点序号不大于n,所以vi到vj的中间顶点序号不大于n的最短路径,已考虑了所有顶点作为中间顶点的可能性,因此,它就是vi到vj的最短路径。
1.2程序流程图
二、详细设计
2.1建立有向图的存储结构
voidCreateMGraph(MGraph*G,intn,inte)
inti,j,k,w;
for(i=1;
i<
=n;
i++)
G->
vexs[i]=(char)i;
for(j=1;
j<
j++)
arcs[i][j]=Maxint;
printf("
输入%d条边的i,j及w:
\n"
e);
for(k=1;
k<
=e;
k++)
{
scanf("
%d,%d,%d"
&
i,&
j,&
w);
arcs[i][j]=w;
}
有向图建立完毕\n"
);
2.2迪杰斯特拉算法
voidDijkstra(MGraph*G,intv1,intn)
intD2[MVNum],P2[MVNum];
intv,i,w,min;
enumbooleanS[MVNum];
for(v=1;
v<
v++)
S[v]=FALSE;
D2[v]=G->
arcs[v1][v];
if(D2[v]<
Maxint)
P2[v]=v1;
else
P2[v]=0;
D2[v1]=0;
for(i=2;
n;
min=Maxint;
for(w=1;
w<
w++)
if(!
S[w]&
&
D2[w]<
min)
v=w;
min=D2[w];
S[v]=TRUE;
(D2[v]+G->
arcs[v][w]<
D2[w]))
D2[w]=D2[v]+G->
arcs[v][w];
P2[w]=v;
printf("
路径长度路径\n"
%5d"
D2[i]);
i);
v=P2[i];
while(v!
=0)
<
-%d"
v);
v=P2[v];
2.3费洛伊德算法
voidFloyd(MGraph*G,intn)
{inti,j,k,v,w;
{if(G->
arcs[i][j]!
=Maxint)
P[i][j]=j;
P[i][j]=0;
D[i][j]=G->
arcs[i][j];
if(D[i][k]+D[k][j]<
D[i][j])
D[i][j]=D[i][k]+D[k][j];
P[i][j]=P[i][k];
;
2.4运行主控程序
voidmain()
{MGraph*G;
intm,n,e,v,w,k;
intxz=1;
G=(MGraph*)malloc(sizeof(MGraph));
输入图中顶点个数和边数n,e:
"
%d,%d"
n,&
e);
CreateMGraph(G,n,e);
while(xz!
{printf("
******求城市间的最短路径******\n"
1.求一个城市到所有城市的最短路径\n"
printf(“2.求任意的两个城市之间的最短路径\n"
请选择:
1或2,选择0退出:
%d"
xz);
if(xz==2)
Floyd(G,n);
输入起点和终点:
v,w:
v,&
k=P[v][w];
if(k==0)
顶点%d到%d无路径!
v,w);
从顶点%d到%d的最短路径是:
:
v,w,v);
while(k!
=w)
{
→¨
²
k);
k=P[k][w];
w);
路径长度:
%d\n"
D[v][w]);
elseif(xz==1)
求单源路径,输入源点v:
v);
Dijkstra(G,v,n);
结束求最短路径"
三、调试分析
编译:
在第一次编译时出现了很多错误,是因为我对C语言的不熟练,比如调用费洛伊德算法时出现了调用的错误,找了好久,才改正过来,还有就是for语句的运用,由于本次程序要用很多for循环,我把一次for循环放到了上面for循环中,导致程序不能正确输出结果。
最后把调到外面才OK。
四、测试结果
交通咨询系统通过迪杰斯特拉算法得出每一个城市到所有城市的最短路径,然后通过费洛伊德算法得出任意两个城市间最短路径。
本程序的运行环境为DOS操作系统,执行文件为交通咨询.exe。
进入程序后,即显示文本的操作界面:
输入定点数和边数后,依次输入有向边,以及两点间的距离,用逗号间隔开,回车换行,当有向图输入完全时,进入菜单
选择你要进行的模式,求一个城市到所有城市的最短路径输入1,回车进入该模式,输入你所要求的城市代表的顶点数,回车确认。
求任意的两个城市之间的最短路径输入2,回车确认进入该模式,输入两座城市代表的顶点数,以逗号间隔开,回车确认。
运行结果:
下面是城市交通图
五、附录
#include<
stdio.h>
stdlib.h>
#defineMVNum100//最大顶点数
#defineMaxint35000
enumboolean{FALSE,TRUE};
typedefcharVertextype;
typedefintAdjmatrix;
typedefstruct
{
Vertextypevexs[MVNum];
//顶点数组,类型假定为char型
Adjmatrixarcs[MVNum][MVNum];
//邻接矩阵,假定为int型
}MGraph;
intD1[MVNum],p1[MVNum];
intD[MVNum][MVNum],p[MVNum][MVNum];
//文件名save.c
voidCreateMGraph(MGraph*G,intn,inte)
{//采用邻接矩阵表示法构造有向图G,n,e表示图的当前顶点数和边数
i++)//输入顶点信息
G->
//初始化邻接矩阵
printf("
\n"
{//读入e条边,建立邻接矩阵
scanf("
有向图的存储结构建立完毕!
//文件名:
dijkstra.c(迪杰斯特拉算法)
voidDijkstra(MGraph*G,intv1,intn)
{//用Dijkstra算法求有向图G的v1顶点到其他顶点v的最短路径p[v]及其权D[v]
//设G是有向图的邻接矩阵,若边<
不存在,则G[i][j]=Maxint
//S[v]为真当且仅当v属于S,及以求的从v1到v的最短路径
intD2[MVNum],p2[MVNum];
intv,i,w,min;
enumbooleanS[MVNum];
for(v=1;
{//初始化S和D
//置空最短路径终点集
D2[v]=G->
//置初始的最短路径值
if(D2[v]<
Maxint)
p2[v]=v1;
//v1是的前趋(双亲)
else
p2[v]=0;
//v无前趋
}//End_for
D2[v1]=0;
//S集初始时只有源点,源点到源点的距离为0
//开始循环,每次求的V1到某个V顶点的最短路径,并加V到S集中
for(i=2;
{//其余n-1个顶点
//当前所知离v1顶点的最近距离,设初值为∞
w++)//对所有顶点检查
if(!
S[w]&
D2[w]<
{//找离v1最近的顶点w,并将其赋给v,距离赋给min
v=w;
//在S集之外的离v1最近的顶点序号
min=D2[w];
//最近的距离
}//W顶点距离V1顶点更近
S[v]=TRUE;
//将v并入S集
for(w=1;
w++)//更新当前最短路径及距离
if(!
{//修改D2[w]和p2[w],w属于V-S
D2[w]=D2[v]+G->
//更新D2[w]
p2[w]=v;
}//End_if
}//End_for
{printf("
D2[i]);
i);
v=p2[i];
=0){
printf("
v);
v=p2[v];
//文件名floyd.c(费洛伊德算法)
voidFloyd(MGraph*G,intn)
inti,j,k;
i++)//设置路径长度D和路径path初值
for(j=1;
if(G->
p[i][j]=j;
//j是i的后继
else
p[i][j]=0;
D[i][j]=G->
}
k++){
{//做K次迭代,每次均试图将顶点K扩充到当前求得的从i到j的最短路径pij上
for(i=1;
for(j=1;
{if(D[i][k]+D[k][j]<
D[i][j]){
D[i][j]=D[i][k]+D[k][j];
//修改长度
p[i][j]=p[i][k];
}
{MGraph*G;
intn,e,v,w,k;
"
&
n,&
CreateMGraph(G,n,e);
//建立图的存储结构
while(xz!
求城市之间的最断路径\n"
printf("
-------------------------------\n"
2.求任意的两个城市之间的最短路径\n"
--------------------------------\n"
1或2,选择0:
退出:
if(xz==2)
{
Floyd(G,n);
//调用费洛伊德算法
printf("
v,w:
w);
k=p[v][w];
//k为起点v的后继顶点
if(k==0)
printf("
顶点%d到 %d无路径!
{
%d"
}
while(k!
-->
//输出后继顶点
k=p[k][w];
//继续找下一个后继顶点
//输出终点w
路径长度:
if(xz==1)
求单源路径,输入起点 v:
Dijkstra(G,v,n);
//调用迪杰斯特拉算法
结束求最短路径,再见!
六、心得体会
课程设计是对自己一个学期以来学习的良好检验,在所有课题中,我选择了交通咨询系统的课题题目。
在一个星期的设计过程中,碰到过许多困难,也从书籍和网络中找过许多相关资料来完善自己的设计。
但结果不是让人很满意,可能是接触时间太短,对数据结构的认识还不深刻。
通过这次课程设计,我对用迪杰斯特拉算法计算最短路径有了更深刻的认识,加强了我对数据结构的了解。
我知道如果要想在C方面有成就,必须多想多练习。
现在我还是缺少练习,看程序可能看得懂。
而当真正的编写程序起来就遇到困难了。
通过这次系统的编写,我明白了编程要细心,有耐心。
有些程序看起来很多很难,但只要抽丝剥茧的慢慢分析,就能明白它是怎么运行的。
其实,只要踏踏实实一步一步来,你会发觉原来编程难度也不过如此。
同时,我懂得了要多思考,遇到问题或不解时不能不懂装懂,要尽量去想请清其缘由,要懂得查阅有用的资料。
我知道在这次的课程设计中自已还有许多的欠缺,在以后的日子中,我一定得更加严格要求自己,改正缺点,不断努力,不断进步。