OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
MAPLE
f:= n -> add(max(d, n/d)^2, d = numtheory:-divisors(n)):
map(f, [$1..100]); # Robert Israel, Jan 10 2018
MATHEMATICA
f[n_] := Block[{d = Divisors@ n}, Plus @@ (Max[#, n/#]^2 & /@ d)]; Array[f, 50] (* Robert G. Wilson v, Jan 07 2018 *)
PROG
(PARI) {a(n) = sumdiv(n, d, max(d, n/d)^2)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 07 2018
STATUS
approved