OFFSET
1,1
COMMENTS
For every k >= 0, 2^k + 1 is in the sequence.
EXAMPLE
MAPLE
rad:= n -> convert(numtheory:-factorset(n), `*`):
count:= 0: lastr:= rad(1):
for n from 2 to 10^7 do
newr:= rad(n);
P[lastr, newr]:= n-1;
if assigned(P[newr, lastr]) then
count:= count+1; A[count]:= n-1; M[count]:= P[newr, lastr];
fi;
lastr:= newr;
od:
seq(M[n], n=1..count); # Robert Israel, Aug 11 2014
CROSSREFS
KEYWORD
more,nonn,less
AUTHOR
David Wasserman, Jan 29 2004
EXTENSIONS
a(13) from Robert Israel, Aug 11 2014
STATUS
approved