|
|
A306645
|
|
a(n) is the least positive multiple of n belonging to A306263 if any, or a(n) = -1 otherwise.
|
|
1
|
|
|
1, 2, 6, 4, 10, 6, 42, 8, 18, 10, 66, 12, 156, 42, 60, 16, 34, 18, 228, 20, 42, 66, 92, 24, 300, 156, 108, 84, 116, 60, 310, 32, 66, 34, 420, 36, 222, 228, 156, 40, 246, 42, 172, 132, 180, 92, 2820, 48, 588, 300, 204, 156, 212, 108, 660, 168, 228, 116, 590, 60
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Is a(n) > 0 for every n > 0?
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
For n = 7:
- the divisors of 7 are: 1, 7,
- the corresponding Hamming weights are: 1, 3,
- 3 does not divide 7,
- the divisors of 3*7 are: 1, 3, 7, 21,
- the corresponding Hamming weights are: 1, 2, 3, 3,
- 2 does not divide 3*7,
- the divisors of 2*3*7 are: 1, 2, 3, 6, 7, 14, 21, 42,
- the corresponding Hamming weights are: 1, 1, 2, 2, 3, 3, 3, 3,
- they all divide 2*3*7,
- hence a(7) = 2*3*7 = 42.
|
|
MATHEMATICA
|
With[{s = Select[Range[3000], With[{k = #}, AllTrue[Divisors@ k, Mod[k, DigitCount[#, 2, 1]] == 0 &]] &]}, Table[SelectFirst[s, Mod[#, n] == 0 &] /. k_ /; MissingQ@ k -> -1, {n, 60}]] (* Michael De Vlieger, Mar 05 2019 *)
|
|
PROG
|
(PARI) a(n) = while (1, my (m=n); fordiv (m, d, m = lcm(m, hammingweight(d)); ); if (n==m, return (n), n = m))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|