迪杰斯特拉算法(dijkstra)

迪杰斯特拉算法求单源点到其他顶点间的最短路径

/* 模板 */
int d[maxn];   //起点到其他顶点的最短路径长度
bool vis[maxn]; //标记数组
void dijkstra(int s){
	fill(d,d+maxn,INF);
	d[s] = 0;
	for(int i=0;i<n;i++){
		int min = INF;
		int temp;
		for(int j=0;j<n;j++){ //找最小保存到temp 
			if(vis[j]==false && d[j]<min){
				min = d[j];
				temp = j;
			} 
		}
		vis[temp] = true;
		for(int v=0;v<n;v++){ //更新 
			if(vis[v]==false && d[v]>d[temp]+road[temp][v]) {
				d[v] = d[temp]+road[temp][v];	
			}
		}
	}
}

题1:旅游规划

有了一张自驾旅游路线图,你会知道城市间的高速公路长度、以及该公路要收取的过路费。现在需要你写一个程序,帮助前来咨询的游客找一条出发地和目的地之间的最短路径。如果有若干条路径都是最短的,那么需要输出最便宜的一条路径。

输入格式:

输入说明:输入数据的第1行给出4个正整数N、M、S、D,其中N(2≤N≤500)是城市的个数,顺便假设城市的编号为0~(N−1);M是高速公路的条数;S是出发地的城市编号;D是目的地的城市编号。随后的M行中,每行给出一条高速公路的信息,分别是:城市1、城市2、高速公路长度、收费额,中间用空格分开,数字均为整数且不超过500。输入保证解的存在。

输出格式:

在一行里输出路径的长度和收费总额,数字间以空格分隔,输出结尾不能有多余空格。

输入样例:

4 5 0 3
0 1 1 20
1 3 2 30
0 3 4 10
0 2 2 20
2 3 1 20

输出样例:

3 40
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn = 501;
const int INF = 501;
int road[maxn][maxn];
int cost[maxn][maxn];
int d[maxn]; //最短距离 
int money[maxn]; //最少花费 
bool vis[maxn]={false}; //结点是否访问过 
int n;
void dijkstra(int s) //s为起点 
{
	fill(d,d+maxn,INF);  //初始化 
	fill(money,money+maxn,INF);
	d[s] = 0;
	money[s] = 0;
	for(int i=0;i<n;i++){   //循环n次 
		int u = -1,MIN = INF;  
		for(int j=0;j<n;j++){  //找到未访问结点中最小的   
			if(vis[j] == false && d[j] < MIN){
				u = j;
				MIN = d[j];
			}
		}
		if(u == -1) return ;
		vis[u] = true;			//标记为已访问 
		for(int v=0;v<n;v++){  //优化的d[v] 
			if(vis[v] == false && road[u][v] != INF && d[u] + road[u][v] < d[v]){
				d[v] = d[u] + road[u][v];
				money[v] = money[u] + cost[u][v];
			}
			else if(vis[v] == false && road[u][v] != INF && d[u] + road[u][v] == d[v]){
				if(money[u] + cost[u][v] < money[v]){
					  money[v] = money[u] + cost[u][v];
				}
			}
		}
	}
}
int main()
{
	int m,s,end;
	int v1,v2,length,c;
	cin>>n>>m>>s>>end;
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			if(i==j) road[i][j] = 0;
			else road[i][j] = maxn; 
		}
	}
	for(int i=0;i<m;i++){
		cin>>v1>>v2>>length>>c;
		road[v1][v2] = road[v2][v1] = length;
		cost[v1][v2] = cost[v2][v1] = c;
	} 
	dijkstra(s);
	cout<<d[end]<<" "<<money[end];
	return 0;	
} 

题2: 紧急救援

题目链接:https://pintia.cn/problem-sets/994805046380707840/problems/994805073643683840

作为一个城市的应急救援队伍的负责人,你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候,你的任务是带领你的救援队尽快赶往事发地,同时,一路上召集尽可能多的救援队。

输入格式:

输入第一行给出4个正整数N、M、S、D,其中N(2≤N≤500)是城市的个数,顺便假设城市的编号为0 ~ (N−1);M是快速道路的条数;S是出发地的城市编号;D是目的地的城市编号。

第二行给出N个正整数,其中第i个数是第i个城市的救援队的数目,数字间以空格分隔。随后的M行中,每行给出一条快速道路的信息,分别是:城市1、城市2、快速道路的长度,中间用空格分开,数字均为整数且不超过500。输入保证救援可行且最优解唯一。

输出格式:

第一行输出最短路径的条数和能够召集的最多的救援队数量。第二行输出从S到D的路径中经过的城市编号。数字间以空格分隔,输出结尾不能有多余空格。

输入样例:

4 5 0 3
20 30 40 10
0 1 1
1 3 2
0 3 3
0 2 2
2 3 2

输出样例:

2 60
0 1 3

注意:输出的是最短路径的条数,而不是最短路径的边数

#include<bits/stdc++.h>
using namespace std;
const int maxn = 501;
int a[maxn][maxn];
int b[maxn];
bool vis[maxn]={false};
int n,M,S,D,d[maxn],j[maxn]={0}; //d数组保存最短路径 ,j数组保存救援队数量 
int pre[maxn],cnt[maxn];   //pre保存是谁到达的该顶点,cnt保存该顶点的最短路径条数 
void init()
{
	fill(d,d+n,0x3f3f3f);
	fill(pre,pre+n,-1);
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			if(i!=j) {
				a[i][j] = a[j][i] = maxn;
			}else{
				a[i][j] = a[j][i] = 0;
			}
		}
	}
} 
void djkstla(){   
	cnt[S] = 1;  //开头默认有一条 
	j[S] = b[S]; //初始化 
	d[S] = 0;    //初始化到S为0  
	for(int i=0;i<n;i++){
		int u = -1,min = maxn;
		for(int j=0;j<n;j++){  //d中找最小 
			if(min>d[j] && vis[j]==false){
				min = d[j];
				u = j; 
			}
		}
		vis[u] = true;
		if(u == -1) return ;
		for(int v=0;v<n;v++){
			if(vis[v]==false && d[v] > d[u] + a[u][v]){ //若大于则更新d数组 
				d[v] = d[u] + a[u][v];
				j[v] = j[u] + b[v];
				pre[v] = u; //保存u到达的该顶点
				cnt[v] = cnt[u]; 
			}
			else if(vis[v]==false && d[v] == d[u] + a[u][v]){  
				cnt[v] = cnt[u] + cnt[v];  //若等于,就更新,更新时也要加上到达中转点v的那几条 
				if(j[v] < j[u] + b[v]){ 
					j[v] = j[u] + b[v];
					pre[v] = u;
				}
			}
		}
	} 
	cout<<cnt[D]<<" "<<j[D]<<endl;
	vector<int> t; //保存路径 
	int v = D; 
	t.push_back(D);
	while(1){  
		v = pre[v];
		t.push_back(v);
		if(v==S) break;
	}
	for(int i=t.size()-1;i>=0;i--){ //打印路径 
		printf("%d",t[i]);
		if(i!=0) printf(" ");
	}
}
int main()
{
	cin >> n >> M >> S >> D;
	init();
	for(int i=0;i<n;i++){
		cin>>b[i];
	}
	int n1,n2,n3;
	for(int i=0;i<M;i++){
		cin>>n1>>n2>>n3;
		a[n1][n2] = a[n2][n1] = n3;
	}
	djkstla();
	return 0;
}

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值