OFFSET
1,2
COMMENTS
If p is prime, a(p) = Sum_{d|p} d * sopf(d) = p * p = p^2.
EXAMPLE
a(10) = Sum_{d|10} d * sopf(d) = 1*0 + 2*2 + 5*5 + 10*7 = 99.
MATHEMATICA
Table[Sum[i*Sum[k (PrimePi[k] - PrimePi[k - 1]) (1 - Ceiling[i/k] + Floor[i/k]), {k, i}] (1 - Ceiling[n/i] + Floor[n/i]), {i, n}], {n, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, May 23 2021
STATUS
approved