distincts=0; N=4000000; evtA=0 V<-sample(1:900,N,replace=TRUE); for (v in V) { ok=FALSE; #cat('v=',v,'\n'); if (v%%2==0){ #cat('2|v',v,evtA,'\n'); if (v%%4!=0){ #cat('4!|v',v,v%%4,'\n') if (v%%6!=0) { ok=TRUE; evtA<-evtA+1;} } } } print(evtA/N);