NOJ数据结构实验009——用迪杰斯特拉算法求赋权图中的最短路径

在这里插入图片描述

本题是以上一道题为基础的提高,但我们需要求出给出的最短路径,需要多设置一个数组TargetPath记录每一点的最优前驱(即在源点到该点的最短路径上该点前面紧邻的那个点)。

代码:

#include <stdio.h>
#include <stdlib.h>

#define MAX_VER_NUM 20


typedef struct {
	int vertex[MAX_VER_NUM];
	int weight[MAX_VER_NUM][MAX_VER_NUM];
	int num_vex;//弧的个数用不上
}AdjMatrix;

typedef struct {
	int visited[MAX_VER_NUM];
	int length[MAX_VER_NUM];
}Path;

int TargetPath[MAX_VER_NUM];//记录所求路径

AdjMatrix* CreateMatrix(int n);   //构建邻接矩阵
void InitPath(AdjMatrix A, Path* dist, int start);   //初始化dist数组
int JudgePath(AdjMatrix A, Path* dist);   //判断是否遍历完成
int FindMinPath(AdjMatrix A, Path* dist);     //寻找最短路径的位置
void UpdatePath(AdjMatrix A, Path* dist, int pos_minpath);   //更新路径
void dijsktra(AdjMatrix A, Path* dist, int start);
//dijsktra算法并将所需路径存放在target中
void PrintDist(AdjMatrix A, int start, int end);   //打印数据

int main()
{
	int n, start, end;
	scanf("%d", &n);
	AdjMatrix* A;
	A = CreateMatrix(n);
	scanf("%d%d", &start, &end);
	Path dist;
	dijsktra(*A, &dist, start);
	//Print(A,start,end);
	PrintDist(*A, start, end);
	return 0;
}

AdjMatrix* CreateMatrix(int n)
{   //构建邻接矩阵
	AdjMatrix* A;
	A = (AdjMatrix*)malloc(sizeof(AdjMatrix));
	if (A == NULL) {
		return NULL;
	}
	A->num_vex = n;
	for (int i = 0; i < A->num_vex; i++) {
		for (int j = 0; j < A->num_vex; j++) {
			scanf("%d", &A->weight[i][j]);
		}
	}
	return A;
}

void InitPath(AdjMatrix A, Path* dist, int start)
{   //初始化dist数组
	for (int i = 0; i < A.num_vex; i++) {
		dist->visited[i] = -1;
		dist->length[i] = 10000;
	}
	for (int i = 0; i < A.num_vex; i++) {
		if (A.weight[start][i] != 10000) {
			dist->visited[i] = 1;
			dist->length[i] = A.weight[start][i];
			TargetPath[i]=start;
		}
	}
	dist->visited[start] = 0;
}

int JudgePath(AdjMatrix A, Path* dist)
{   //判断是否遍历完成
	for (int i = 0; i < A.num_vex; i++) {
		if (dist->visited[i] == 1) {
			return 1;
		}
	}
	return 0;
}

int FindMinPath(AdjMatrix A, Path* dist)
{   //寻找最短路径的位置
	int min = 10000, tmp = -1;
	for (int i = 0; i < A.num_vex; i++) {
		if (dist->visited[i] == 1) {
			if (dist->length[i] < min) {
				min = dist->length[i];
				tmp = i;
			}
		}
	}
	return tmp;
}

void UpdatePath(AdjMatrix A, Path* dist, int pos_minpath)
{   //更新路径
	int minpath = dist->length[pos_minpath];
	dist->visited[pos_minpath] = 0;
	for (int i = 0; i < A.num_vex; i++) {
		if (A.weight[pos_minpath][i] != 10000) {
			if (dist->visited[i] == 1) {
				if (minpath + A.weight[pos_minpath][i] < dist->length[i]) {
					dist->length[i] = minpath + A.weight[pos_minpath][i];
					TargetPath[i] = pos_minpath;
				}
			}
			else if (dist->visited[i] == -1) {
				dist->visited[i] = 1;
				dist->length[i] = A.weight[pos_minpath][i] + minpath;
				TargetPath[i] = pos_minpath;
			}
		}
	}
}

void dijsktra(AdjMatrix A, Path* dist, int start)
{   //dijsktra算法并将所需路径存放在a中
	int pos_path;
	InitPath(A, dist, start);
	while (JudgePath(A, dist)) {
		pos_path = FindMinPath(A, dist);
		UpdatePath(A, dist, pos_path);

	}
}

void PrintDist(AdjMatrix A, int start, int end)
{   //打印数据
	int i, mid = end, tmp[MAX_VER_NUM];
	for (i = 0;; i++) {
		if (mid == start) {
			tmp[i] = start;
			break;
		}
		else {
			tmp[i] = mid;
			mid = TargetPath[mid];
		}
	}
	for (; i >= 0; i--) {
		printf("%d\n", tmp[i]);
	}
}
  • 4
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值