OFFSET
1,1
COMMENTS
EXAMPLE
16560 is a term because sigma(16560) = 58032 and psi(16560) = 41472; sigma(16560) - psi(16560) = 16560.
PROG
(PARI) a001615(n) = n*sumdivmult(n, d, issquarefree(d)/d);
isok(n) = sigma(n)==a001615(n)+n; \\ after Charles R Greathouse IV at A001615
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Aug 21 2017
EXTENSIONS
a(9)-a(11) from Giovanni Resta, Aug 21 2017
STATUS
approved