OFFSET
1,2
COMMENTS
This sequence is a permutation of the natural numbers, with inverse A297499; as a(1) = 1, for any k > 1, a(2^k) has only to be a multiple of 1, and so a(2^k) will be the least unused value, and eventually any number will appear in the sequence.
Prime numbers can only appear at positions that are powers of 2.
For any n > 1, a(n) is a multiple of m(n) = lcm(a(e_1), ..., a(e_h)) where the list (e_1, ..., e_h) corresponds to the ones in the binary expansion of n (in particular, e_1 = 1 and h = A000120(n)); the lines and dashed lines visible in the logarithmic scatterplot of the first terms correspond to sets of terms a(n) where m(n) has the same value (see Links section).
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..50000
Rémy Sigrist, Colored logarithmic scatterplot of the first 150000 terms (where the color is function of m(n))
Rémy Sigrist, PARI program for A295029
EXAMPLE
The first terms, alongside the binary expansion of n and m(n), are:
n a(n) bin(n) m(n)
-- ---- ------ -------------
1 1 1 1 = lcm(a(1))
2 2 10 1 = lcm(a(1))
3 4 11 2 = lcm(a(1), a(2))
4 3 100 1 = lcm(a(1))
5 8 101 4 = lcm(a(1), a(3))
6 6 110 2 = lcm(a(1), a(2))
7 12 111 4 = lcm(a(1), a(2), a(3))
8 5 1000 1 = lcm(a(1))
9 9 1001 3 = lcm(a(1), a(4))
10 16 1010 4 = lcm(a(1), a(3))
11 24 1011 12 = lcm(a(1), a(3), a(4))
12 10 1100 2 = lcm(a(1), a(2))
13 18 1101 6 = lcm(a(1), a(2), a(4))
14 20 1110 4 = lcm(a(1), a(2), a(3))
15 36 1111 12 = lcm(a(1), a(2), a(3), a(4))
16 7 10000 1 = lcm(a(1))
17 32 10001 8 = lcm(a(1), a(5))
18 15 10010 3 = lcm(a(1), a(4))
19 48 10011 24 = lcm(a(1), a(4), a(5))
20 28 10100 4 = lcm(a(1), a(3))
PROG
(PARI) See Links section.
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Rémy Sigrist, Dec 30 2017
STATUS
approved