%I #8 Aug 11 2014 05:09:10
%S 2,3,5,9,17,33,65,35,129,257,513,1025,2049
%N Values of m corresponding to members of A088966.
%C For every k >= 0, 2^k + 1 is in the sequence.
%e a(2) = 3, corresponding to A088966(2) = 8 because A007947(8) = A007947(3+1) and A007947(3) = A007947(8+1).
%p rad:= n -> convert(numtheory:-factorset(n), `*`):
%p count:= 0: lastr:= rad(1):
%p for n from 2 to 10^7 do
%p newr:= rad(n);
%p P[lastr, newr]:= n-1;
%p if assigned(P[newr, lastr]) then
%p count:= count+1; A[count]:= n-1; M[count]:= P[newr, lastr];
%p fi;
%p lastr:= newr;
%p od:
%p seq(M[n], n=1..count); # _Robert Israel_, Aug 11 2014
%K more,nonn,less
%O 1,1
%A _David Wasserman_, Jan 29 2004
%E a(13) from _Robert Israel_, Aug 11 2014