A063749 a(n) = floor((A000005(n)*(n+1)/2) - A000203(n)).

0, 0, 0, 0, 0, 2, 0, 3, 2, 4, 0, 11, 0, 6, 8, 11, 0, 18, 0, 21, 12, 10, 0, 40, 8, 12, 16, 31, 0, 52, 0, 36, 20, 16, 24, 75, 0, 18, 24, 74, 0, 76, 0, 51, 60, 22, 0, 121, 18, 60, 32, 61, 0, 100, 40, 108, 36, 28, 0, 198, 0, 30, 88, 100, 48, 124, 0, 81, 44, 140, 0, 243, 0, 36, 104

Offset

1,6

Comments

a(n) = 0 if n is prime, otherwise >= 0; first part of formula provides estimation of sigma(n).

Links

Harry J. Smith, Table of n, a(n) for n = 1..1000

Example

n = 4: tau(4) = 3, (3*(4+1)/2) - 7 = 0.5 so a(4) = 0;
n = 12: tau(12) = 6, (6*(12+1)/2) - 28 = 39 - 28 = 11 = a(12).

Prog

(PARI) { for (n=1, 1000, a = numdiv(n)*(n + 1)\2 - sigma(n); write("b063749.txt", n, " ", a) ) } \\ Harry J. Smith, Aug 29 2009

Crossrefs

Cf. A000005, A000203.

Sequence in context: A140944 A057860 A092915 * A231333 A357866 A331622

Adjacent sequences: A063746 A063747 A063748 * A063750 A063751 A063752

Keyword

nonn

Author

Labos Elemer, Aug 14 2001

Status

approved