RosettaCodeData/Task/Repunit-primes/PARI-GP/repunit-primes.parigp
2024-03-06 22:25:12 -08:00

18 lines
386 B
Text

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 */
)
}