|
%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,changed
%O 0,6
%A _Omar E. Pol_, Sep 14 2019
|
