A221715 Start of n-th doubling run in A114183.

1, 5, 3, 9, 33, 11, 13, 7, 21, 25, 20, 17, 23, 19, 49, 39, 35, 47, 27, 29, 15, 43, 37, 97, 55, 41, 51, 57, 85, 73, 193, 111, 59, 61, 31, 89, 53, 329, 145, 192, 221, 237, 87, 105, 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

Offset: 1

Author

N. J. A. Sloane, Jan 27 2013

Keywords

• nonn

Links

• N. J. A. Sloane, Table of n, a(n) for n = 1..7328

Crossrefs

Cf. A114183, A189419, A221716, A213220.

Programs

Maple

# A114181
M:=10000; M2:=1000;
s1:={1}; v0:=[1]; v1:=[1]; v2:=[]; vi:=Array(1..M2);
t1:=1; r1:=1; vi[1]:=1;
for n from 2 to M do
t2:=floor(sqrt(t1));
if t2 in s1 then
   v0:=[op(v0),2*t1]; s1:={op(s1),2*t1}; r1:=r1+1; t1:=2*t1;
   if t1<=M2 then vi[t1]:=n; fi;
else
   v0:=[op(v0),t2]; s1:={op(s1),t2}; v1:=[op(v1),t2]; v2:=[op(v2),r1]; r1:=1; t1:=t2;
   if t1<=M2 then vi[t1]:=n; fi;
fi;
od:
# A114183:
[seq(v0[i],i=1..nops(v0))];
# A221715:
[seq(v1[i],i=1..nops(v1))];
# A221716:
[seq(v2[i],i=1..nops(v2))];
# A189419:
[seq(vi[i],i=1..M2)];