注意:==超级容易被卡!==没有负权边千万千万别用,用堆优化的DJ!!!
思路
优化Bellman-Ford
对Bellman-Ford算法 中的dist[b] = min(dist[b], dist[a] + w )
做优化。
因为dist[a]不一定会更新dist[b]。具体而言就是,当dist[a]
变小了, 与其相连的dist[b]才会变小。
用宽搜做优化。
大致模板
*1 定义一个
queue q
, 用来储存dist变小了的点;
*2 while q 不空
*3 取出队头 t,pop(t);
*4 更新 t 的所有出边。
*5 如果更新成功 且 队列中没有b, 将 b 入队;
写法和Dijkstra算法很像,不过是把针对“点”的操作 改为 针对“边”的操作。
注:很多正权图可以用SPFA。
例题
#include<bits/stdc++.h>
using namespace std;
const int N = 100010;
int n, m;
int dist[N];
int ne[N], h[N], e[N], w[N], idx = 0;
bool st[N];
void add(int a, int b, int c){
e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
}
int spfa(){
memset(dist, 0x3f, sizeof dist);
dist[1] = 0;
st[1] = 1;
queue<int> q;
q.push(1);
while(q.size()){
int t = q.front();
q.pop();
st[t] = 0;
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]){
st[j] = 1;
q.push(j);
}
}
}
}
if(dist[n] == 0x3f3f3f3f) return -1;
return dist[n];
}
int main(){
cin >> n >> m;
memset(h, -1, sizeof h);
for(int i = 0; i < m; i ++){
int a, b, c;
cin >> a >> b >> c;
add(a, b, c);
}
int t = spfa();
if(t == -1) cout << "impossible";
else cout << t;
return 0;
}
————————————————————————————————————