set.seed(10415)
Card<-c(0:51)
n=100000
result<-matrix(nrow=n,ncol=3)
for(i in 1:n){
result[i,]<-sample(x=Card,size=3)
}
result_rem<-result%%13
for(i in 1:n){
result_rem[i,]<-sort(result_rem[i,])
}
result_rem<-result_rem-apply(result_rem,1,min)
result_flo<-floor(result/13)
n3=0
for(i in 1:n){
if(all(result_rem[i,]==c(0:2))&length(unique(result_flo[i,]))==1)
n3=n3+1
}
n3/n
4*11/choose(52,3)