-
-
Notifications
You must be signed in to change notification settings - Fork 2.4k
/
Copy pathburrows_wheeler_transform.rs
117 lines (106 loc) · 3.21 KB
/
burrows_wheeler_transform.rs
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
106
107
108
109
110
111
112
113
114
115
116
117
pub fn burrows_wheeler_transform(input: &str) -> (String, usize) {
let len = input.len();
let mut table = Vec::<String>::with_capacity(len);
for i in 0..len {
table.push(input[i..].to_owned() + &input[..i]);
}
table.sort_by_key(|a| a.to_lowercase());
let mut encoded = String::new();
let mut index: usize = 0;
for (i, item) in table.iter().enumerate().take(len) {
encoded.push(item.chars().last().unwrap());
if item.eq(&input) {
index = i;
}
}
(encoded, index)
}
pub fn inv_burrows_wheeler_transform<T: AsRef<str>>(input: (T, usize)) -> String {
let len = input.0.as_ref().len();
let mut table = Vec::<(usize, char)>::with_capacity(len);
for i in 0..len {
table.push((i, input.0.as_ref().chars().nth(i).unwrap()));
}
table.sort_by(|a, b| a.1.cmp(&b.1));
let mut decoded = String::new();
let mut idx = input.1;
for _ in 0..len {
decoded.push(table[idx].1);
idx = table[idx].0;
}
decoded
}
#[cfg(test)]
mod tests {
use super::*;
#[test]
//Ensure function stand-alone legitimacy
fn stand_alone_function() {
assert_eq!(
burrows_wheeler_transform("CARROT"),
("CTRRAO".to_owned(), 1usize)
);
assert_eq!(inv_burrows_wheeler_transform(("CTRRAO", 1usize)), "CARROT");
assert_eq!(
burrows_wheeler_transform("THEALGORITHMS"),
("EHLTTRAHGOMSI".to_owned(), 11usize)
);
assert_eq!(
inv_burrows_wheeler_transform(("EHLTTRAHGOMSI".to_string(), 11usize)),
"THEALGORITHMS"
);
assert_eq!(
burrows_wheeler_transform("!.!.!??.=::"),
(":..!!?:=.?!".to_owned(), 0usize)
);
assert_eq!(
inv_burrows_wheeler_transform((":..!!?:=.?!", 0usize)),
"!.!.!??.=::"
);
}
#[test]
fn basic_characters() {
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("CARROT")),
"CARROT"
);
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("TOMATO")),
"TOMATO"
);
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("THISISATEST")),
"THISISATEST"
);
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("THEALGORITHMS")),
"THEALGORITHMS"
);
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("RUST")),
"RUST"
);
}
#[test]
fn special_characters() {
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("!.!.!??.=::")),
"!.!.!??.=::"
);
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("!{}{}(((&&%%!??.=::")),
"!{}{}(((&&%%!??.=::"
);
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("//&$[]")),
"//&$[]"
);
}
#[test]
fn empty() {
assert_eq!(
inv_burrows_wheeler_transform(burrows_wheeler_transform("")),
""
);
}
}