spfa 判负环
基于抽屉原理,如果一条正在搜索的最短路径上的点的个数大于总共点的个数,则说明路径上一定有至少重复的两个点,走了回头路。又因为满足最短路径的性质,所以一定存在负环。
可以处理负权边的最短路算法就可以找环。通常用spfa,记录每个点的cnt,当cnt >= 点的个数说明有负环。
一个不稳定优化:仅仅找负环时,将spfa用的队列换成栈可能会加速。因为用栈的spfa搜索类似于深搜,在环上更容易一次把环跑完,所以更容易搜到负环。
另,找负环不需要初始化dist数组,因为只要有负环,某个点的距离就能被一直更新下降,cnt总会大于点的个数。
再另,如果建立了虚拟源点连接一条0权边到所有的点,则总共的点数n要 + 1。但cnt的判断要看初始化的情况。这里一次性将所有点入队后仍然只需要cnt >= n,因为虚拟源点到到j的边数 >= n + 1,一次性入队已经扣除一条边了。但如果没有一次性全部入队则需要 cnt >= n + 1。(切记)
bool spfa(){
int hh = 0, tt = 0;
memset(st, 0, sizeof st);
memset(cnt, 0, sizeof cnt);
memset(dist, 0, sizeof dist);// 可以不初始化dist
for(int i = 1; i <= n; ++ i){// 虚拟源点连的边全部入队
// 全部点入队是因为负环不一定带起点
q[tt ++] = i;
st[i] = true;
}
while(hh != tt){
int t = q[hh ++];
if(hh == N) hh = 0;
st[t] = false;
for(int i = h[t]; ~ i; 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; // 虚拟源点到到j的边数 >= n + 1(已经算过一次边了),则说明跑了n + 2个点,有环
if(!st[j]){
q[tt ++] = j;
if(tt == N) tt = 0;
st[j] = true;
}
}
}
}
return false;
}