OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Donovan Johnson)
EXAMPLE
54 is not a square, sigma(54)=120, 6 is the squarefree part of 54 and 6 divides 120. Hence a(5)=54 is in the sequence.
MATHEMATICA
f[p_, e_] := p^Mod[e, 2]; sqfp[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[2000], (s = sqfp[#]) > 1 && Divisible[DivisorSigma[1, #], s] &] (* Amiram Eldar, May 10 2022 *)
PROG
(PARI) c=0; for(n=2, 244783, if(issquare(n)==0, if(sigma(n)%core(n)==0, c++; write("b069235.txt", c " " n)))) /* Donovan Johnson, Aug 18 2013 */
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 13 2002
STATUS
approved