N = 15
Mi = c(23,30,41,26,53,60,28,52,113,72,80,35,42,38,52)
sum(Mi)

n = 5
K = max(Mi)
K

# 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 njih 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 <= Mi[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)
  # ako R<=Mi i i-ta jedinka nije do sada izabrana u uzorak, dodajemo je u uzorak
  if(R <= Mi[i] && !(i %in% uzorak_bp)) {
    uzorak_bp = c(uzorak_bp, i)
  }
}
uzorak_bp