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”).
%I #15 Feb 01 2019 05:04:06
%S 1,2,3,6,4,5,20,8,9,72,16,7,14,21,78,32,11,352,64,10,40,15,24,12,30,
%T 35,390,48,96,51,102,60,13,832,117,144,18,168,42,28,39,180,56,84,63,
%U 70,780,120,26,128,19,504,36,288,126,45,112,151,896,156,720,224
%N Lexicographically earliest sequence of distinct positive terms such that for any n > 0 and any k > 0, floor((2^k) / a(n)) AND floor((2^k) / a(n+1)) = 0 (where AND denotes the bitwise AND operator).
%C In other words, for any n > 0, the binary expansions of 1/a(n) and of 1/a(n+1) have no common one bit; in this sense, this sequence is similar to A109812.
%C This sequence is a permutation of the natural numbers, with inverse A306233 (we can first prove that all the powers of 2 appear in the sequence and then that every natural number appear in the sequence).
%H Rémy Sigrist, <a href="/A306231/b306231.txt">Table of n, a(n) for n = 1..2000</a>
%H Rémy Sigrist, <a href="/A306231/a306231.gp.txt">PARI program for A306231</a>
%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
%F For any n > 0, if A000120(a(n)) <> 1 and A000120(a(n+1)) <> 1, then gcd(A007733(a(n)), A007733(a(n+1))) > 1.
%e The first terms, alongside A007733(a(n)) and the binary representation of 1/a(n) with periodic part in parentheses, are:
%e n a(n) period bin(1/a(n))
%e -- ---- ------ -------------------
%e 1 1 1 1.(0)
%e 2 2 1 0.1(0)
%e 3 3 2 0.(01)
%e 4 6 2 0.0(01)
%e 5 4 1 0.01(0)
%e 6 5 4 0.(0011)
%e 7 20 4 0.00(0011)
%e 8 8 1 0.001(0)
%e 9 9 6 0.(000111)
%e 10 72 6 0.000(000111)
%e 11 16 1 0.0001(0)
%e 12 7 3 0.(001)
%e 13 14 3 0.0(001)
%e 14 21 6 0.(000011)
%e 15 78 12 0.0(000001101001)
%e 16 32 1 0.00001(0)
%e 17 11 10 0.(0001011101)
%e 18 352 10 0.00000(0001011101)
%e 19 64 1 0.000001(0)
%e 20 10 4 0.0(0011)
%o (PARI) See Links section.
%Y Cf. A000120, A007733, A109812, A306233 (inverse).
%K nonn,base
%O 1,2
%A _Rémy Sigrist_, Jan 30 2019