cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A066523 Crowded numbers: for any k in the sequence, d(k)/k is larger than d(m)/m for all m > k.

Original entry on oeis.org

2, 4, 6, 12, 24, 30, 36, 48, 60, 72, 84, 120, 144, 180, 240, 252, 360, 420, 480, 504, 540, 720, 840, 900, 1008, 1080, 1260, 1440, 1680, 1800, 2520, 2640, 2880, 3360, 3780, 3960, 5040, 5280, 5400, 5460, 5544, 6300, 7560, 7920, 10080, 10920, 12600
Offset: 1

Views

Author

Roy Maulbogat (maulbogat(AT)gmail.com), Jan 05 2002

Keywords

Comments

Since d(m) < 2*sqrt(m), we need only test values of m < (2k/d(k))^2.
It was briefly conjectured that this sequence was the same as the highly composite numbers (A002182) larger than 1, but this is false: 30 is crowded but not highly composite and 50400 is highly composite but not crowded. Is every super-abundant number (A004394) crowded?
Additional comments from Roy Maulbogat, Jan 22 2008: (Start)
It can easily be shown that all crowded numbers are even and that there is always a crowded number between N and 2N. This allows us to improve the algorithm as follows:
crowded[k_] := Module[{},
* If[OddQ[k], Return [False]];*
div = DivisorSigma[0,k]/k;
For [ *m=k+2, m<=2k, m+=2*, If[
DivisorSigma [0, m] / m<=div, Return [False]]];True];
numlist = Select[Range[1,10^7],crowded]
On second thought, it might be wise to use Min[2k, stop] as the stopping condition of the loop ("stop" being the variable defined in the original algorithm). (End)

Crossrefs

Programs

  • Maple
    Tau := n -> NumberTheory:-tau(n): t := (n, k) -> k*Tau(n) < n*Tau(k):
    isCrowded := proc(k) local n; if k::odd then return false fi;
    andmap(j -> t(j, k), [seq(k+2..2*k,2)]) end:
    aList := n -> select(isCrowded, [seq(1..n)]): aList(1100);  # Peter Luschny, Jan 24 2025
  • Mathematica
    crowded[n_] := Module[{}, stop=(2/(dovern=DivisorSigma[0, n]/n))^2; For[m=n+1, m=dovern, Return[False]]]; True]; Select[Range[1, 13000], crowded]

Extensions

Edited by Dean Hickerson, Jan 07 2002