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