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