|
%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
| 