OFFSET
1,2
COMMENTS
The prime terms are given by A143853.
PROG
(PARI) s=0; pp=0; for(n=2, 10^8, p=factor(n)[, 1]; s += (n-2) - pp - sigma(n) + sum(i=1, #p, p[i]) + if(!ispseudoprime(n), n, pp++; 0) + 1; if(s%n==0, print1(n, ", ")) )
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Max Alekseyev, Dec 07 2013
STATUS
approved