read DET: a:=binomial(m+n,n)*binomial(2*p-m-n,p-n): RproofP(a,m,n,N,20,R,p,100,120); time(); quit: