最短路学习总结(一)

图论这个东东怪难的,写个博客总结一下吧。
首先先贴一个图论的学习路线吧:https://blog.csdn.net/Fire_to_cheat_/article/details/80028763

再贴上从别处收集来的题:https://blog.csdn.net/cxb569262726/article/details/7859515(把这些刷完你就无敌了)


部分图片和文字来自:https://www.bilibili.com/video/av64263196/?p=1

最短路问题就是给定起点和终点,在途中找到一条路径使得沿此路径上各边的权值总和最小。

解决最短路问题我们有四种算法:

Dijkstra算法,Bellman Ford算法,SPFA算法,Floyd算法

在这里插入图片描述

Dijkstra算法:

Dijstra算法思想:设有两个顶点集合S和T,集合S中存放的图中已经找到最短路径的顶点,集合T中存放图中剩余顶点。初始状态的时候,集合S中只包含源点v0,然后不断从集合T中选取到顶点v0路径长度最短的顶点v并入到集合S中去,集合S每并入一个新的顶点,都要修改顶点v0到集合T中顶点的最短路径长度值。

算法实现:
1.设置两个顶点的集合S和T
a) S中存放已找 到最短路径的顶点,初始时,集合S只有 一个顶点,
即源点v0.
b)下中存放当前还未找到最短路径的顶点
2.在工集合中选取当前最短的-条最短路径(v0, … vk), 从而将
vk加入到顶点集合S中,并且修改源点v0到T中各顶点的最短路径长度。
3.重复步骤2,直到所有的顶点都被加入集合S中。

我们根据图来解释如何进行操作:

在这里插入图片描述
1.首先求出长度最短的一条最短路径: 0 -> 2,长度为5 (2号点)。
2.顶点2的最短路求出来之后,0到其他顶点的最短路径长度有可能要改变,比如0->1从原来的∞变成了20,0->5从原来的∞变成了12.这样次短的最短路就是12 (5号点) .
3.5的最短路求出后,0到其他点的最短路还会发生改变。比如0 ->3从原来的30变为了22,0->4也从∞变为了28。这时第三短的最短路就确定了,是20 (1号点) .
4.之后再依次确定第四短的最短路是22 (3号点) .
5.第五短的最短路是28 (4号点).

下面我们做个板子题练练手:
传送门:https://vjudge.net/problem/HDU-2544

AC代码如下:

//Dijstra算法不加堆优化,时间复杂度O(n^2),利用邻接矩阵存图
#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<list>
#include<cmath>
#include<algorithm>
#define lson node<<1,st,mid
#define rson node<<1|1,mid+1,ed
#define mem(a,x) memset(a,x,sizeof(a))
#define me(a) memset(a,0,sizeof(a))
#define IOS ios::sync_with_stdio(false)
#define lowbit(x) x&(-x)
#define up(i,x,y) for(long long i=x;i<y;i++)
#define down(i,x,y) for(long long i=x;i>=y;i--)
#define in freopen("in.txt","r",stdin)
#define out freopen("out.txt","w",stdout)  	
typedef long long ll;
const ll mod=1e9+7;
const ll INF=0x3f3f3f3f;
const int maxn=10005;
const double pi = acos(-1.0);
using namespace std; 
ll n,m;
ll edge[maxn][maxn],dist[maxn];
bool vis[maxn];
void dijstra(ll u){
	mem(dist,INF);
	for(ll i=1;i<=n;i++){
		vis[i]=false;
		dist[i]=edge[u][i];
	} 
	vis[u]=true;
	for(ll i=1;i<n;i++){
		ll min=INF,v=u;
		for(ll j=2;j<=n;j++){
			if(!vis[j]&&dist[j]<min){
				min=dist[j];
				v=j;
			}
		}
		vis[v]=true;
		for(ll j=2;j<=n;j++){
			if(!vis[j]&&dist[v]+edge[v][j]<dist[j]){
				dist[j]=dist[v]+edge[v][j];
			}
		}
	}
}
int main(){
	while(~scanf("%lld%lld",&n,&m),n,m){
		for(ll i=1;i<=n;i++){
			for(ll j=1;j<=n;j++){
				if(i!=j) edge[i][j]=INF;
				else edge[i][j]=0;
			}
		}
		while(m--){
			ll u,v,w;
			scanf("%lld%lld%lld",&u,&v,&w);
			edge[u][v]=min(w,edge[u][v]);
			edge[v][u]=min(w,edge[v][u]);
		}
		dijstra(1);
		printf("%lld\n",dist[n]);
	}
	
}

顺便贴一个单源最短路二重权值的问题:
传送门:HDU 3790

简单思路:二重权值和一重其实本质上区别不大,这个题是先看最短路,再看花费时间,我们可以先求最短路,再取里面时间最小的。

AC代码

//注意有没有重边问题

