OFFSET
1,1
COMMENTS
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
8 is in the sequence because 8 has 3 proper divisors {1,2,4}, 1^2 + 2^2 + 4^2 = 21 and 3 divides 21.
MATHEMATICA
Select[Range[150], Mod[DivisorSigma[2, #1] - #1^2, DivisorSigma[0, #1] - 1] == 0 &]
Select[Range[200], IntegerQ[Mean[Most[Divisors[#]]^2]]&] (* Harvey P. Dale, Jul 26 2019 *)
PROG
(PARI) is(n)=my(d=divisors(n)); d=apply(k->k^2, d[1..#d-1]); n>1 && sum(i=1, #d, d[i])%#d==0 \\ Charles R Greathouse IV, Dec 19 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ilya Gutkovskiy, Dec 19 2016
STATUS
approved