|
|
|
|
0, 0, 1, 2, 5, 8, 12, 16, 23, 30, 38, 46, 56, 66, 77, 88, 103, 118, 134, 150, 168, 186, 205, 224, 246, 268, 291, 314, 339, 364, 390, 416, 447, 478, 510, 542, 576, 610, 645, 680, 718, 756, 795, 834, 875, 916, 958, 1000, 1046, 1092, 1139, 1186, 1235, 1284, 1334
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
a(n) ~ (n^2)/2 + O(n*log_2(n)). [Lagarias and Mehta, theorem 4.2 with p=2]
a(n) = ( (n+1)^2 - Sum_{i=1..k} (e[i]+2*i-1) * 2^e[i] )/2, where binary expansion n+1 = 2^e[1] + ... + 2^e[k] with descending exponents e[1] > e[2] > ... > e[k] (A272011).
(End)
|
|
MAPLE
|
a:= proc(n) option remember; `if`(n<1, 0,
a(n-1)+n-add(i, i=Bits[Split](n)))
end:
|
|
MATHEMATICA
|
Accumulate[Table[n-DigitCount[n, 2, 1], {n, 0, 130}]] (* Harvey P. Dale, Feb 26 2015 *)
|
|
PROG
|
(PARI) a(n) = n++; my(v=binary(n), t=#v-1); for(i=1, #v, if(v[i], v[i]=t++, t--)); (n^2 - fromdigits(v, 2))>>1; \\ Kevin Ryde, Oct 29 2021
|
|
CROSSREFS
|
Cf. A249152 (hyperfactorial valuation), A187059 (binomial valuation), A173345 (superfactorial 10-valuation).
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|