|
%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
|
