use super::{Flags, Graph, VertexId};
use crate::vector2::*;
use std::cmp::min;
impl Graph {
pub(crate) fn tarjan(&mut self) -> Vector2<VertexId> {
if self.lines.len() <= 1 {
let mut sccs = Vector2::with_capacities(self.lines.len(), self.lines.len());
sccs.push();
sccs.push_to_last(VertexId(0));
return sccs;
}
let mut call_stack = vec![(VertexId(1), 0, true)];
let mut index = 0;
let mut stack = Vec::new();
let mut scc = Vector2::new();
'recursion: while let Some((n_l, i, first_visit)) = call_stack.pop() {
if first_visit {
let l = &mut self[n_l];
l.index = index;
l.lowlink = index;
l.flags = l.flags | Flags::ONSTACK | Flags::VISITED;
stack.push(n_l);
index += 1;
} else {
let &(_, n_child) = self.child(n_l, i);
self[n_l].lowlink = self[n_l].lowlink.min(self[n_child].lowlink);
}
for j in i..self[n_l].n_children {
let &(_, n_child) = self.child(n_l, j);
if !self[n_child].flags.contains(Flags::VISITED) {
call_stack.push((n_l, j, false));
call_stack.push((n_child, 0, true));
continue 'recursion;
} else if self[n_child].flags.contains(Flags::ONSTACK) {
self[n_l].lowlink = min(self[n_l].lowlink, self[n_child].index)
}
}
if self[n_l].index == self[n_l].lowlink {
let n_scc = scc.len();
scc.push();
loop {
match stack.pop() {
None => break,
Some(n_p) => {
self[n_p].scc = n_scc;
self[n_p].flags ^= Flags::ONSTACK;
scc.push_to_last(n_p);
if n_p == n_l {
break;
}
}
}
}
}
}
scc
}
}