The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A019278 Numbers n such that sigma(sigma(n)) = k*n for some k. 37
 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 k = sigma(sigma(n))/n is an integer (see A098223). Then also k = (sigma(s)/s)*(sigma(n)/n) with s = sigma(n). That is, k = 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 k = 3,4,7,8,9,10; 512, 1023, 1536, 4092, 10752, 32736 with k = 3,4,7,8,9,10; 29127, 47360, 116508, 331520, 932064, 2983680 with k = 4,7,8,9,13,14; 1556480, 3932040, 14008320 with k = 9,13,14; 106151936, 251650560, 955367424 with k = 9,13,14; 312792480, 1505806848 with k = 19,20; 6604416000, 30834059256 with k = 19,20; 9623577600, 46566269568 with k = 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 k = 19,20; 355744082763, 572941926400, 1422976331052, 4010593484800, 11383810648416, 36095341363200 with k = 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 k = 13,14; 343976, 710400 with k = 7,8; 1980342, 5621760 with k = 10,14; 4404480, 14913024 with k = 11,12; 238608384, 775898880 with k = 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 For m in A090748 = A000043 - 1 and c in A205597 (= odd a(n)), c*2^m is in the sequence, unless 2^(m+1)-1 | sigma(c). Indeed, from sigma(xy) = sigma(x)*sigma(y) for gcd(x,y) = 1, we get sigma(sigma(c*2^m)) = sigma(sigma(c))*2^(m+1), so c*2^m is in the sequence if sigma(sigma(c))/c = k/2 (where k can't be odd: A330598 has no odd c). - M. F. Hasler, Jan 06 2020 LINKS Jud McCranie and Giovanni Resta, Table of n, a(n) for n = 1..145 (first 130 terms from Jud McCranie) 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 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, 10000) if D(D(i, 1), 1)%i==0]) # Indranil Ghosh, Mar 17 2017 (MAGMA) [m: m in [1..560000]| IsIntegral(DivisorSigma(1, DivisorSigma(1, m))/m)]; // Marius A. Burtea, Nov 16 2019 CROSSREFS Cf. A098219, A098220, A098221, A098222, A098223, A008333, A051027, A019276. For sigma see A000203 and A007691. Cf. A205597 (odd terms), A330598 (half-integer ratio). Sequence in context: A098056 A097100 A002954 * A267894 A084345 A084561 Adjacent sequences:  A019275 A019276 A019277 * A019279 A019280 A019281 KEYWORD nonn AUTHOR 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.

Last modified June 1 19:32 EDT 2020. Contains 334762 sequences. (Running on oeis4.)