OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
FORMULA
From Michel Marcus, Jan 09 2018: (Start)
a(p) = 2, for prime p.
a(2p) = 10, for odd prime p. (End)
a(r*p) = 2 + 2*r^2, where r, p are primes and p > r. - Bruno Berselli, Jan 09 2018
MATHEMATICA
a[n_] := Block[{d = Divisors@n}, Plus @@ (Min[#, n/#]^2 & /@ d)]; Array[a, 70] (* Robert G. Wilson v, Jan 06 2018 *)
PROG
(PARI) {a(n) = sumdiv(n, d, min(d, n/d)^2)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 06 2018
STATUS
approved