velicine = c(32, 21, 24, 81, 29, 22, 62, 3, 56, 65, 100, 87,
             68, 83, 77, 51, 76, 67, 38, 9)
N = length(velicine)
sum(velicine)
n = 10
(K = max(velicine))

# u svakom koraku biramo uredjeni par (i, R)
# ako R <= Mi - biramo i-tu jedinku u uzorak
# ako R > Mi - nista, nastavljamo postupak
# uzorak sp - dok ne izaberemo n
# uzorak bp - dok ne izaberemo n razlicitih
uzorak_sp = c()
while(length(uzorak_sp) < n) {
  i = sample(1:N, 1) 
  R = sample(1:K, 1)
  if(R <= velicine[i]) {
    uzorak_sp = c(uzorak_sp, i)
  }
}
uzorak_sp

uzorak_bp = c()
while(length(uzorak_bp) < n) {
  i = sample(1:N, 1)
  R = sample(1:K, 1)
  if(R <= velicine[i] & !(i %in% uzorak_bp)) {
    uzorak_bp = c(uzorak_bp, i)
  }
}
uzorak_bp