read DET: RproofP(binomial(2*p+2*m+2*n,p+m+n),m,n,N,20,R,p,1,10); time(); quit: