By Eric Burden | December 24, 2020

In the spirit of the holidays (and programming), I’ll be posting my solutions to the Advent of Code 2020 puzzles here, at least one day after they’re posted (no spoilers!). I’ll be implementing the solutions in R because, well, that’s what I like! What I won’t be doing is posting any of the actual answers, just the reasoning behind them.

Also, as a general convention, whenever the puzzle has downloadable input, I’m
saving it in a file named `input.txt`

.

# Day 23 - Crab Cups

Find the problem description HERE.

## Part One - I’m Not Trapped on This Raft With You

This puzzle presents us with a lovely cup game, administered by a small crab on
a raft. Who’d have ever imagined we end up here? The crab, that’s who. To help
myself out, we’ve got a couple of functions to plot a nice network diagram to
visualize the state. I won’t tell you exactly how many times I used that, mostly
because I’m not sure I can count that high. The ’three vector’ approach to
keeping track of the game state wasn’t actually my first direction, but it was
by far the *fastest*, which you know is going to matter. This is the second (
I think) puzzle where we’ve made good use of a closure. This first part required
a good bit of code, but it could have been accomplished in a less round-about
manner just by following the process outlined in the part one instructions.

```
# Setup ------------------------------------------------------------------------
test_input <- c(389125467)
real_input <- c(198753462)
library(igraph) # For plotting networks, debugging
library(magrittr) # Provides the `%>%` pipe operator
# Functions --------------------------------------------------------------------
# Given a vector of cup labels `cups`, a vector indicating the next cup in
# sequence `next_i`, the index (in `cups`) of the current cup `current_i`,
# the index (in `cups`) of the destination cup `destination_i`, and the indices
# (in `cups`) of any cups that have been picked up `picked_up_i`, plots
# a network diagram of the current state of the cups game
pplot <- function(cups, next_i, current_i, destination_i, picked_up_i) {
g <- graph_from_edgelist(cbind(cups, cups[next_i])) %>%
set_vertex_attr('current', value = F) %>%
set_vertex_attr('current', index = cups[current_i], value = T) %>%
set_vertex_attr('destination', value = F) %>%
set_vertex_attr('destination', index = cups[destination_i], value = T) %>%
set_vertex_attr('picked_up', value = F) %>%
set_vertex_attr('picked_up', index = cups[picked_up_i], value = T)
V(g)$color <- ifelse(V(g)$destination, 'coral', 'lightgreen')
V(g)$color <- ifelse(V(g)$current, 'lightblue', V(g)$color)
V(g)$color <- ifelse(V(g)$picked_up, 'purple', V(g)$color)
plot(g)
}
# Given the closure containing all the elements of the cup game in its
# enclosing environment, extracts those elements from the enclosing environemnt
# and passes them to `pplot()`
pplot_func <- function(f) {
e <- environment(f)
pplot(e$cups, e$next_i, e$current_i, e$destination_i, e$picked_up_i)
}
# Given a number representing cup labels in order, returns a numeric vector
# representing cup labels, in order
parse_input <- function(input) {
input %>%
as.character() %>%
strsplit('') %>%
unlist %>%
as.numeric()
}
# A closure whose enclosing environment will contain the current state of the
# cups game, including:
# - a vector representing the cup labels in the initial order `cups`
# - a vector where the value at each index represents the index of the next
# cup (in `cups`) in series `next_i` (Example, `cups[next_i[3]]` is the next
# cup in sequence after cup `3`)
# - a vector where the value at each index `i` represents the index in `cups`
# containing the value `i`, essentially a reverse lookup for `cups` (Example,
# `val_map[3]` returns the index in `cups` for the value `3`)
# - the maximum value of a cup label `max_cup`
# - a number representing the index in `cups` for the 'current' cup `current_i`
# - a number representing the index in `cups` for the `destination` cup
# `destination_i`
# - a vector representing the indices in `cups` for the cups that have been
# picked up `picked_up_i`
#
# Returns a function that advances the game state according to the rules of
# the game as described in the puzzle instructions
crab_mover <- function(input) {
# Values representing the current game state, see comment for `crab_mover()`
# for explanations of each
cups <- parse_input(input)
next_i <- c(2:length(cups), 1)
val_map <- sort.list(cups)
max_cup <- max(cups)
current_i <- 1
destination_i <- 0
picked_up_i <- 0
function() {
# Pick up three cups. Note the use of the `<<-` assignment operator,
# assigning the values to `picked_up_i` in the enclosing environment
picked_up_i <<- c(
next_i[current_i],
next_i[next_i[current_i]],
next_i[next_i[next_i[current_i]]]
)
# Close the circle, setting the `next_i` for the current cup to the cup
# after the last cup picked up
next_i[current_i] <<- next_i[picked_up_i[3]]
# Identify the 'destination' cup, the next cup in reverse label order from
# the current cup. Can't be one of the cups picked up. If there are no cups
# with labels lower than the current cup, loop back to the top of the cup
# label range (`max_cup`)
destination_value <- cups[current_i]
invalid_values <- c(destination_value, cups[picked_up_i])
while (destination_value %in% invalid_values) {
destination_value <- if (destination_value == 1) {
max_cup
} else {
destination_value - 1
}
}
destination_i <<- val_map[destination_value]
# Add the picked up cups back to the circle. Set the next cup after the
# third picked up cup to the cup following the destination cup. Set the
# next cup after the destination cup to the first picked up cup.
next_i[picked_up_i[3]] <<- next_i[destination_i]
next_i[destination_i] <<- picked_up_i[1]
# Get the next current cup
current_i <<- next_i[current_i]
}
}
# Processing -------------------------------------------------------------------
move_cups <- crab_mover(real_input) # Set up the game
for (i in 1:10) { move_cups() } # Advance 10 times
e <- environment(move_cups) # Extract the game state
# Get cup labels, in order, following the cup labeled '1'
current_cup <- 1
cups_in_order <- character(0)
for (i in 1:8) {
# Get the cup following the cup with the value `current_cup`
next_cup <- e$cups[e$next_i[e$val_map[current_cup]]]
cups_in_order <- c(cups_in_order, next_cup)
current_cup <- next_cup
}
answer1 <- paste(cups_in_order, collapse = '')
```

