This function computes all sample vectors of size size composed of the elements 1, ..., n, either with or without replacement
of elements. If size = n and replace = TRUE then the list of all sample vectors corresponds to a list of all permutations of
the integers 1, ..., n.