login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A019278 Numbers n such that sigma(sigma(n)) = k*n for some k. 25
1, 2, 4, 8, 15, 16, 21, 24, 42, 60, 64, 84, 160, 168, 240, 336, 480, 504, 512, 960, 1023, 1344, 1536, 4092, 4096, 10752, 13824, 16368, 29127, 32256, 32736, 47360, 57120, 58254, 61440, 65472, 65536, 86016, 116508, 217728, 262144, 331520, 343976, 466032, 550095 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Let sigma_m (n) be result of applying sum-of-divisors function (A000203) m times to n; call n (m,k)-perfect if sigma_m (n) = k*n; then this is the sequence of (2,k)-perfect numbers.

From Michel Marcus, May 14 2016: (Start)

For these numbers, the quotient q = sigma(sigma(n))/n is an integer (see A098223). Then also q = (sigma(s)/s)*(sigma(n)/n) with s = sigma(n). That is, q = abundancy(s)*abundancy(n).

So looking at the abundancy of these terms may be interesting. Indeed we see that 459818240 and 51001180160 are actually 3-perfect numbers (A005820), and the reason they are here is that they are coprime to 3. So their sums of divisors are 4-perfect numbers (A027687), yielding q=12.

In a similar way, we can see that the 5-perfect numbers (A046060) that are coprime to 5 will be terms of this sequence with q=30. There are 20 such numbers, the smallest being 13188979363639752997731839211623940096. (End)

From Michel Marcus, May 15 2016: (Start)

It is also interesting to note that for a(2)=8, s=sigma(8)=15 is also a term. This happens to be the case for chains of several terms in a row:

8, 15, 24, 60, 168, 480 with q = 3,4,7,8,9,10;

512, 1023, 1536, 4092, 10752, 32736 with q = 3,4,7,8,9,10;

29127, 47360, 116508, 331520, 932064, 2983680 with q = 4,7,8,9,13,14;

1556480, 3932040, 14008320 with q = 9,13,14;

106151936, 251650560, 955367424 with q = 9,13,14;

312792480, 1505806848 with q = 19,20;

6604416000, 30834059256 with q = 19,20;

9623577600, 46566269568 with q = 19,20.

When n is a term, we can test if s=sigma(n) is also a term; this way we get 6 more terms: 572941926400, 845734196736, 1422976331052, 4010593484800, 11383810648416, 36095341363200.

And the corresponding chains are:

173238912000, 845734196736 with q = 19,20;

355744082763, 572941926400, 1422976331052, 4010593484800, 11383810648416, 36095341363200 with q = 4,7,8,9,13,14. (End)

From Altug Alkan, May 17 2016: (Start)

Here are additional chains for the above list:

57120, 217728 with q = 13,14;

343976, 710400 with q = 7,8;

1980342, 5621760 with q = 10,14;

4404480, 14913024 with q = 11,12;

238608384, 775898880 with q = 11,12. (End)

Currently, the coefficient pairs are [1, 1], [3, 4], [4, 7], [7, 8], [8, 9], [9, 10], [9, 13], [10, 14], [11, 12], [13, 14], [16, 17], [16, 21], [17, 18], [19, 20], [23, 24], [25, 26], [25, 31], [27, 28], [29, 30], [31, 32], [32, 33], [37, 38]. It is interesting to note that for some of them, the pair (s,t) also satisfies t=sigma(s). - Michel Marcus, Jul 03 2016; Sep 06 2016

Using these empirical pairs of coefficients in conjunction with first comment allows us to find if some term is the sum of divisors of another yet unknown smaller term. - Michel Marcus, Jul 04 2016

LINKS

Jud McCranie, Table of n, a(n) for n = 1..130

G. L. Cohen and H. J. J. te Riele, Iterating the sum-of-divisors function, Experimental Mathematics, 5 (1996), pp. 93-100.

Michel Marcus, Unexhaustive list of terms, grouped by quotient q

MATHEMATICA

Select[Range[100000], Mod[DivisorSigma[1, DivisorSigma[1, #]], #] == 0 &] (* Carl Najafi, Aug 22 2011 *)

PROG

(PARI) is_A019278(n)=sigma(sigma(n))%n==0 \\ M. F. Hasler, Jul 02 2016

(Python)

from sympy.ntheory import divisor_sigma as D

print [i for i in range(1, 100000) if D(D(i, 1), 1)%i==0] # Indranil Ghosh, Mar 17 2017

CROSSREFS

Cf. A098219, A098220, A098221, A098222, A098223, A008333, A051027, A019276.

For sigma see A000203.

Sequence in context: A098056 A097100 A002954 * A267894 A084345 A084561

Adjacent sequences:  A019275 A019276 A019277 * A019279 A019280 A019281

KEYWORD

nonn

AUTHOR

N. J. A. Sloane

EXTENSIONS

Simpler definition from M. F. Hasler, Jul 02 2016

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 22 05:35 EDT 2019. Contains 326172 sequences. (Running on oeis4.)