OFFSET
1,2
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..71 (terms < 2*10^9)
EXAMPLE
Divisors of 3145 are 1, 5, 17, 37, 85, 185, 629, 3145.
sigma(1) = 1, sigma(5) = 6, sigma(17) = 18, sigma(37) = 38, sigma(85) = 108, sigma(185) = 228, sigma(629) = 684, sigma(3145) = 4104.
(1^2 + 6^2 + 18^2 + 38^2 + 108^2 + 228^2 + 684^2 + 4104^2) / 3145 = 5525.
MAPLE
with(numtheory); ListA226563:=proc(q) local a, b, k, n;
for n from 1 to q do a:=[op(divisors(n))]; b:=add(sigma(a[k])^2, k=1..nops(a));
if type(b/n, integer) then print(n); fi; od; end: ListA226563(10^6);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Jun 11 2013
EXTENSIONS
a(14)-a(33) from Giovanni Resta, Jun 11 2013
STATUS
approved