login
a(n) = A006218(n) - A005187(n).
1

%I #13 Apr 18 2024 04:22:18

%S 0,0,0,1,1,2,4,5,5,7,9,10,13,14,16,19,19,20,24,25,28,31,33,34,38,40,

%T 42,45,48,49,55,56,56,59,61,64,70,71,73,76,80,81,87,88,91,96,98,99,

%U 104,106,110,113,116,117,123,126,130,133,135,136,145,146,148,153,153,156,162,163,166,169,175,176,184

%N a(n) = A006218(n) - A005187(n).

%C Also, 0 together with the partial sums of A326987.

%F a(n) ~ n * (log(n) + 2*gamma - 3), where gamma is Euler's constant (A001620). - _Amiram Eldar_, Apr 18 2024

%t With[{nmax = 100}, Join[{0}, Accumulate[Table[DivisorSigma[0, n], {n, 1, nmax}]]] - Table[2*n - DigitCount[n, 2, 1], {n, 0, nmax}]] (* _Amiram Eldar_, Apr 18 2024 *)

%Y Cf. A000005, A001620, A005187, A006218, A326987, A327327.

%K nonn

%O 0,6

%A _Omar E. Pol_, Sep 14 2019