OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..300
PROG
(PARI) {for(n=1, 22000000, d=numdiv(n); if(d==sigma(d)-d, s=0; F=factor(n); for(i=1, matsize(F)[1], s=s+F[i, 1]*F[i, 2]); if(s==sigma(s)-s, print1(n, ", "))))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, May 10 2002
EXTENSIONS
More terms from Klaus Brockhaus, May 15 2002
STATUS
approved