OFFSET
1,2
COMMENTS
LINKS
Michel Marcus, Table of n, a(n) for n = 1..5000
EXAMPLE
a(15) = 24 - 6 = 18 because the sum of the divisors of 15 is 1 + 3 + 5 + 15 = 24 and there are 6 primes not exceeding 15 (2, 3, 5, 7, 11 and 13). - Emeric Deutsch, Dec 29 2008
MAPLE
with(numtheory): seq(sigma(n)-pi(n), n = 1 .. 80); # Emeric Deutsch, Dec 29 2008
MATHEMATICA
Table[DivisorSigma[1, n]-PrimePi[n], {n, 80}] (* Harvey P. Dale, Oct 20 2021 *)
PROG
(PARI) a(n) = sigma(n) - primepi(n); \\ Michel Marcus, Jun 18 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Omar E. Pol, Dec 19 2008
EXTENSIONS
Corrected and extended by Emeric Deutsch, Dec 29 2008
STATUS
approved