OFFSET
1,2
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
EXAMPLE
MATHEMATICA
binPalQ[n_] := PalindromeQ[IntegerDigits[n, 2]]; Select[Range[10^6], AllTrue[{#, Plus @@ IntegerDigits[#], DivisorSigma[1, #] - #}, binPalQ] &] (* Amiram Eldar, Jul 28 2023 *)
PROG
(Python)
from sympy import divisor_sigma as s1
def bp(n): return (b:=bin(n)[2:])==b[::-1]
def ok(n): return n and bp(n) and bp(sum(map(int, str(n)))) and bp(s1(n)-n)
print([k for k in range(10**6) if ok(k)]) # Michael S. Branicky, Jul 24 2023
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Aakash Gautam, Jul 12 2023
EXTENSIONS
a(6) and beyond from Michael S. Branicky, Jul 24 2023
STATUS
approved