数据结构—图的操作与应用—第9关：求图（邻接矩阵存储）最短路径的狄克斯特拉算法

熙云Sacire

已于 2024-06-02 23:07:46 修改

阅读量991

点赞数 22

文章标签：算法数据结构 c++ 图论

于 2024-06-02 23:06:34 首次发布

本文链接：https://blog.csdn.net/xiyunSacire/article/details/139399396

版权

任务描述

本关任务：图的存储结构为邻接矩阵，要求编写函数实现狄克斯特拉算法。

编程要求

网G的存储结构为邻接矩阵，编写函数利用狄克斯特拉（Dijkstra））求图的单源最短路径:

void Dijkstra(MGraph g,int v); //求从v到其他顶点的最短路径

测试说明

平台会对你编写的代码进行测试：

测试输入： 1 lt4.txt

输入说明：第一行输入1，表示输入图的类型为有向网。第二行输入文件名，该文件里保存了图的数据信息，内容如下： 7 12 0 1 2 3 4 5 6 0 1 4 0 2 6 0 3 6 1 2 1 1 4 7 2 4 6 2 5 4 3 2 2 3 5 5 4 6 6 5 4 1 5 6 8 第1行为图的顶点的个数n；第2行为图的边的条数m；第3行至第n+2行是n个顶点的数据；第n+3行至第n+m+2行是m条边的数据；

预期输出： 有向网 7个顶点12条边。顶点依次是: 0 1 2 3 4 5 6 图的邻接矩阵: ∞ 4 6 6 ∞ ∞ ∞ ∞ ∞ 1 ∞ 7 ∞ ∞
∞ ∞ ∞ ∞ 6 4 ∞
∞ ∞ 2 ∞ ∞ 5 ∞
∞ ∞ ∞ ∞ ∞ ∞ 6
∞ ∞ ∞ ∞ 1 ∞ 8 ∞ ∞ ∞ ∞ ∞ ∞ ∞
dist: ∞ 4 6 6 ∞ ∞ ∞
path: -1 0 0 0 -1 -1 -1
dist: ∞ 4 5 6 11 ∞ ∞
path: -1 0 1 0 1 -1 -1
dist: ∞ 4 5 6 11 9 ∞
path: -1 0 1 0 1 2 -1
dist: ∞ 4 5 6 11 9 ∞
path: -1 0 1 0 1 2 -1
dist: ∞ 4 5 6 10 9 17 path: -1 0 1 0 5 2 5
dist: ∞ 4 5 6 10 9 16 path: -1 0 1 0 5 2 4 dist: ∞ 4 5 6 10 9 16 path: -1 0 1 0 5 2 4
从0到1最短路径长度为:4 0→1 从0到2最短路径长度为:5 0→1→2 从0到3最短路径长度为:6 0→3 从0到4最短路径长度为:10 0→1→2→5→4 从0到5最短路径长度为:9 0→1→2→5 从0到6最短路径长度为:16 0→1→2→5→4→6

输出说明：第一行输出图的类型。第二部分起输出图的顶点和边的数据信息。第三部分输出辅助数组的变化过程。第四部分输出从起点到其余各顶点的最短路径。

#include<stdio.h> 
#include<stdlib.h> 
#include<string.h>
#include<limits.h>  
#include<iostream>
using namespace std;
#include"MGraph.h"


void Dijkstra(MGraph g,int v);		//求从v到其他顶点的最短路径
void DispAllPath(MGraph &g,int dist[],int path[],int S[],int v) ;//输出从顶点v出发的所有最短路径
void Dispdistpath(int dist[],int path[],int n);	//输出dist数组和path数组

int main()
{
	MGraph g;
	int i,j,n;
	CreateGraphF(g);           /* 利用数据文件创建有向图*/
	Display(g);	/* 输出有向图*/  
	Dijkstra(g,0);	
	return 0;
}

void Dijkstra(MGraph g,int v)		
{	
    //求从v到其他顶点的最短路径
    /********** Begin **********/
    int dist[MAX_VERTEX_NUM];
    int path[MAX_VERTEX_NUM];
    int S[MAX_VERTEX_NUM];
    int mindis,i,j,u=0;
    for (i=0;i<g.vexnum;i++)
    {
        dist[i]=g.arcs[v][i].adj;
        S[i]=0;
        if (g.arcs[v][i].adj < INFINITY)
            path[i]=v;
        else
            path[i]=-1;
    }
    Dispdistpath(dist,path,g.vexnum);
    S[v]=1;
    for (i=0;i<g.vexnum-1;i++)
    {
        mindis=INFINITY;
        for(j=0;j<g.vexnum;j++)
            if(S[j]==0&&dist[j]<mindis){
                u=j;
                mindis=dist[j];
            }
        S[u]=1;
        for(j=0;j<g.vexnum;j++)
            if(S[j]==0)
                if(g.arcs[u][j].adj < INFINITY &&dist[u]+g.arcs[u][j].adj<dist[j]){
                    dist[j]=dist[u]+g.arcs[u][j].adj;
                    path[j]=u;
                }
        Dispdistpath(dist,path,g.vexnum);
    }
    DispAllPath(g,dist,path,S,v);
    /********** End **********/
}

void DispAllPath(MGraph &g,int dist[],int path[],int S[],int v) //输出从顶点v出发的所有最短路径
{
	int i,j,k,count=0;
	int apath[MAX_VERTEX_NUM],d;				//存放一条最短路径(逆向)及其顶点个数
	for (i=0;i<g.vexnum;i++)
		if (path[i]!=-1)
			count++;
	if (count==1)						//path中只有一个不为-1时表示没有路径
	{	printf("从指定的顶点到其他顶点都没有路径!!!\n");
		return;
	}
	for (i=0;i<g.vexnum;i++)					//循环输出从顶点v到i的路径
		if (S[i]==1 && i!=v)
		{
			//printf("从%s到%s最短路径长度为:%s\t路径:",g.vexs [v],g.vexs[i],dist[i]);
			cout<<"从"<<g.vexs [v]<<"到"<<g.vexs[i]<<"最短路径长度为:"<<dist[i]<<"\t";
			d=0; 			apath[d]=i;			//添加路径上的终点
			k=path[i];
			if (k==-1)					//没有路径的情况
				printf("无路径\n");
			else						//存在路径时输出该路径
			{	while (k!=v)
				{	d++; apath[d]=k;
					k=path[k];
				}
				d++; apath[d]=v;		//添加路径上的起点
				//printf("%d",apath[d]);	//先输出起点
				cout<<g.vexs [ apath[d] ];
				for (j=d-1;j>=0;j--)	//再输出其他顶点
					//printf("→%d",apath[j]);
					cout<<"→"<<g.vexs [ apath[j] ];
				printf("\n");
			}
		}
}


void Dispdistpath(int dist[],int path[],int n)	//输出dist数组和path数组
{
	int i;
	printf("dist:\t");
	for (i=0;i<n;i++)
		if (dist[i]==INFINITY)
			printf("%s\t","~");
		else
			printf("%d\t",dist[i]);
	printf("\n");
	printf("path:\t");
	for (i=0;i<n;i++)
		printf("%d\t",path[i]);
	printf("\n");
}