OFFSET
1,1
COMMENTS
MATHEMATICA
q[k_] := Module[{h = DigitCount[k, 2, 1]}, DivisorSum[k, # &, DigitCount[#, 2, 1] == h &] == 4*k]; Select[Range[3*10^8], q]
PROG
(PARI) isok(k) = {my(h = hammingweight(k)); sumdiv(k, d, d*(hammingweight(d) == h)) == 4*k; }
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amiram Eldar, Feb 05 2025
EXTENSIONS
a(19)-a(22) from Jinyuan Wang, Feb 12 2025
STATUS
approved