So, there you have it. I actually implemented it first using a graph and switching the edges about, but that turned out to be way too slow for…

## Part Two - You’re Trapped on This Raft With Me

This crab is a maniac. Part two is the exact same process as part one, just…
way…WAY…**WAY** more. This was the reason for the three vectors, referenced
by index and in place without any copying (I don’t think) or iterating through.
Having got most of the operations down to O(1), we’re able to get the time to
process this down to under a minute. Yay!

```
# Setup ------------------------------------------------------------------------
source('exercise_1.R')
# Functions --------------------------------------------------------------------
# Now add the range from 10 to one million to our cup labels. Yeesh.
parse_input <- function(input) {
input %>%
as.character() %>%
strsplit('') %>%
unlist %>%
as.numeric() %>%
c(10:1000000)
}
# Processing -------------------------------------------------------------------
move_cups <- crab_mover(real_input) # Set up the game
# Progress the game 10 MILLION times, with a progress bar to keep tabs
rounds <- 10000000
pb <- txtProgressBar(min = 0, max = rounds, style = 3)
for (i in 1:rounds) {
if (i %% 100000 == 0) { setTxtProgressBar(pb, i)} # Update `pb` at each %
move_cups()
}
close(pb)
e <- environment(move_cups) # Extract the environment
next_one <- e$cups[e$next_i[e$val_map[1]]] # First cup after '1'
next_two <- e$cups[e$next_i[e$val_map[next_one]]] # Second cup after '1'
answer2 <- next_one * next_two # The answer!
```

I swear to God, if this is that ‘Ferris’ crab trying to make a point about interpreted languages versus compiled languages, I will flip my lid.

## Wrap-Up

So, this is way outside the scope of things I normally think about, but optimizing the heck out of the run time was a really good experience in thinking through how to make the code run as efficiently as I could. There are probably still optimizations to be made here, but the feeling of triumph when that progress bar started moving across my screen at a readily observable rate was a thing a beauty.

If you found a different solution (or spotted a mistake in one of mine), please drop me a line!