x <- 1:5
dis <- c(0, 0, 0)
n <- 1000000
for(i in 1:n){
    samps <- sample(x, 3)
    m <- max(samps)
    dis[m - 2] <- dis[m - 2] + 1
}
dis <- dis / n


蒙特卡洛模拟1000000次