OFFSET
1,1
COMMENTS
EXAMPLE
a(9) = 89 is a term because 2^11 - 1 = 23 * 89 whose greatest prime factor is 89.
MAPLE
A:= NULL: count:= 0: P:= {}: p:= 2:
while count < 40 do
p:= nextprime(p);
m:= NumberTheory:-MultiplicativeOrder(2, p);
t:= 2^m-1;
t:= t/p^padic:-ordp(t, p);
for q in P do if t mod q = 0 then
t:= t/q^padic:-ordp(t, q);
if t = 1 then break fi;
fi od;
if t = 1 then count:= count+1; A:= A, p fi;
P:= P union {p};
od:
A;
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, May 11 2026
STATUS
approved