#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<list>
#include<cmath>
#include<algorithm>
#define lson node<<1,st,mid
#define rson node<<1|1,mid+1,ed
#define mem(a,x) memset(a,x,sizeof(a))
#define me(a) memset(a,0,sizeof(a))
#define IOS ios::sync_with_stdio(false)
#define lowbit(x) x&(-x)
#define up(i,x,y) for(long long i=x;i<y;i++)
#define down(i,x,y) for(long long i=x;i>=y;i--)
#define in freopen("in.txt","r",stdin)
#define out freopen("out.txt","w",stdout)  	
typedef long long ll;
const ll mod=1e9+7;
const ll INF=0x3f3f3f3f;
const int maxn=10005;
const double pi = acos(-1.0);
using namespace std; 
bool vis[maxn];
struct code{
	int w;
	int t;
};
int n,m;
code edge[maxn][maxn],dist[maxn];
void dijstra(int u){
	mem(vis,false);
	dist[u].w=0;
	dist[u].t=0;
	for(int i=1;i<=n;i++){
		dist[i].w=edge[u][i].w;
		dist[i].t=edge[u][i].t; 
	}
	vis[u]=true;
	for(int i=2;i<=n;i++){
		int min=INF,v=u;
		for(int j=1;j<=n;j++){
			if(!vis[j]&&dist[j].w<min){
				min=dist[j].w;
				v=j;
			}
		}
		vis[v]=1;
		for(int j=1;j<=n;j++){
			if(!vis[j]&&dist[v].w+edge[v][j].w<dist[j].w){
				dist[j].w=dist[v].w+edge[v][j].w;
				dist[j].t=dist[v].t+edge[v][j].t;
			}
			else if(!vis[j]&&dist[j].w==dist[v].w+edge[v][j].w){
				if(dist[j].t>dist[v].t+edge[v][j].t)
				dist[j].t=dist[v].t+edge[v][j].t;
			}
		}
	}
}
int main(){
	while(~scanf("%d%d",&n,&m),n,m){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				if(i!=j) edge[i][j].w=INF,edge[i][j].t=INF;
				else edge[i][j].w=0,edge[i][j].t=0;
			}
		}
		while(m--){
			int u,v,w,t;
			scanf("%d%d%d%d",&u,&v,&w,&t);
			if(edge[u][v].w>w){
				edge[u][v].w=edge[v][u].w=w;
				edge[u][v].t=edge[v][u].t=t;
			}
		} 
		int l,r;
		scanf("%d%d",&l,&r);
		dijstra(l);
		printf("%d %d\n",dist[r].w,dist[r].t);
	}
}


Dijkstra算法的堆优化:

这里采用了优先队列和链式前向星对dijkstra算法进行优化。
优先队列插入操作的时间复杂度是O(logn),所以直接把直接复杂度降下来了,可是存图的时候有个问题,如果我们采用邻接矩阵存图,那么空间复杂度就是O(n^2),对其优化就失去了意义,所以我们采用链式前向星存图。
优先队列的基本操作:点我🤪
链式前向星:点我😀
先贴个题:HDU2544

//堆优化, 复杂度O(nlogn)
#include<iostream>
#include<cstdio> 
#include<cstring>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<list>
#include<cmath>
#include<algorithm>
#include<fstream>
#include<sstream>
#define lson node<<1,st,mid
#define rson node<<1|1,mid+1,ed
#define mem(a,x) memset(a,x,sizeof(a))
#define me(a) memset(a,0,sizeof(a))
#define IOS ios::sync_with_stdio(false)
#define lowbit(x) x&(-x)
typedef long long ll;
const ll mod=1e9+7;
const ll INF=0x3f3f3f3f;
const ll maxn=1e5+5;
const double pi = acos(-1.0);
using namespace std;
ll n,m;
int dist[maxn],head[maxn];
bool vis[maxn];
struct code{
	int to,w,next;
}edge[maxn];
struct node{//优先队列的重载
	int to,w;
	bool operator <(const node &x) const{
		return w>x.w;//最小值优先
		//return w<x.w//最大值优先
	}
};
priority_queue<node>q;
int cnt;
void intn(){//初始化
	cnt=0;
	mem(head,-1);
	mem(vis,false);
	mem(dist,INF);
	
}
void add(int u,int v,int w){//链式前向星加边操作
	cnt++;
	edge[cnt].to=v;
	edge[cnt].w=w;
	edge[cnt].next=head[u];
	head[u]=cnt;
}
void dijkstra(int u){//优化dijkstrea算法
	dist[u]=0;
	q.push(node {u,dist[u]});
	
	while(!q.empty()){
		node x=q.top();
		q.pop();
		if(vis[x.to]) continue;
		vis[x.to]=true;
		for(int i=head[x.to];i!=-1;i=edge[i].next){
			int v=edge[i].to;
			if(dist[v]>dist[x.to]+edge[i].w){
				dist[v]=dist[x.to]+edge[i].w;
				q.push(node {v,dist[v]});
			}
		}
	}
}
int main(){
	while(~scanf("%lld%lld",&n,&m),n,m){
		intn();
		while(m--){
			ll u,v,w;
			scanf("%lld%lld%lld",&u,&v,&w);
			add(u,v,w);
			add(v,u,w);
		}
		dijkstra(1);
		printf("%lld\n",dist[n]);
	}
}
  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值