最短路问题 - Dijkstra,Bellman_ford,SPFA与SPFA判断负环,Floyd

Dijkstra

朴素Dijkstra
图中可能存在重边和自环,所有边权均为非负值

#include<bits/stdc++.h>

//朴素版Dijkstra适合稀疏图
//稠密图:邻接矩阵

using namespace std;

const int N=510;
int n,m;
int g[N][N];
int dist[N];
bool st[N];//初始为 false 

int dijkstra()
{
	memset(dist,0x3f3f3f3f,sizeof dist);
	dist[1]=0;
	for(int i=1;i<=n;i++)
	{
		int t=-1;
		for(int j=1;j<=n;j++)
			if(!st[j]&&(t==-1||dist[t]>dist[j])) t=j;//确定中间节点 
			
		st[t]=true;
		
		for(int j=1;j<=n;j++)
			dist[j]=min(dist[j],dist[t]+g[t][j]);
	}
	if(dist[n]==0x3f3f3f3f) return -1;
	return dist[n];
}

int main()
{
	scanf("%d%d",&n,&m);
	
	memset(g,0x3f,sizeof g);
	
	while(m--)//记录最短边 
	{
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		g[a][b]=min(g[a][b],c);
	}
	int t=dijkstra();
	printf("%d",t);
	return 0;
}

堆优化版Dijkstra
图中可能存在重边和自环,所有边权均为非负值

#include<bits/stdc++.h>

//堆优化Dijkstra适合稀疏图
//稀疏图:邻接表 
 
using namespace std;
typedef pair<int,int> PII;
const int N=1000010;

int n,m;
//邻接表 
int h[N],w[N],e[N],ne[N],idx;
//vector<vector<PII>> gra;  邻接矩阵


int dist[N];
bool st[N];//初始为 false 

void add(int a,int b,int c)
{
	e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

int dijkstra()
{
	memset(dist,0x3f3f3f3f,sizeof dist);
	dist[1]=0;
	
	priority_queue<PII,vector<PII>,greater<PII>> heap;//小根堆 
	heap.push({0,1});
	//first :距离
	//second:点的位置
	
	while(heap.size())
	{
		auto t=heap.top();
		heap.pop();
		
		int ver=t.second,distance=t.first;
		if(st[ver]) continue;
		st[ver]=true;
		
		for(int i=h[ver];i!=-1;i=ne[i])
		{
			int j=e[i];//此时的 j 为 ver 的下一个点 
			if(dist[j]>w[i]+distance)
			{
				dist[j]=w[i]+distance;
				heap.push({dist[j],j});
			}
		}
	}
	if(dist[n]==0x3f3f3f3f) return -1;
	return dist[n];
}

int main()
{
	scanf("%d%d",&n,&m);
	
	memset(h,-1,sizeof h);
	
	while(m--)//记录最短边 
	{
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		add(a,b,c);
	}
	int t=dijkstra();
	printf("%d",t);
	return 0;
}
#include<bits/stdc++.h>

using namespace std;
const int N=100010,M=2*N;
int n;
int h[N],e[M],ne[M],idx;


void add(int a,int b)
{
	e[idx]=b,ne[idx]=h[a],h[a]=idx++;
} 

int main()
{
    memset(h,-1,sizeof h);
    cin>>n;
    for(int i=1;i<=n-1;i++)
    {
        int a,b;cin>>a>>b;
        add(a,b);
        //add(b,a);
    }
    
    for(int u=1;u<=n;u++)
    {
        cout<<u<<endl;
        for(int i=h[u];i!=-1;i=ne[i])
	    {
	    	int j=e[i];
	    	cout<<i<<" "<<j<<" "<<endl;
 	    }
 	    cout<<endl;
    }
    
    return 0;
}
/*
数据:
9
1 2
1 7
1 4
2 8
2 5
4 3
3 9
4 6

*/

在这里插入图片描述

Bellman_ford
有边数限制的的最短路
图中可能存在重边和自环, 边权可能为负数

#include<bits/stdc++.h>

using namespace std;
const int N=510,M=10010;

int n,m,k;
int dist[N],backup[N];

struct Edge
{
    int a,b,w;
}edges[M];

int bellman_ford()
{
    memset(dist,0x3f3f3f3f,sizeof dist);
    dist[1]=0;
    
    for(int i=1;i<=k;i++)
    {
        memcpy(backup,dist,sizeof dist);
        for(int j=1;j<=m;j++)
        {
            int a=edges[j].a,b=edges[j].b,w=edges[j].w;
            dist[b]=min(dist[b],backup[a]+w);
        }
    }
    if(dist[n]>0x3f3f3f3f/2) return -1e8;
    //返回一个不可能出现的值,特判一下
    return dist[n];
}

int main()
{
    scanf("%d%d%d",&n,&m,&k);
    
    for(int i=1;i<=m;i++)
    {
        int a,b,w;
        scanf("%d%d%d",&a,&b,&w);
        edges[i]={a,b,w};
    }
    if(bellman_ford()==-1e8) puts("impossible");
    else printf("%d",bellman_ford());
    return 0;
}

SPFA
SPFA是在Bellman_ford上做的一点优化

(对于有些题spfa可能会被卡掉,谨慎使用,如有需要请自行百度)
Bellman_ford算法可以存在负权回路,
因为其存在循环次数限制,最终也不会进入死循环;
但是SPFA算法不可以,
由于用了队列来存储,只要发生了更新就会不断的入队,从而进入死循环,
因此存在负权回路,不要使用SPFA。

图中可能存在重边和自环, 边权可能为负数
使用spfa时,需要保证不存在负环
负环是路径之和是负数的环,不是指包含负权边的环。

#include<bits/stdc++.h>

using namespace std;
const int N=100010;
int h[N],e[N],ne[N],w[N],idx;
int dist[N];
bool st[N];
int n,m;

void add(int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}

int spfa()
{
    memset(dist,0x3f,sizeof dist);
    dist[1]=0;
    
    queue<int> q;
    q.push(1);
    st[1]=true;
    
    while(q.size())
    {
        int t=q.front();
        q.pop();
        //将t弹出来了,重置t为false
        st[t]=false;
        
        for(int i=h[t];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                if(!st[j])
                {
                    q.push(j);
                    st[j]=true;
                }
            }
        }
    }
    //为什么在Bellman_ford中特判的是 dist[n]>0x3f3f3f3f/2
    //而在spfa中特判的是 dist[n]==0x3f3f3f3f;
    /*
	  最主要的原因就是,
	  Bellman_ford中对于每一个点都进行的遍历,
	  spfa只遍历了需要的点
	*/
    if(dist[n]==0x3f3f3f3f) return -1e9;
    return dist[n];
}

