login
A066961
Numbers k such that sigma(k) divides sigma(sigma(k)).
9
1, 5, 12, 54, 56, 87, 95, 276, 308, 427, 429, 446, 455, 501, 581, 611, 9120, 9180, 9504, 9720, 9960, 10296, 10620, 10740, 10824, 11070, 11310, 11480, 11484, 11556, 11628, 11748, 11934, 11960, 12024, 12036, 12072, 12084, 12376, 12460, 12510, 12570
OFFSET
1,2
COMMENTS
Is this sequence finite?
These are numbers k such that sigma(k) is a multiply-perfect number (A007691). - Ivan N. Ianakiev, Sep 13 2016
EXAMPLE
12 is in the sequence since sigma(12) = 1 + 2 + 3 + 4 + 6 + 12 = 28 divides sigma(28) = 1 + 2 + 4 + 7 + 14 + 28 = 56. - Michael B. Porter, Sep 22 2016
MAPLE
with(numtheory): A066961:=n->`if`(sigma(sigma(n)) mod sigma(n) = 0, n, NULL): seq(A066961(n), n=1..2*10^4); # Wesley Ivan Hurt, Sep 22 2016
MATHEMATICA
Select[Range[30000], Divisible[DivisorSigma[1, DivisorSigma[1, #]], DivisorSigma[1, #]] &] (* Ivan N. Ianakiev, Sep 13 2016 *)
PROG
(PARI) { n=0; for (m=1, 10^10, if (sigma(sigma(m)) % sigma(m) == 0, write("b066961.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, Apr 11 2010
(PARI) isok(n) = (s=sigma(n)) && ((sigma(s) % s) == 0); \\ Michel Marcus, Sep 17 2016
(Magma) [n: n in [1..13000] | (SumOfDivisors(SumOfDivisors(n)) mod SumOfDivisors(n) eq 0)]; // Vincenzo Librandi, Sep 13 2016
CROSSREFS
Subsequences: A323653 (intersection with A007691, or equally, with A019278), A353365 (where the quotient is a power of 2).
Sequence in context: A188118 A062978 A189421 * A179994 A131549 A342549
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Jan 26 2002
EXTENSIONS
More terms from Lior Manor, Feb 06 2002
STATUS
approved