OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..250 (terms 1..100 from Harvey P. Dale)
MAPLE
q:= n-> is(irem(numtheory[sigma](n^2+1), n)=0):
select(q, [$1..100000])[]; # Alois P. Heinz, Jan 26 2023
MATHEMATICA
Select[Range[50000], Divisible[DivisorSigma[1, #^2+1], #]&] (* Harvey P. Dale, Nov 04 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 05 2002
STATUS
approved