https://twitter.com/fermatslibrary/status/1449027186267197463

this-is-pretty-cool.R

library(tidyverse)
library(stevemisc)

jenny()

num_nums <- tibble()

for (i in 1:100) {
  tibble(x = runif(50),
         runsum = cumsum(x),
         lag = lag(runsum),
         lag1 = ifelse(is.na(lag), x, lag)) %>%
    mutate(seq = 1:n()) %>%
    filter(runsum > 1 & lag1 < 1) -> aa
  
  bind_rows(num_nums, aa) -> num_nums
}

mean(num_nums$seq)