Primality by trial division: Difference between revisions

m
No edit summary
Line 675:
 
=={{header|PARI/GP}}==
<lang parigp>trial(n)={
if(n%2 ==< 04, return(n==2 > 1)); /* Handle negatives */
ifforprime(n < 9p=2, returnsqrt(n > 1));,
forprime(p=3,sqrt(n),
if(n%p == 0, return(0))
);