Repunit primes: Difference between revisions

=={{header|PARI/GP}}==

{{trans|Julia}}

<syntaxhighlight lang="PARI/GP">

default(parisizemax, "128M")

default(parisize, "64M");

repunitprimeinbase(n, base) = {

repunit = sum(i=0, n-1, base^i); /* Construct the repunit */

return(isprime(repunit)); /* Check if it's prime */

}

{

for(b=2, 40,

print("Base ", b, ": ");

for(n=1, 2700,

if(repunitprimeinbase(n, b), print1(n, " "))

);

print(""); /* Print a newline */

)

}

</syntaxhighlight>