estpremc(n):={ //utilise la fonction crible pour tester si n est premier local prem,s; prem:=crible(floor(sqrt(n))); s:=size(prem); if (n==1) return(false); for (k:=0;k<s;k++){ if (irem(n,prem[k])==0){ return(false); } } return(true); }