OFFSET
1,4
LINKS
MAPLE
for n from 1 to 1000 do
L:= convert(n, base, 2);
P:= add(L[i]*X^(i-1), i=1..nops(L));
R:= Factors(P) mod 2;
a[n]:= add(r[2], r=R[2]);
od:
seq(a[n], n=1..1000); # Robert Israel, Jun 07 2015
PROG
(PARI) a(n)=my(fm=factor(Pol(binary(n))*Mod(1, 2))); sum(k=1, #fm~, fm[k, 2]) \\ Franklin T. Adams-Watters, Jun 07 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 03 2004
STATUS
approved