最短路最常用的两个算法,网上也有很多不同的写法,总结了下,最后参考了qsc B站视频的代码,因为比较精简.
Dijkstra+优先队列 O(V+E)logV:
Ps:我曾经思考过在稀疏图上优化过的版本比普通的Dijkstra上快是肯定的,但是在完全图上,E就是V2,那复杂度岂不是(V2+V)LogV 比V2还大,那优化什么?.不过查了资料说这只是理论上分析。在实际使用中由于他的入队条件往往得不到满足,所以实际的使用效率会大大的好于O(V^2)的版本,所以可放心使用。
参考http://blog.sina.com.cn/s/blog_476a25110100mag0.html
const int maxn = 205; //顶点数自己定
const int INF = 0x3f3f3f3f;
vector<pair<int, int> >E[maxn];//建图
int d[maxn]; // d[i]保存源点到i点的最短距离.
inline void Init() {for(int i = 0; i < maxn; i++) E[i].clear(), dis[i] = INF;}
void dijkstra(int s,int t) { //s源点,t终点
priority_queue<pair<int ,int > >q;
d[s] = 0;
q.push(make_pair(-d[s], s));//优先队列每次取最大,故-dis[s]保证每次取到最小.
while(!q.empty()) {
int now = q.top().second;
if(now == t) return;
q.pop();
for(int i = 0; i < (int)E[now].size(); i++) {
int next = E[now][i].first;
int t = E[now][i].second + d[now];
if(d[next] > t) {
dis[next] = temp;
q.push(make_pair(-d[next], next));
}
}
}
}
SPFA(解决负权) O(kE)E为边数,k是期望,一般为2或3。
const int maxn=205; //顶点数自己定
vector <pair<int,int> >E[maxn];//存图
int dis[maxn];
bool vis[maxn];
void Init()
{
for(int i=0;i<maxn;i++)
E[i].clear(),vis[i]=false,dis[i]=1e9;
}
void SPFA(int s,int t)//源点,终点.
{
queue<int>q;
q.push(s);
dis[s]=0;
vis[s]=true;
while(!q.empty())
{
int now=q.front();
q.pop();
vis[now]=false;
for(int i=0;i<E[now].size();i++)
{
int next=E[now][i].first;
if(dis[next]>dis[now]+E[now][i].second)
{
dis[next]=dis[now]+E[now][i].second;
if(vis[next])continue; //已经在队列里的就跳过
vis[next]=true;
q.push(next);
}
}
}
}
练习:HDU 1874 http://acm.hdu.edu.cn/showproblem.php?pid=1874
dijkstra+优先队列版
#include <bits/stdc++.h>
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define rep(i,n) for(int i=0;i<n;i++)
const int INF = 0x3f3f3f3f;
const int maxn = 205;
using namespace std;
//head
int d[maxn];
vector<pair<int,int> > E[maxn];
inline void init(int n) {
rep(i, n) E[i].clear(), d[i] = INF;
}
void dijkstra(int s, int t) {
d[s] = 0;
priority_queue<pair<int,int> >q;
q.push(mp(-d[s],s));
while(!q.empty()) {
int now = q.top().second;
if(now == t) return;
q.pop();
for(int i = 0; i < (int)E[now].size(); i++) {
int next = E[now][i].fi;
int temp = d[now] + E[now][i].se;
if(d[next] > temp) {
d[next] = temp;
q.push(mp(-d[next],next));
}
}
}
}
int main() {
int n, m;
while(cin >> n >> m) {
init(n);
rep(i, m) {
int u, w, v;
cin >> u >> w >> v;
E[u].pb(mp(w,v));
E[w].pb(mp(u,v));
}
int s,t;
cin >> s >> t;
dijkstra(s,t);
int ans = (d[t] == INF ? -1 : d[t]);
cout << ans << endl;
}
return 0;
}
spfa版
#include <bits/stdc++.h>
using namespace std;
const int maxn=205; //顶点数自己定
const int INF = 0x3f3f3f3f;
vector <pair<int,int> >E[maxn];//存图
int d[maxn];
bool vis[maxn]; //判断是否在队列里
void init() {
for(int i=0;i<maxn;i++) E[i].clear(), vis[i]=false,d[i]=INF;
}
void spfa(int s,int t) {
queue<int>q;
q.push(s);
d[s] = 0;
vis[s] = true;
while(!q.empty()) {
int now = q.front();
q.pop();
vis[now] = false;
for(int i = 0; i < (int)E[now].size(); i++) {
int next = E[now][i].first;
if(d[next] > d[now] + E[now][i].second) {
d[next] = d[now] + E[now][i].second;
if(vis[next]) continue; //已经在队列里的就跳过
vis[next] = true;
q.push(next);
}
}
}
}
int main() {
int n,m;
while(cin >> n >> m) {
init();
for(int i = 0; i < m; i++) {
int x,y,z;
cin>>x>>y>>z;
E[x].push_back(make_pair(y,z));
E[y].push_back(make_pair(x,z));
}
int s,t;
cin >> s >> t;
spfa(s,t);
int ans = (d[t] == INF ? -1 : d[t]);
cout << ans << endl;
}
return 0;
}
效率比较参考: http://blog.sina.com.cn/s/blog_476a25110100mafq.html