OFFSET
1,1
EXAMPLE
For k=10, sigma(10)^2 / (10-1) = 18^2 / 9 = 324 / 9 = 36.
MATHEMATICA
Select[Range[2, 10^5], Divisible[DivisorSigma[1, #]^2, # - 1] &] (* Amiram Eldar, May 15 2021 *)
PROG
(PARI)
list(nn) = for(n=2, nn, if (sigma(n)^2 % (n-1) == 0, print1(n, ", ")))
list(100000)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zdenek Cervenka, May 15 2021
STATUS
approved