Simulação do problema da caverna

testes = 10000

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.0414

hist(dias)

max(dias)

## [1] 31