编辑代码

set.seed(10415)
Card<-c(0:51)
n=100000
result<-matrix(nrow=n,ncol=5)
for(i in 1:n){
    result[i,]<-sample(x=Card,size=5)
}
result_rem<-result%%13
for(i in 1:n){
    result_rem[i,]<-sort(result_rem[i,])
}
result_rem1<-result_rem-apply(result_rem,1,min)
n2=0
for(i in 1:n){
    if(all(result_rem1[i,]==c(0:4)))
    n2=n2+1
}
n2/n
(9*4^5)/choose(52,5)