proposed
approved
proposed
approved
editing
proposed
(Python)
from sympy import divisors
def A000120(n): return bin(n).count('1')
def aupto(limit):
alst = []
for m in range(1, limit+1):
if A000120(m) == sum(A000120(d) for d in divisors(m)[:-1]): alst += [m]
return alst
print(aupto(1563)) # Michael S. Branicky, Feb 25 2021
approved
editing
Amiram Eldar, <a href="/A175522/b175522_1.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Reinhard Zumkeller)
reviewed
approved
proposed
reviewed
editing
proposed
binw[n_] := DigitCount[n, 2, 1]; Select[Range[1500], binw[#] == DivisorSum[#, binw[#1] &]/2 &] (* Amiram Eldar, Dec 14 2020 *)
Amiram Eldar, <a href="/A175522/b175522_1.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Reinhard Zumkeller)
Reinhard Zumkeller, Amiram Eldar, <a href="/A175522/b175522_1.txt">Table of n, a(n) for n = 1..100010000</a>
approved
editing