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
pub mod bellman_ford {
pub fn shortest_path(
graph: &Vec<Vec<(usize, i64)>>,
start: usize,
inf: i64,
) -> (Vec<i64>, Vec<bool>) {
let n = graph.len();
let mut dist = vec![inf; n];
dist[start] = 0;
for _ in 0..n {
for v in 0..n {
for &(to, cost) in &graph[v] {
if dist[v] == inf || dist[to] <= dist[v] + cost {
continue;
}
dist[to] = dist[v] + cost;
}
}
}
let mut negative = vec![false; n];
for _ in 0..n {
for v in 0..n {
for &(to, cost) in &graph[v] {
if dist[v] == inf {
continue;
}
if dist[to] > dist[v] + cost {
dist[to] = dist[v] + cost;
negative[to] = true;
}
if negative[v] {
negative[to] = true;
}
}
}
}
return (dist, negative);
}
}
#[cfg(test)]
mod tests {
use super::*;
use std;
use utils::test_helper::TestCaseProducer;
#[test]
fn solve_grl_1_b() {
let mut input = TestCaseProducer::new_from_directory("./assets/GRL_1_B/in/");
let mut output = TestCaseProducer::new_from_directory("./assets/GRL_1_B/out/");
while !input.is_empty() {
let v: usize = input.next();
let e: usize = input.next();
let r: usize = input.next();
let mut graph = vec![vec![]; v];
for _ in 0..e {
let s: usize = input.next();
let t: usize = input.next();
let d: i64 = input.next();
graph[s].push((t, d));
}
let inf = std::i64::MAX;
let (dist, negative) = bellman_ford::shortest_path(&graph, r, inf);
let mut neg = false;
for &b in &negative {
neg = neg || b;
}
if neg {
let out1: String = output.next();
let out2: String = output.next();
assert_eq!(out1, "NEGATIVE");
assert_eq!(out2, "CYCLE");
} else {
for i in 0..v {
if dist[i] == inf {
let out: String = output.next();
assert_eq!(out, "INF");
} else {
let out: i64 = output.next();
assert_eq!(dist[i], out);
}
}
}
}
}
}