Function reach_index 

pub fn reach_index(adj: &[Vec<usize>], sccs: &[Vec<usize>]) -> ReachIndex

Build a ReachIndex for pairwise connectivity queries. Same forward-closure construction as reachability, but retains the reach-sets instead of collapsing them to counts.