OFFSET
0,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 0..99
EXAMPLE
MATHEMATICA
q[k_] := DivisorSum[k, 1 &, DigitCount[#, 2, 1] == DigitCount[k, 2, 1] &] == 1; seq[len_] := Module[{s = Table[0, {len}], c = 0, k = 1, i}, While[c < len, i = DivisorSigma[0, k]; If[i <= len && s[[i]] == 0 && q[k], c++; s[[i]] = k]; k++]; s]; seq[16]
PROG
(PARI) is1(k) = {my(h = hammingweight(k)); sumdiv(k, d, hammingweight(d) == h) == 1};
list(len) = {my(s = vector(len), c = 0, k = 1, i); while(c < len, i = numdiv(k); if(i <= len && s[i] == 0 && is1(k), c++; s[i] = k); k++); s; }
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amiram Eldar, Apr 24 2025
STATUS
approved
