login
a(n) is the least k with sigma(k) >= n*k.
24

%I #35 Jan 21 2022 03:45:55

%S 1,6,120,27720,122522400,130429015516800,1970992304700453905270400,

%T 1897544233056092162003806758651798777216000,

%U 4368924363354820808981210203132513655327781713900627249499856876120704000

%N a(n) is the least k with sigma(k) >= n*k.

%C Following a suggestion from _Ed Pegg Jr_, the sequence can be written in a more readable form as: 1!, 3!, 5!, 11# * 3! * 2, 17# * 5! * 2, 29# * 7! * 4, 53# * 7! * 12, 89# * 11! * 2, 157# * 17# * 8! * 6, 271# * 23# * 10!, 487# * 29# * 10!, 857# * 37# * 11! * 42, 1487# * 53# * 15! * 2, ..., where p# = primorial(p) = A034386.

%C From _T. D. Noe_, Jul 06 2005: (Start)

%C Let c(p) be the smallest colossally-abundant number having the prime factor p. See A073751 for info about computing these numbers.

%C Then the terms of this sequence can be expressed as

%C a(2) = c(3)

%C a(3) = c(5) * 2

%C a(4) = c(11) / 2

%C a(5) = c(17) / 3

%C a(6) = c(29) * 14

%C a(7) = c(53)

%C a(8) = c(89) * 4

%C a(9) = c(157) * 34

%C a(10) = c(271) * 23

%C a(11) = c(487) / 2

%C a(12) = c(857) / 2

%C a(13) = c(1487) * 212

%C a(14) = c(2621) * 710

%C a(15) = c(4567) * 2/21

%C a(16) = c(8011) / 2

%C a(17) = c(13999) * 1630. (End)

%C Initially, each term is divisible by the previous one. Is there a reason this should always be true? - _Santi Spadaro_, Aug 13 2002

%C The conjecture a(n)|a(n+1) holds out to n=10. - Devin Kilminster (devin(AT)maths.uwa.edu.au), Mar 10 2003

%C The conjecture a(n)|a(n+1) fails for n=15. - _T. D. Noe_, Jul 08 2005

%C We have a(n) = min{A007539(n), A134716(n)}, and clearly A007539(n) != A134716(n) for every n. For what values of n is the former less than the latter? - _Jeppe Stig Nielsen_, Jun 16 2015

%H Walter Nissen, <a href="http://upforthecount.com/math/abundance.html">Abundancy : Some Resources </a>

%H T. D. Noe, <a href="http://www.sspectra.com/math/A023199.pdf">An algorithm for finding the least k with sigma(k) >= nk</a>

%o (PARI) a(n) = my(k=1); while (sigma(k)/k < n, k++); k; \\ _Michel Marcus_, Oct 07 2019

%Y A subsequence of A004394. The dominating primes are in A108402.

%K nonn

%O 1,2

%A _David W. Wilson_

%E More terms from _Walter Nissen_, Apr 15 1997

%E Further terms from Devin Kilminster (devin(AT)maths.uwa.edu.au), Mar 10 2003

%E The term a(10) = 271#23#10! was apparently found independently by _Bodo Zinser_ and _Don Reble_, circa Jul 05 2005

%E The next term, a(11) = 487#29#10!, was corrected by _Don Reble_, Jul 06 2005

%E a(12) = 857#37#11!42 from _Don Reble_, Jul 06 2005

%E a(13) = 1487#53#15!2 found by _T. D. Noe_ and confirmed by _Don Reble_, Jul 07 2005

%E a(14)-a(17) found by _T. D. Noe_ and rechecked by him Oct 11 2005

%E a(15) corrected. The conjecture still fails at n=15. - _T. D. Noe_, Oct 13 2009