A318734 a(n) = Sum_{k=1..n} (-1)^(k + 1) * d(2*k - 1), where d(k) is the number of divisors of k (A000005).
1, -1, 1, -1, 2, 0, 2, -2, 0, -2, 2, 0, 3, -1, 1, -1, 3, -1, 1, -3, -1, -3, 3, 1, 4, 0, 2, -2, 2, 0, 2, -4, 0, -2, 2, 0, 2, -4, 0, -2, 3, 1, 5, 1, 3, -1, 3, -1, 1, -5, -3, -5, 3, 1, 3, -1, 1, -3, 3, -1, 2, -2, 2, 0, 4, 2, 6, -2, 0, -2, 2, -2
Offset: 1
Links
- Hugo Pfoertner, Table of n, a(n) for n = 1..10000
- Hugo Pfoertner, Records of A318734, showing negative and positive records.
Crossrefs
Programs
-
Mathematica
a[n_] := Sum[(-1)^(k + 1) DivisorSigma[0, 2 k - 1], {k, 1, n}]; Array[a, 100] (* Jean-François Alcover, Sep 17 2018 *) -
PARI
s=0;j=-1;forstep(k=1,141,2,j=-j;s=s+j*numdiv(k);print1(s,", "))
-
PARI
a(n) = sum(k=1, n, (-1)^(k+1)*numdiv(2*k-1)); \\ Michel Marcus, Sep 08 2018
Comments