dsalgo 0.3.10

A package for Datastructures and Algorithms.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105

use crate::strongly_connected_components_topological_sort::toposort;

pub fn scc(g: &[Vec<usize>]) -> Vec<usize> {
    fn labeling(
        g: &[Vec<usize>],
        preorder: &mut Vec<usize>,
        low_order: &mut Vec<usize>,
        order: &mut [usize],
        labels: &mut [usize],
        ord: &mut usize,
        label: &mut usize,
        u: usize,
    ) {
        order[u] = *ord;

        *ord += 1;

        preorder.push(u);

        let n = g.len();

        for &v in g[u].to_owned().iter() {
            if order[v] == n {
                labeling(g, preorder, low_order, order, labels, ord, label, v);

                low_order[u] = low_order[u].min(low_order[v]);
            } else if labels[v] == n {
                low_order[u] = low_order[u].min(order[v]);
            }
        }

        if low_order[u] < order[u] {
            return;
        }

        loop {
            let v = preorder.pop().unwrap();

            labels[v] = *label;

            if v == u {
                break;
            }
        }

        *label += 1;
    }

    let n = g.len();

    let mut preorder = Vec::with_capacity(n);

    let mut low_order = vec![n; n];

    let mut order = vec![n; n];

    let mut labels = vec![n; n];

    let mut ord = 0;

    let mut label = 0;

    for i in 0..n {
        if order[i] == n {
            labeling(
                &g,
                &mut preorder,
                &mut low_order,
                &mut order,
                &mut labels,
                &mut ord,
                &mut label,
                i,
            );
        }
    }

    toposort(labels)
}

#[cfg(test)]

mod tests {

    use super::*;

    #[test]

    fn test() {
        let cases = vec![(
            (6, vec![(1, 4), (5, 2), (3, 0), (5, 5), (4, 1), (0, 3), (4, 2)]),
            vec![3, 1, 2, 3, 1, 0],
        )];

        for ((n, edges), ans) in cases {
            let mut g = vec![vec![]; n];

            for (u, v) in edges {
                g[u].push(v);
            }

            assert_eq!(scc(&g), ans);
        }
    }
}