|
|
A297841
|
|
a(n) = Sum_{d|n} max(d, n/d)^2.
|
|
4
|
|
|
1, 8, 18, 36, 50, 90, 98, 160, 171, 250, 242, 392, 338, 490, 500, 656, 578, 882, 722, 1050, 980, 1210, 1058, 1640, 1275, 1690, 1620, 2058, 1682, 2522, 1922, 2688, 2420, 2890, 2548, 3726, 2738, 3610, 3380, 4328, 3362, 4900, 3698, 5082, 4662, 5290, 4418, 6688
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
|
|
MAPLE
|
f:= n -> add(max(d, n/d)^2, d = numtheory:-divisors(n)):
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|