here is one way
set.seed(999)
B<- A <-replicate(
200,
sample(x=1:49,size=6,replace=FALSE)
)
winvector <- c(33, 41,59,02,20,71)
payout <- function(x){
ifelse(x>=3,10^x,0)
}
B[,] <- purrr::map_lgl(A,~ . %in% winvector)
scorevec <- colSums(B) %>% setNames(1:ncol(A))
table(scorevec)
payoutvec <- payout(scorevec)
table(payoutvec)
(sum(payoutvec))