login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A309668
a(n) is the least positive number of the form floor(2^k/n) for some k >= 0 not yet in the sequence.
2
1, 2, 5, 4, 3, 10, 9, 8, 7, 6, 11, 21, 19, 18, 17, 16, 15, 14, 13, 12, 24, 23, 22, 42, 20, 39, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 52, 25, 49, 48, 47, 46, 45, 44, 43, 85, 41, 40, 80, 78, 38, 75, 74, 73, 71, 70, 69, 68, 67, 66, 65, 64, 63, 62, 61
OFFSET
1,2
COMMENTS
The sequence is well defined as for any n > 0, there are infinitely many positive numbers of the form floor(2^k/n) with k >= 0.
The sequence is a permutation of the natural numbers, with inverse A309734:
- for any m > 0, floor(2^k/A300475(m)) = m for some k,
- also, for any u > 0, floor(2^(k-u)/(A300475(m)*2^u)) = m,
- so the set S_m = { v such that floor(2^k/v) = m for some k >= 0 } is infinite
- and eventually a(n) = m for some n in S_m, QED.
EXAMPLE
The first terms, alongside the binary representations of a(n) and of 1/n (with that of a(n) in parentheses), are:
-- ---- --------- ---------------------
1 1 1 (1).00000000000000...
2 2 10 0.(10)000000000000...
3 5 101 0.0(101)0101010101...
4 4 100 0.0(100)0000000000...
5 3 11 0.00(11)0011001100...
6 10 1010 0.00(1010)10101010...
7 9 1001 0.00(1001)00100100...
8 8 1000 0.00(1000)00000000...
9 7 111 0.000(111)00011100...
10 6 110 0.000(110)01100110...
11 11 1011 0.000(1011)1010001...
12 21 10101 0.000(10101)010101...
13 19 10011 0.000(10011)101100...
14 18 10010 0.000(10010)010010...
15 17 10001 0.000(10001)000100...
PROG
(PARI) s=1; for (n=1, 67, q=1/n; while (bittest(s, f=floor(q)), q*=2); print1 (f ", "); s+=2^f)
CROSSREFS
Cf. A070939, A300475, A309734 (inverse).
Sequence in context: A132664 A072029 A309734 * A238758 A354182 A065652
KEYWORD
nonn,look,base
AUTHOR
Rémy Sigrist, Aug 11 2019
STATUS
approved