set.seed(10415)
Card<-c(0:51)
n=100000
result<-matrix(nrow=n,ncol=4)
for(i in 1:n){
    result[i,]<-sample(x=Card,size=4)
}
result_rem<-result%%13  #余数
n1=0
for(i in 1:n){
    if(length(unique(result_rem[i,]))==1)
    n1<-n1+1
}
n1/n
13/choose(52,4)