OFFSET
0,3
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 0..1000
Project Euler, Problem 759, sequence f(n).
FORMULA
a(2*n) = 2*a(n).
a(2*n+1) = 2*n + 1 + (2+1/n)*a(n). - Robert Israel, Aug 01 2014
G.f.: x * (d/dx) (1/(1 - x))*Sum_{k>=0} x^(2^k)/(1 + x^(2^k)). - Ilya Gutkovskiy, Mar 27 2018
EXAMPLE
G.f. = x + 2*x^2 + 6*x^3 + 4*x^4 + 10*x^5 + 12*x6 + 21*x^7 + 8*x^8 + 18*x^9 + ...
MAPLE
a:= n -> n * convert(convert(n, base, 2), `+`):
seq(a(n), n=0..100); # Robert Israel, Aug 01 2014
MATHEMATICA
Table[n*DigitCount[n, 2, 1], {n, 0, 100}] (* Harvey P. Dale, Dec 16 2014 *)
PROG
(PARI) sumbit(n) = my(r); while(n>0, r+=n%2; n\=2); r
a(n) = n*sumbit(n)
(PARI) {a(n) = if( n<0, 0, n * sumdigits(n, 2))}; /* Michael Somos, Aug 05 2014 */ /* since version 2.6.0 */
(Python) [n*bin(n)[2:].count('1') for n in range(1000)] # Chai Wah Wu, Aug 03 2014
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Franklin T. Adams-Watters, Aug 01 2014
STATUS
approved