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
Keywords
Links
- N. J. A. Sloane, Table of n, a(n) for n = 1..7328
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)];
Comments