# 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
# 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)
IyBUd28gYm9udXMKYm9udXNfQSA8LSBjKDI1MDAsIDIwMDAsIDE1MDAsIAogICAgICAgICAgICAgcmVwKDI1MCwgMTApLCByZXAoMTUwLCAxMCkpCmJvbnVzX0IgPC0gYygzNTAwLCAyMDAwLCA2MDAsIAogICAgICAgICAgICAgcmVwKDI1MCwgMTApLCByZXAoNzAsIDEwKSkKCiMjIFNpbXVsYXRlIDEwMDAwMCBiZXRzIC0tLS0Kc2V0LnNlZWQoMjAyMDA3MjIpCiMgVGltZXMgb2Ygc2ltdWxhdGlvbgpOc2ltIDwtIDFlKzUKIyAxMDAwMDAgb3V0Y29tZQpPUF9BIDwtIE9QX0IgPC0gdmVjdG9yKGxlbmd0aCA9IE5zaW0pCiMgcG9vbApwb29sIDwtIHNwcmludGYoIiUwNC5mIiwwOjk5OTkpCiMgcHJpbmNpcGFsCnByaW5jaXBhbCA8LSAxCmZvcihuIGluIDE6TnNpbSl7CiAgIyBjYXQobiwgIi8iLE5zaW0pCiAgYmV0X251bSA8LSBzYW1wbGUocG9vbCwgMSkKICBMb3R0b19yZXN1bHQgPC0gc2FtcGxlKHBvb2wsIDIzLCByZXBsYWNlID0gVCkKICBPUF9BW25dIDwtIHN1bSgoTG90dG9fcmVzdWx0PT1iZXRfbnVtKSpib251c19BKSAtIHByaW5jaXBhbAogIE9QX0Jbbl0gPC0gc3VtKChMb3R0b19yZXN1bHQ9PWJldF9udW0pKmJvbnVzX0IpIC0gcHJpbmNpcGFsCiAgIyBjYXQoIlwwMTQiKQp9CgojIyBPdXRwdXQgLS0tLQojIGNvdW50cwp0YWJsZShPUF9BKQp0YWJsZShPUF9CKQ==