OFFSET
1,1
COMMENTS
Odd terms in A094694.
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000
MAPLE
select(t -> convert(convert(t^2, base, 2), `+`) < convert(convert(t, base, 2), `+`), [seq(i, i=1..10^4, 2)]); # Robert Israel, Mar 13 2022
MATHEMATICA
Select[Range[1, 6000, 2], Greater @@ DigitCount[{#, #^2}, 2, 1] &] (* Amiram Eldar, Mar 11 2022 *)
PROG
(PARI) forstep(n=1, 10^4, 2, if(hammingweight(n^2)<hammingweight(n), print1(n, ", ")));
(Python)
def ok(n): return n%2 == 1 and bin(n).count('1') > bin(n**2).count('1')
print([k for k in range(6000) if ok(k)]) # Michael S. Branicky, Mar 11 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Joerg Arndt, Mar 11 2022
STATUS
approved