Simulação do problema da caverna

testes = 1000

dias = numeric(testes)

for (i in 1:testes) {
  soma = 0
  
  repeat {
    resultado = sample(c(0, 1, 2), size = 1, prob = c(1/3, 1/3, 1/3))
    
    if (resultado == 0) {
      break
    }
    
    soma = soma + resultado
  }
  
  dias[i] = soma
}

mean(dias)      

## [1] 3.143

hist(dias)