nums <- 1:100
probs <- rep(1 / 100, 100)
numList <- NULL
for(i in 1:100) {
    num <- sample(nums, 1, prob = probs)
    numList <- c(numList, num)
    probs[num] <- probs[num] / 2;
    probs[-num] <- probs[-num] + probs[num] / 99
}
numList

(nums <- 1:100)
(probs <- rep(1 / 100, 100))
(numList <- NULL)
for(i in 1:100) {
  num <- sample(nums, 20, prob = probs)
  numList <- c(numList, num)
  probs[num] <- probs[num] / 2
  probs[-num] <- probs[-num] + sum(probs[num])/ 80
  sum(probs[-num])+sum(probs[num])
}
numList

谢谢大神，懂啦~