OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
6 is in this sequence because 6 divides sigma(6^2) - 1 = 90.
MAPLE
select(n -> type((numtheory:-sigma(n^2)-1)/n, integer), [$1..1000]); # Robert Israel, Jun 20 2018
MATHEMATICA
Select[Range@ 270, Divisible[DivisorSigma[1, #^2] - 1, #] &] (* Michael De Vlieger, Mar 15 2017 *)
PROG
(Magma) [n: n in [1..200] | Denominator(((SumOfDivisors(n^2))-1)/n) eq 1];
(PARI) for(n=1, 269, if((sigma(n^2) - 1)%n==0, print1(n, ", "))) \\ Indranil Ghosh, Mar 14 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Mar 14 2017
STATUS
approved