RosettaCodeData/Task/Repunit-primes/Maxima/repunit-primes.maxima
2026-04-30 12:34:36 -04:00

10 lines
307 B
Text

repunit(n, base) := (base^n-1)/(base-1)$
repunit_primes(base) := sublist_indices(makelist(repunit(i, base), i, 1, 1000), primep)$
for i from 2 thru 16 do(
sprint("Base", i),
newline(),
rup: repunit_primes(i),
if length(rup) > 0 then apply(sprint, rup) else sprint("None"),
newline()
)$