OFFSET
1,1
FORMULA
PROG
(PARI) for(n=1, 1000000, if(sumdiv(n, d, d)-n+1-sumdiv(n, d, isprime(d)*d^2)==0, print1(n, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, May 07 2002
STATUS
approved