OFFSET
1,2
COMMENTS
All elements except the first, a(1)=1, are of the form 2*p, where p is a prime and p == 1 (mod 12). Also, sigma(2*p) = (1+2)*(1+p) = 6m where m = (1+p)/2 and m == 1 (mod 6). A squarefree composite number not of the form 2*p cannot be in the sequence since sigma is multiplicative. For example, sigma(p*q) = (1+p)*(1+q) is divisible by 4 for p,q > 2. - Walter Kehowski, Mar 21 2007
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
PROG
(PARI) { n=0; for (m = 1, 10^9, if (!isprime(m) && moebius(m) && moebius(sigma(m)), write("b065302.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, Oct 16 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Oct 29 2001
STATUS
approved