int main()
{
    memset(h,-1,sizeof h);
    cin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        int a,b,c;cin>>a>>b>>c;
        add(a,b,c);
    }
    if(spfa()==-1e9) cout<<"impossible";
    else cout<<spfa();
    return 0;
}

SPFA判断负环

//负环是路径之和是负数的环,不是指包含负权边的环。
#include<bits/stdc++.h>

using namespace std;
const int N=10010;
int n,m;
int h[N],e[N],ne[N],w[N],idx;
int dist[N],cnt[N];
bool st[N];

void add(int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}

bool spfa()
{
	//只要负环存在必定存在某些点的距离最终被更新为负无穷,因此我们不需要对距离进行初始化
    queue<int> q;
    for(int i=1;i<=n;i++)
    {
        st[i]=true;
        q.push(i);
    }
    
    while(q.size())
    {
        int t=q.front();
        q.pop();
        
        st[t]=false;
        
        for(int i=h[t];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                cnt[j]=cnt[t]+1;
                if(cnt[j]>=n) return true;
                if(!st[j])
                {
                    q.push(j);
                    st[j]=true;
                }
            }
        }
        
    }
    return false;
}

int main()
{
    cin>>n>>m;
    memset(h,-1,sizeof h);
    for(int i=1;i<=m;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    if(spfa()) cout<<"Yes";
    else cout<<"No";
    return 0;
}

Floyd
基于动态规划的思想,求多源最短路
可以存在负权边,不能存在负权回路

#include<bits/stdc++.h>

using namespace std;
const int N=210,INF=1e9;
int d[N][N];//邻接矩阵
int n,m,Q;

void floyd()//n^3  8000000次 时间复杂度高
{
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
}

int main()
{
    scanf("%d%d%d",&n,&m,&Q);
    //初始化
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(i!=j) d[i][j]=INF;
    while(m--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        d[a][b]=min(d[a][b],c);//初始化并更新从a到b的最短距离
    }
    floyd();
    while(Q--)
    {
        int a,b;scanf("%d%d",&a,&b);
        if(d[a][b]>INF/2) puts("impossible");
        //d[a][b] 为什么>INF/2,见下方代码与图
        else printf("%d\n",d[a][b]);
    }
    return 0;
}

错误代码:

#include<bits/stdc++.h>

using namespace std;
const int N=210,INF=1e9;
int d[N][N];
int n,m,k;

void floyd()
{
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
}

int main()
{
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(i!=j) d[i][j]=INF;
    while(m--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        d[a][b]=min(d[a][b],c);
    }
    
    floyd();
    
    while(k--)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        if(d[x][y]==INF) puts("impossible");
        //
        else printf("%d\n",d[x][y]);
    }
    
    return 0;
}

在这里插入图片描述

//1~4 存在   4~1 不存在
输入
5 10 10
1 5 10
3 4 9
3 2 10
3 3 7
3 1 5
4 5 7
1 2 -5
5 1 6
4 3 8
5 5 6
4 1
1 4
5 2
3 2
5 4
4 5
2 2
3 5
3 1
1 1

输出
13
999999995
1
0
impossible
7
0
15
5
0

标准答案
13
impossible
1
0
impossible
7
0
15
5
0
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值