OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harvey P. Dale)
EXAMPLE
1369 has 3 divisors which sum to 1407; 1407 is divisible by 3; the sum of the squares of the divisors of 1369 is 1875531 which is divisible by 1407; so 1369 is a term of the sequence.
MATHEMATICA
Select[Range[50000], Divisible[DivisorSigma[1, #], DivisorSigma[0, #]] && Divisible[ DivisorSigma[2, #], DivisorSigma[1, #]]&]
PROG
(PARI) isok(k) = {my(f = factor(k), d = numdiv(f), s = sigma(f), s2 = sigma(f, 2)); !(s % d) && !(s2 % s); } \\ Amiram Eldar, Jan 25 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Harvey P. Dale, Oct 19 2016
STATUS
approved
