c(1, 2, 3)
x <- c(1, 2, 3)
1:5
rep(1:5, 3)

n <- 5 # pocet parov
ponozky <- rep(1:n, 2)
k <- 7 # pocet vybranych ponoziek
vyber <- sample(ponozky,  # vyberame z vektora ponozky
                size = k, # vyberame k hodnot
                replace = FALSE)  # bez navratu
unique(vyber)
pocetParov <- k - length(unique(vyber)) # length = dlzka vektora

vyberaniePonoziek <- function(n, k){
  ponozky <- rep(1:n, 2)
  vyber <- sample(ponozky,  # vyberame z vektora ponozky
                  size = k, # vyberame k hodnot
                  replace = FALSE)  # bez navratu
  pocetParov <- k - length(unique(vyber)) # length = dlzka vektora
  return(pocetParov)
}

vyberaniePonoziek(5, 7)

set.seed(123) # kvoli reprodukovatelnosti, tie iste nahodne cisla pri novom spusteni kodu
pocty <- replicate(10^4, vyberaniePonoziek(5, 7))
mean(pocty)

n  # 5
k  # 7
k*(k-1)/(2*(2*n-1))

x <- 1:5
cumsum(x)  # kumulativne sucty

priebehOdhadov <- cumsum(pocty)/(1:10^4)  # podiel bude po zlozkach
plot(priebehOdhadov, type = "l") # l = line, typ grafu
abline(h = k*(k-1)/(2*(2*n-1)), col = "blue") # h = horizontalna ciara