OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
7 is in the sequence because 7 has 2 divisors {1,7} therefore 5 non-divisors {2,3,4,5,6}, 2^2 + 3^2 + 4^2 + 5^2 + 6^2 = 90 and 5 divides 90.
MATHEMATICA
Select[Range[800], Mod[#1 (#1 + 1) ((2 #1 + 1)/6) - DivisorSigma[2, #1], #1 - DivisorSigma[0, #1]] == 0 & ]
PROG
(PARI) is(n)=my(f=factor(n)); n>2 && ((2*n^3+3*n^2+n)/6-sigma(f, 2))%(n-numdiv(f))==0 \\ Charles R Greathouse IV, Dec 19 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ilya Gutkovskiy, Dec 19 2016
STATUS
approved