# Two bonus
bonus_A <- c(2500, 2000, 1500, 
             rep(250, 10), rep(150, 10))
bonus_B <- c(3500, 2000, 600, 
             rep(250, 10), rep(70, 10))

## Simulate 100000 bets ----
set.seed(20200722)
# Times of simulation
Nsim <- 1e+5
# 100000 outcome
OP_A <- OP_B <- vector(length = Nsim)
# pool
pool <- sprintf("%04.f",0:9999)
# principal
principal <- 1
for(n in 1:Nsim){
  # cat(n, "/",Nsim)
  bet_num <- sample(pool, 1)
  Lotto_result <- sample(pool, 23, replace = T)
  OP_A[n] <- sum((Lotto_result==bet_num)*bonus_A) - principal
  OP_B[n] <- sum((Lotto_result==bet_num)*bonus_B) - principal
  # cat("\014")
}

## Output ----
# counts
table(OP_A)
table(OP_B)