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,2

Links

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

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)];

Crossrefs

Cf. A114183, A189419, A221716, A213220.

Sequence in context: A082454 A108245 A228789 * A248660 A141620 A195140

Adjacent sequences: A221712 A221713 A221714 * A221716 A221717 A221718

Keyword

nonn

Author

N. J. A. Sloane, Jan 27 2013

Status

approved