1,1
(PARI) u(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d)); z(n) = sigma(n)-u(n); for(n=1, 10^7, if(z(sigma(n))==2*n, print(n)))
Cf. A000203, A034448, A048146, A019284.
Sequence in context: A099060 A035174 A203973 * A019284 A112389 A187634
Adjacent sequences: A063882 A063883 A063884 * A063886 A063887 A063888
more,nonn
Jason Earls (zevi_35711(AT)yahoo.com), Aug 28 2001