OFFSET
1,1
COMMENTS
EXAMPLE
a(4) = 13 is a term because 13 is the greatest prime factor of 2^6 + 1 = 65 = 5 * 13.
MAPLE
A:= 2: count:= 1: P:= {}: p:= 2:
while count < 37 do
p:= nextprime(p);
m:= NumberTheory:-MultiplicativeOrder(2, p);
if m::odd then next fi;
m:= m/2;
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 13 2026
EXTENSIONS
a(40)-a(43) from Jinyuan Wang, May 23 2026
STATUS
approved
