login
Values of m corresponding to members of A088966.
1

%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