A221715 - OEIS

%I #9 Mar 02 2013 20:31:53

%S 1,5,3,9,33,11,13,7,21,25,20,17,23,19,49,39,35,47,27,29,15,43,37,97,

%T 55,41,51,57,85,73,193,111,59,61,31,89,53,329,145,192,221,237,87,105,

%U 81,101,113,481,175,149,69,93,77,99,79,71,67,65,45,75,195,157,141,189,109,83,103,229,121,352,849,233,345,297,137,187,309,281,379,155,199,159

%N Start of n-th doubling run in A114183.

%H N. J. A. Sloane, <a href="/A221715/b221715.txt">Table of n, a(n) for n = 1..7328</a>

%p # A114181

%p M:=10000; M2:=1000;

%p s1:={1}; v0:=[1]; v1:=[1]; v2:=[]; vi:=Array(1..M2);

%p t1:=1; r1:=1; vi[1]:=1;

%p for n from 2 to M do

%p t2:=floor(sqrt(t1));

%p if t2 in s1 then

%p    v0:=[op(v0),2*t1]; s1:={op(s1),2*t1}; r1:=r1+1; t1:=2*t1;

%p    if t1<=M2 then vi[t1]:=n; fi;

%p else

%p    v0:=[op(v0),t2]; s1:={op(s1),t2}; v1:=[op(v1),t2]; v2:=[op(v2),r1]; r1:=1; t1:=t2;

%p    if t1<=M2 then vi[t1]:=n; fi;

%p fi;

%p od:

%p # A114183:

%p [seq(v0[i],i=1..nops(v0))];

%p # A221715:

%p [seq(v1[i],i=1..nops(v1))];

%p # A221716:

%p [seq(v2[i],i=1..nops(v2))];

%p # A189419:

%p [seq(vi[i],i=1..M2)];

%Y Cf. A114183, A189419, A221716, A213220.

%K nonn

%O 1,2

%A _N. J. A. Sloane_, Jan 27 2013