OFFSET
0,3
MAPLE
a := n -> n * binomial(ilog2(n) + 1, 2):
seq(a(n), n = 0..56); # Peter Luschny, Jul 19 2024
MATHEMATICA
Table[n * Binomial[BitLength[n], 2], {n, 0, 60}] (* Amiram Eldar, Jul 23 2024 *)
PROG
(Python)
def a(n):
q = n.bit_length() - 1
return n * q * (q + 1) >> 1
print([a(n) for n in range(61)])
(PARI) a(n) = n * binomial(logint(n, 2) + 1, 2); \\ Michel Marcus, Jul 23 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Darío Clavijo, Jul 17 2024
STATUS
approved