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
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
|
use std::fs;
use std::path::Path;
use std::collections::HashSet;
use std::ops::RangeInclusive;
type Instruction = (bool,(isize,isize),(isize,isize),(isize,isize));
/* AOC21 Day 22: https://adventofcode.com/2021/day/22 */
fn main() {
let input = Path::new("resources").join("input.txt");
let content = fs::read_to_string(input).expect("Unable to read input file");
let instructions = parse_input(&content);
println!("Ex1: the result is {}", ex1(&instructions));
println!("Ex2: the result is {}", ex2(&instructions));
}
fn parse_input(input: &str) -> Vec<Instruction> {
input.lines().map(|l| {
let (action, ranges) = l.split_once(" ").unwrap();
let ranges: Vec<_> = ranges.split(',').map(|r| {
let (a,b) = r.split_once('=').unwrap().1.split_once("..").unwrap();
(a.parse().unwrap(),b.parse().unwrap())
}).collect();
(action == "on", ranges[0], ranges[1], ranges[2])
}).collect()
}
trait Clamp {
fn clamp(&self, a: isize, b:isize) -> Self;
}
impl Clamp for RangeInclusive<isize> {
fn clamp(&self, a: isize, b:isize) -> Self {
isize::max(*self.start(),a)..=isize::min(*self.end(),b)
}
}
/// Exercise 1 bruteforces the problem
fn ex1(instructions: &[Instruction]) -> usize {
let mut reactor = HashSet::new();
for instruction in instructions {
let &(act,(x1,x2),(y1,y2),(z1,z2)) = instruction;
for x in (x1..=x2).clamp(-50,50) {
for y in (y1..=y2).clamp(-50,50) {
for z in (z1..=z2).clamp(-50,50) {
if act {
reactor.insert((x,y,z));
} else {
reactor.remove(&(x,y,z));
}
}
}
}
}
reactor.len()
}
type Point = (isize,isize,isize);
#[derive(Debug,Clone)]
struct Cuboid { p1: Point, p2: Point }
impl Cuboid {
fn new(x1: isize, y1: isize, z1: isize, x2: isize, y2: isize, z2: isize) -> Self {
Cuboid { p1: (x1,y1,z1), p2: (x2,y2,z2) }
}
fn len(&self) -> usize {
(self.p2.0 - self.p1.0) as usize *
(self.p2.1 - self.p1.1) as usize *
(self.p2.2 - self.p1.2) as usize
}
fn intersects(&self, other: &Cuboid) -> bool {
let (c1x1,c1y1,c1z1) = self.p1;
let (c1x2,c1y2,c1z2) = self.p2;
let (c2x1,c2y1,c2z1) = other.p1;
let (c2x2,c2y2,c2z2) = other.p2;
c1x1 < c2x2 && c2x1 < c1x2 &&
c1y1 < c2y2 && c2y1 < c1y2 &&
c1z1 < c2z2 && c2z1 < c1z2
}
fn minus(&self, other: &Cuboid) -> Vec<Cuboid> {
let (c1x1,c1y1,c1z1) = self.p1;
let (c1x2,c1y2,c1z2) = self.p2;
let (c2x1,c2y1,c2z1) = other.p1;
let (c2x2,c2y2,c2z2) = other.p2;
let mut xs = vec![c1x1,c1x2,c2x1,c2x2]; xs.sort();
let mut ys = vec![c1y1,c1y2,c2y1,c2y2]; ys.sort();
let mut zs = vec![c1z1,c1z2,c2z1,c2z2]; zs.sort();
let mut diff: Vec<Cuboid> = vec![];
xs.windows(2).for_each(|xs|
ys.windows(2).for_each(|ys|
zs.windows(2).for_each(|zs| {
let cuboid = Cuboid::new(xs[0],ys[0],zs[0],xs[1],ys[1],zs[1]);
if self.intersects(&cuboid) && !other.intersects(&cuboid) {
diff.push(cuboid)
}
})
)
);
diff
}
}
/// Exercise 2 uses a more clever approach keeping track of cuboids and splitting
/// them only when necessary.
fn ex2(instructions: &[Instruction]) -> usize {
let mut reactor = vec![];
for instruction in instructions {
let &(add,(x1,x2),(y1,y2),(z1,z2)) = instruction;
let cuboid1 = Cuboid::new(x1,y1,z1,x2+1,y2+1,z2+1);
reactor = reactor.into_iter().flat_map(|cuboid2|
if cuboid1.intersects(&cuboid2) {
cuboid2.minus(&cuboid1)
} else {
vec![cuboid2]
}
).collect();
if add { reactor.push(cuboid1) }
}
reactor.iter().map(|c| c.len()).sum()
}
#[cfg(test)]
mod tests {
use super::*;
const TEST1: &str = "on x=10..12,y=10..12,z=10..12
on x=11..13,y=11..13,z=11..13
off x=9..11,y=9..11,z=9..11
on x=10..10,y=10..10,z=10..10";
const TEST2: &str = "on x=-20..26,y=-36..17,z=-47..7
on x=-20..33,y=-21..23,z=-26..28
on x=-22..28,y=-29..23,z=-38..16
on x=-46..7,y=-6..46,z=-50..-1
on x=-49..1,y=-3..46,z=-24..28
on x=2..47,y=-22..22,z=-23..27
on x=-27..23,y=-28..26,z=-21..29
on x=-39..5,y=-6..47,z=-3..44
on x=-30..21,y=-8..43,z=-13..34
on x=-22..26,y=-27..20,z=-29..19
off x=-48..-32,y=26..41,z=-47..-37
on x=-12..35,y=6..50,z=-50..-2
off x=-48..-32,y=-32..-16,z=-15..-5
on x=-18..26,y=-33..15,z=-7..46
off x=-40..-22,y=-38..-28,z=23..41
on x=-16..35,y=-41..10,z=-47..6
off x=-32..-23,y=11..30,z=-14..3
on x=-49..-5,y=-3..45,z=-29..18
off x=18..30,y=-20..-8,z=-3..13
on x=-41..9,y=-7..43,z=-33..15";
const TEST3: &str = "on x=-5..47,y=-31..22,z=-19..33
on x=-44..5,y=-27..21,z=-14..35
on x=-49..-1,y=-11..42,z=-10..38
on x=-20..34,y=-40..6,z=-44..1
off x=26..39,y=40..50,z=-2..11
on x=-41..5,y=-41..6,z=-36..8
off x=-43..-33,y=-45..-28,z=7..25
on x=-33..15,y=-32..19,z=-34..11
off x=35..47,y=-46..-34,z=-11..5
on x=-14..36,y=-6..44,z=-16..29
on x=-57795..-6158,y=29564..72030,z=20435..90618
on x=36731..105352,y=-21140..28532,z=16094..90401
on x=30999..107136,y=-53464..15513,z=8553..71215
on x=13528..83982,y=-99403..-27377,z=-24141..23996
on x=-72682..-12347,y=18159..111354,z=7391..80950
on x=-1060..80757,y=-65301..-20884,z=-103788..-16709
on x=-83015..-9461,y=-72160..-8347,z=-81239..-26856
on x=-52752..22273,y=-49450..9096,z=54442..119054
on x=-29982..40483,y=-108474..-28371,z=-24328..38471
on x=-4958..62750,y=40422..118853,z=-7672..65583
on x=55694..108686,y=-43367..46958,z=-26781..48729
on x=-98497..-18186,y=-63569..3412,z=1232..88485
on x=-726..56291,y=-62629..13224,z=18033..85226
on x=-110886..-34664,y=-81338..-8658,z=8914..63723
on x=-55829..24974,y=-16897..54165,z=-121762..-28058
on x=-65152..-11147,y=22489..91432,z=-58782..1780
on x=-120100..-32970,y=-46592..27473,z=-11695..61039
on x=-18631..37533,y=-124565..-50804,z=-35667..28308
on x=-57817..18248,y=49321..117703,z=5745..55881
on x=14781..98692,y=-1341..70827,z=15753..70151
on x=-34419..55919,y=-19626..40991,z=39015..114138
on x=-60785..11593,y=-56135..2999,z=-95368..-26915
on x=-32178..58085,y=17647..101866,z=-91405..-8878
on x=-53655..12091,y=50097..105568,z=-75335..-4862
on x=-111166..-40997,y=-71714..2688,z=5609..50954
on x=-16602..70118,y=-98693..-44401,z=5197..76897
on x=16383..101554,y=4615..83635,z=-44907..18747
off x=-95822..-15171,y=-19987..48940,z=10804..104439
on x=-89813..-14614,y=16069..88491,z=-3297..45228
on x=41075..99376,y=-20427..49978,z=-52012..13762
on x=-21330..50085,y=-17944..62733,z=-112280..-30197
on x=-16478..35915,y=36008..118594,z=-7885..47086
off x=-98156..-27851,y=-49952..43171,z=-99005..-8456
off x=2032..69770,y=-71013..4824,z=7471..94418
on x=43670..120875,y=-42068..12382,z=-24787..38892
off x=37514..111226,y=-45862..25743,z=-16714..54663
off x=25699..97951,y=-30668..59918,z=-15349..69697
off x=-44271..17935,y=-9516..60759,z=49131..112598
on x=-61695..-5813,y=40978..94975,z=8655..80240
off x=-101086..-9439,y=-7088..67543,z=33935..83858
off x=18020..114017,y=-48931..32606,z=21474..89843
off x=-77139..10506,y=-89994..-18797,z=-80..59318
off x=8476..79288,y=-75520..11602,z=-96624..-24783
on x=-47488..-1262,y=24338..100707,z=16292..72967
off x=-84341..13987,y=2429..92914,z=-90671..-1318
off x=-37810..49457,y=-71013..-7894,z=-105357..-13188
off x=-27365..46395,y=31009..98017,z=15428..76570
off x=-70369..-16548,y=22648..78696,z=-1892..86821
on x=-53470..21291,y=-120233..-33476,z=-44150..38147
off x=-93533..-4276,y=-16170..68771,z=-104985..-24507";
#[test]
fn test1() {
let instructions = parse_input(TEST1);
assert_eq!(39, ex2(&instructions));
}
#[test]
fn test2() {
let instructions = parse_input(TEST2);
assert_eq!(590784, ex2(&instructions));
}
#[test]
fn test3() {
let instructions = parse_input(TEST3);
assert_eq!(2758514936282235, ex2(&instructions));
}
}
|
