A063880 Numbers k such that sigma(k) = 2*usigma(k).

108, 540, 756, 1188, 1404, 1836, 2052, 2484, 3132, 3348, 3780, 3996, 4428, 4644, 5076, 5724, 5940, 6372, 6588, 7020, 7236, 7668, 7884, 8316, 8532, 8964, 9180, 9612, 9828, 10260, 10476, 10908, 11124, 11556, 11772, 12204, 12420, 12852, 13716, 14148

Offset

1,1

Comments

Numbers so far are all == 108 (mod 216). - Ralf Stephan, Jul 07 2003 [Confirmed up to 10^7 by Robert G. Wilson v.]

Also numbers whose unitary and nonunitary divisors have equal sum. - Amiram Eldar, Sep 30 2019

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)

Mathematica

usigma[n_] := Block[{d = Divisors[n]}, Plus @@ Select[d, GCD[ #, n/# ] == 1 &]]; Select[ Range[14363], DivisorSigma[1, # ] == 2 usigma[ # ] &] (* Robert G. Wilson v, Aug 28 2004 *)

Prog

(PARI) u(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d));
j=[];  for(n=1, 30000, if(sigma(n) == 2*u(n), j=concat(j, n))); j

Crossrefs

Cf. A000203, A048146, A034448, A097702.

Sequence in context: A224526 A251270 A255083 * A059436 A233685 A160919

Adjacent sequences: A063877 A063878 A063879 * A063881 A063882 A063883

Keyword

easy,nonn,changed

Author

Jason Earls, Aug 27 2001

Status

approved