OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
core[n_] := Times @@ (First[#]^Mod[Last[#], 2]& /@ FactorInteger[n]); b[n_] := Times @@ DigitCount[n, 2, {0, 1}]; Select[Range[10^5], core[#] == b[#] &] (* Amiram Eldar, Sep 03 2020 *)
PROG
(PARI) for(s=1, 500000, b=binary(s); l=length(b); if(sum(i=1, l, if(component(b, i)-1, 0, 1))*sum(i=1, l, if(component(b, i), 0, 1))==core(s), print1(s, ", ")))
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Benoit Cloitre, Jun 22 2002
EXTENSIONS
More terms from Amiram Eldar, Sep 03 2020
STATUS
approved