A261100 -id:A261100 - cp's OEIS frontend

A002182 Highly composite numbers: numbers n where d(n), the number of divisors of n (A000005), increases to a record.

1, 2, 4, 6, 12, 24, 36, 48, 60, 120, 180, 240, 360, 720, 840, 1260, 1680, 2520, 5040, 7560, 10080, 15120, 20160, 25200, 27720, 45360, 50400, 55440, 83160, 110880, 166320, 221760, 277200, 332640, 498960, 554400, 665280, 720720, 1081080, 1441440, 2162160

Offset: 1

N. J. A. Sloane

Where record values of d(n) occur: d(n) > d(k) for all k < n.
A002183 is the RECORDS transform of A000005, i.e., lists the corresponding values d(n) for n in A002182.
Flammenkamp's page also has a copy of the missing Siano paper.
Highly composite numbers are the product of primorials, A002110. See A112779 for the number of primorial terms in the product of a highly composite number. - Jud McCranie, Jun 12 2005
Sigma and tau for highly composite numbers through the 146th entry conform to a power fit as follows: log(sigma)=A*log(tau)^B where (A,B) =~ (1.45,1.38). - Bill McEachen, May 24 2006
a(n) often corresponds to P(n,m) = number of permutations of n things taken m at a time. Specifically, if start=1, pointers 1-6, 9, 10, 13-15, 17-19, 22, 23, 28, 34, 37, 43, 52, ... An example is a(37)=665280, which is P(12,6)=12!/(12-6)!. - Bill McEachen, Feb 09 2009
Concerning the previous comment, if m=1, then P(n,m) can represent any number. So let's assume m > 1. Searching the first 1000 terms, the only indices of terms of the form P(n,m) are 4, 5, 6, 9, 10, 12, 13, 14, 15, 16, 17, 18, 19, 22, 23, 27, 28, 31, 34, 37, 41, 43, 44, 47, 50, 52, and 54. Note that a(44) = 4324320 = P(2079,2). See A163264. - T. D. Noe, Jun 10 2009
A large number of highly composite numbers have 9 as their digit root. - Parthasarathy Nambi, Jun 07 2009
Because 9 divides all highly composite numbers greater than 1680, those numbers have digital root 9. - T. D. Noe, Jul 24 2009
See A181309 for highly composite numbers that are not highly abundant.
a(n) is also defined by the recurrence: a(1) = 1, a(n+1)/sigma(a(n+1)) < a(n) / sigma(a(n)). - Michel Lagneau, Jan 02 2012 [NOTE: This "definition" is wrong (a(20)=7560 does not satisfy this inequality) and incomplete: It does not determine a sequence uniquely, e.g., any subsequence would satisfy the same relation. The intended meaning is probably the definition of the (different) sequence A004394. - M. F. Hasler, Sep 13 2012]
Up to a(1000), the terms beyond a(5) = 12 resp. beyond a(9) = 60 are a multiples of these. Is this true for all subsequent terms? - M. F. Hasler, Sep 13 2012 [Yes: see EXAMPLE in A199337! - M. F. Hasler, Jan 03 2020]
Differs from the superabundant numbers from a(20)=7560 on, which is not in A004394. The latter is not a subsequence of A002182, as might appear from considering the displayed terms: The two sequences have only 449 terms in common, the largest of which is A002182(2567) = A004394(1023). See A166735 for superabundant numbers that are not highly composite, and A004394 for further information. - M. F. Hasler, Sep 13 2012
Subset of A067128 and of A025487. - David A. Corneth, May 16 2016, Jan 03 2020
It seems that a(n) +- 1 is often prime. For n <= 1000 there are 210 individual primes and 17 pairs of twin primes. See link to Lim's paper below. - Dmitry Kamenetsky, Mar 02 2019
There are infinitely many numbers in this sequence and a(n+1) <= 2*a(n), because it is sufficient to multiply a(n) by 2 to get a number having more divisors. (This proves Guess 0 in the Lim paper.) For n = (1, 2, 4, 5, 9, 13, 18, ...) one has equality in this bound, but asymptotically a(n+1)/a(n) goes to 1, cf. formula due to Erdős. See A068507 for the terms such that a(n)+-1 are twin primes. - M. F. Hasler, Jun 23 2019
Conjecture: For n > 7, a(n) is a Zumkeller number (A083207). Verified for n up to and including 48. If this conjecture is true, one may base on it an alternative proof of the fact that for n>7 a(n) is not a perfect square (see Fact 5, Rao/Peng arXiv link at A083207). - Ivan N. Ianakiev, Jun 29 2019
The conjecture above is true (see the proof in the "Links" section). - Ivan N. Ianakiev, Jan 31 2020
The first instance of omega(a(n)) < omega(a(n-1)) (omega = A001221: number of prime divisors) is at a(26) = 45360. Up to n = 10^4, 1759 terms have this property, but omega decreases by 2 only at indices n = 5857, 5914 and 5971. - M. F. Hasler, Jan 02 2020
Inequality (54) in Ramanujan (1915) implies that for any m there is n* such that m | a(n) for all n > n*: see A199337 for the proof. - M. F. Hasler, Jan 03 2020

			a(5) = 12 is in the sequence because A000005(12) is larger than any earlier value in A000005. - _M. F. Hasler_, Jan 03 2020

CRC Press Standard Mathematical Tables, 28th Ed, p. 61.
J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 180, p. 56, Ellipses, Paris 2008.
L. E. Dickson, History of Theory of Numbers, I, p. 323.
Ross Honsberger, An introduction to Ramanujan's Highly Composite Numbers, Chap. 14 pp. 193-200 Mathematical Gems III, DME no. 9 MAA 1985
Jean-Louis Nicolas, On highly composite numbers, pp. 215-244 in Ramanujan Revisited, Editors G. E. Andrews et al., Academic Press 1988
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 88.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, 128.

Michael De Vlieger, Table of n, a(n) for n = 1..10000 (obtained from A. Flammenkamp's data; first 1000 terms from T. D. Noe)
Yu. Bilu, P. Habegger, and L. Kühne, Effective bounds for singular units, arXiv:1805.07167 [math.NT], 2018.
Benjamin Braun and Brian Davis, Antichain Simplices, arXiv:1901.01417 [math.CO], 2019.
Harold W. Ellingsen, Jr., A Fresh Look at Highly Composite Numbers, The American Mathematical Monthly, Vol. 126, No. 8 (2019), pp. 740-741.
Paul Erdős, On Highly composite numbers, J. London Math. Soc., Vol. 19 (1944), pp. 130-133, MR7,145d; Zentralblatt 61,79.
Achim Flammenkamp, Highly composite numbers.
Achim Flammenkamp, List of the first 1200 highly composite numbers.
Achim Flammenkamp, List of the first 779,674 highly composite numbers.
James Grime and Brady Haran, 5040 and other Anti-Prime Numbers, Numberphile video (2016).
Bob Hinman, Letter to N. J. A. Sloane, Aug. 1980.
Ivan N. Ianakiev, On the question "Which Highly composite numbers (A002182) are Zumkeller numbers (A083207)?".
Stepan Kochemazov, Oleg Zaikin, Eduard Vatutin, and Alexey Belyshev, Enumerating Diagonal Latin Squares of Order Up to 9, J. Int. Seq., Vol. 23 (2020), Article 20.1.2.
Aneesh M. Koya and P. P. Deepthi, Plug and play self-configurable IoT gateway node for telemonitoring of ECG, Computers in Biology and Medicine, Vol. 112 (2019), 103359.
Jeffrey C. Lagarias, An elementary problem equivalent to the Riemann hypothesis, Am. Math. Monthly 109 (#6, 2002), 534-543; arXiv:math/0008177 [math.NT], 2000-2001.
Bill Lauritzen, Versatile Numbers -Versatile Economics.
Benny Lim, Prime Numbers Generated From Highly Composite Numbers, Parabola Magazine, Volume 54, Issue 3, (2018).
R. J. Mathar, Maple program to convert the Flammenkamp file to an OEIS b-file.
R. J. Mathar, Output of above Maple program. [Uncompresses to 9.1 MB]
Graeme McRae, Highly Composite Numbers.
Jean-Louis Nicolas, Ordre maximal d'un élément du groupe S_n de permutations et 'highly composite numbers' (Text in French).
Jean-Louis Nicolas and Guy Robin, Highly Composite Numbers by Srinivasa Ramanujan, The Ramanujan Journal, Vol. 1(2), pp. 119-153, Kluwer Academics Pub.
Kevin O'Bryant, PlanetMath.org, Highly composite number.
S. Ramanujan, Highly composite numbers, Proceedings of the London Mathematical Society, Ser. 2, Vol. XIV, No. 1 (1915), pp. 347-409. (DOI: 10.1112/plms/s2_14.1.347, also available with an additional footnote in the PDF at http://ramanujan.sirinudi.org/Volumes/published/ram15.html)
Steven Ratering, An interesting subset of the highly composite numbers, Math. Mag., Vol. 64, No. 5 (1991), pp. 343-346.
Guy Robin, Méthodes d'optimisation pour un problème de théorie des nombres, RAIRO Informatique Théorique, Vol. 17, No. 3 (1983), pp. 239-247.
Vladimir Shevelev, On Erdős constant, arXiv:1605.08884 [math.NT], 2016.
D. B. Siano and J. D. Siano, An Algorithm for Generating Highly Composite Numbers, 1994.
N. J. A. Sloane, Transforms.
Michel Waldschmidt, From highly composite numbers to transcendental number theory, 2013.
Eric Weisstein's World of Mathematics, Highly Composite Number.
Wikipedia, Highly composite number.

Cf. A000005 (number of divisors), A002110, A002183, A002473, A004394, A025487, A106037, A108602, A112778, A112779, A112780, A112781, A006218, A126098, A002201, A072938, A094348, A003418, A161184, A037992 (infinitary analog), A108951, A329902, A352418.
Cf. A261100 (a left inverse).
Cf. A002808. - Peter J. Marko, Aug 16 2018
Cf. A279930 (highly composite and highly Brazilian).
Cf. A068507 (terms such that a(n)+-1 are twin primes).
Cf. A199337 (number of terms not divisible by n).

a = 0; Do[b = DivisorSigma[0, n]; If[b > a, a = b; Print[n]], {n, 1, 10^7}]
(* Convert A. Flammenkamp's 779674-term dataset; first, decompress, rename "HCN.txt": *)
a = Times @@ {Times @@ Prime@ Range@ ToExpression@ First@ #1, If[# == {}, 1, Times @@ MapIndexed[Prime[First@ #2]^#1 &, #]] &@ DeleteCases[-1 + Flatten@ Map[If[StringFreeQ[#, "^"], ToExpression@ #, ConstantArray[#1, #2] & @@ ToExpression@ StringSplit[#, "^"]] &, #2], 0]} & @@ TakeDrop[StringSplit@ #, 1] & /@ Import["HCN.txt", "Data"] (* Michael De Vlieger, May 08 2018 *)
DeleteDuplicates[Table[{n,DivisorSigma[0,n]},{n,2163000}],GreaterEqual[ #1[[2]],#2[[2]]]&] [[All,1]] (* Harvey P. Dale, May 13 2022 *)
NestList[Function[last,
  Module[{d = DivisorSigma[0, last]},
   NestWhile[# + 1 &, last, DivisorSigma[0, #] <= d &]]], 1, 40] (* Steven Lu, Mar 30 2023 *)

print1(r=1); forstep(n=2,1e5,2, if(numdiv(n)>r, r=numdiv(n); print1(", "n))) \\ Charles R Greathouse IV, Jun 10 2011

v002182 = [1]/*vector for memoization*/; A002182(n, i = #v002182) ={ if(n > i, v002182 = Vec(v002182, n); my(k = v002182[i], d, s=1); until(i == n, d = numdiv(k); s<60 && k>=60 && s=60; until(numdiv(k += s) > d,); v002182[i++] = k); k, v002182[n])} \\ Antti Karttunen, Jun 06 2017; edited by M. F. Hasler, Jan 03 2020 and Jun 20 2022

is_A002182(n, a=1, b=1)={while(n>A002182(b*=2), a*=2); until(a>b, my(m=(a+b)\2, t=A002182(m)); if(tn, b=m-1, return(m)))} \\ Also used in other sequences. - M. F. Hasler, Jun 20 2022

from sympy import divisor_count
A002182_list, r = [], 0
for i in range(1,10**4):
    d = divisor_count(i)
    if d > r:
        r = d
        A002182_list.append(i) # Chai Wah Wu, Mar 23 2015

Also, for n >= 2, smallest values of p for which A006218(p)-A006318(p-1) = A002183(n). - Philippe LALLOUET (philip.lallouet(AT)wanadoo.fr), Jun 23 2007
a(n+1) < a(n) * (1+log(a(n))^-c) for some positive c (see Erdős). - David A. Corneth, May 16 2016
a(n) = A108951(A329902(n)). - Antti Karttunen, Jan 08 2020
a(n+1) <= 2*a(n). For cases where the equal sign holds, see A072938. - A.H.M. Smeets, Jul 10 2021
Sum_{n>=1} 1/a(n) = A352418. - Amiram Eldar, Mar 24 2022

Jun 19 1996: Changed beginning to start at 1.
Jul 10 1996: Matthew Conroy points out that these are different from the super-abundant numbers - see A004394. Last 8 terms sent by J. Lowell; checked by Jud McCranie.
Description corrected by Gerard Schildberger and N. J. A. Sloane, Apr 04 2001
Additional references from Lekraj Beedassy, Jul 24 2001

A002183 Number of divisors of n-th highly composite number.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 30, 32, 36, 40, 48, 60, 64, 72, 80, 84, 90, 96, 100, 108, 120, 128, 144, 160, 168, 180, 192, 200, 216, 224, 240, 256, 288, 320, 336, 360, 384, 400, 432, 448, 480, 504, 512, 576, 600, 640, 672, 720, 768, 800, 864, 896

Offset: 1

N. J. A. Sloane

Record values of tau(n).
RECORDS transform of A000005.
All powers of 2 are present through 2^17. No power of 2 above that is present at least through 2^51. - Comment from Robert G. Wilson v, modified by Ray Chandler, Nov 10 2005
No power of 2 above 2^17 is contained in this sequence - see McRae link for proof. - Graeme McRae, Apr 27 2006
All numbers of the form 9*2^n are present for n=0 through n=30. - Richard Peterson, Sep 07 2024

S. Ramanujan, Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, p. 87.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
A. Flammenkamp, First 1200 highly composite numbers.
Graeme McRae, Highly Composite Numbers.
S. Ramanujan, Table of First 103 Highly Composite Numbers.
N. J. A. Sloane, Transforms.
Eric Weisstein's World of Mathematics, Highly Composite Number.

Cf. A000005, A002182, A002201, A006218, A061799, A070319, A243220, A261100, A329605, A329902.

import Data.List (nub)
a002183 n = a002183_list !! (n-1)
a002183_list = nub $ map (a000005 . a061799) [1..]
-- Reinhard Zumkeller, Apr 01 2011

Reap[ For[ record = 0; n = 1, n <= 10^9, n = If[n < 60, n+1, n+60], tau = DivisorSigma[0, n]; If[tau > record, record = tau; Print[tau]; Sow[tau]]]][[2, 1]] (* Jean-François Alcover, Aug 13 2013 *)
DeleteDuplicates[DivisorSigma[0,Range[3*10^6]],GreaterEqual] (* The program generates the first 42 terms of the sequence. *) (* Harvey P. Dale, Aug 12 2025 *)

a(n) = A000005(A002182(n)).
Also record values of differences A006218(p)-A006218(p-1). These record values occur for any p = A002182(q) where q>=2. - Philippe LALLOUET (philip.lallouet(AT)wanadoo.fr), Jun 23 2007
a(A261100(n)) = A070319(n). - Antti Karttunen, Jun 06 2017
a(n) = A329605(A329902(n)). - Antti Karttunen, Jan 14 2020

More terms from Robert G. Wilson v, Jul 24 2002

A060990 Number of solutions to x - d(x) = n, where d(n) is the number of divisors of n (A000005).

Original entry on oeis.org

2, 2, 1, 1, 1, 1, 3, 0, 0, 1, 1, 3, 1, 0, 1, 1, 1, 2, 1, 0, 0, 1, 4, 1, 0, 0, 1, 2, 0, 2, 1, 1, 1, 0, 2, 2, 0, 0, 2, 2, 0, 1, 1, 0, 1, 1, 3, 1, 2, 0, 0, 2, 0, 1, 1, 0, 0, 3, 2, 1, 1, 1, 2, 0, 0, 2, 0, 0, 0, 2, 4, 1, 1, 1, 0, 0, 1, 1, 2, 0, 1, 2, 1, 1, 1, 0, 1, 2, 0, 1, 1, 2, 1, 1, 1, 1, 2, 1, 0, 1, 0, 1, 3, 0, 1, 1

Offset: 0

Labos Elemer, May 11 2001

nonn

If x-d(x) is never equal to n, then n is in A045765 and a(n) = 0.

Number of solutions to A049820(x) = n. - Jaroslav Krizek, Feb 09 2014

			a(11) = 3 because three numbers satisfy equation x-d(x)=11, namely {13,15,16} with {2,4,5} divisors respectively.

Antti Karttunen, Table of n, a(n) for n = 0..110880

Cf. A000005, A002183, A049820, A049816, A082284, A155043, A236561, A236565, A259934, A261100, A262507, A262513, A262686.
Cf. A045765 (positions of zeros), A236562 (positions of nonzeros), A262511 (positions of ones).
Cf. A263087 (computed for squares).

lim = 105; s = Table[n - DivisorSigma[0, n], {n, 2 lim + 3}]; Length@ Position[s, #] & /@ Range[0, lim] (* Michael De Vlieger, Sep 29 2015, after Wesley Ivan Hurt at A049820 *)

allocatemem(123456789);
uplim = 2162160; \\ = A002182(41).
v060990 = vector(uplim);
for(n=3, uplim, v060990[n-numdiv(n)]++);
A060990 = n -> if(!n,2,v060990[n]);
uplim2 = 110880; \\ = A002182(30).
for(n=0, uplim2, write("b060990.txt", n, " ", A060990(n)));
\\ Antti Karttunen, Sep 25 2015

(define (A060990 n) (if (zero? n) 2 (add (lambda (k) (if (= (A049820 k) n) 1 0)) n (+ n (A002183 (+ 2 (A261100 n)))))))
;; Auxiliary function add implements sum_{i=lowlim..uplim} intfun(i)
(define (add intfun lowlim uplim) (let sumloop ((i lowlim) (res 0)) (cond ((> i uplim) res) (else (sumloop (1+ i) (+ res (intfun i)))))))
;; Proof-of-concept code for the given formula, by Antti Karttunen, Sep 25 2015

a(0) = 2; for n >= 1, a(n) = Sum_{k = n .. n+A002183(2+A261100(n))} [A049820(k) = n]. (Here [...] denotes the Iverson bracket, resulting 1 when A049820(k) is n and 0 otherwise.) - Antti Karttunen, Sep 25 2015, corrected Oct 12 2015.
a(n) = Sum_{k = A082284(n) .. A262686(n)} [A049820(k) = n] (when tacitly assuming that A049820(0) = 0.) - Antti Karttunen, Oct 12 2015
Other identities and observations. For all n >= 0:
a(A045765(n)) = 0. a(A236562(n)) > 0. - Jaroslav Krizek, Feb 09 2014

Offset corrected by Jaroslav Krizek, Feb 09 2014

A045765 k - d(k) never takes these values, where d(k) = A000005(k).

Original entry on oeis.org

7, 8, 13, 19, 20, 24, 25, 28, 33, 36, 37, 40, 43, 49, 50, 52, 55, 56, 63, 64, 66, 67, 68, 74, 75, 79, 85, 88, 98, 100, 103, 108, 109, 112, 113, 116, 117, 123, 124, 126, 131, 132, 133, 134, 136, 140, 143, 145, 150, 153, 156, 159, 160, 163, 164, 167, 168

Offset: 1

David W. Wilson

nonn

Complement of A236562. - Jaroslav Krizek, Feb 09 2014
Positions of zeros in A060990, leaf-nodes in the tree generated by edge-relation A049820(child) = parent. - Antti Karttunen, Oct 06 2015
Since A000005(x) <= 1 + x/2, k is in the sequence if there are no x <= 2*(k+1) with k = x - d(x). - Robert Israel, Oct 12 2015
This can be improved as: k is in the sequence if there are no x <= k + A002183(2+A261100(k)) with k = x - d(x). Cf. also A070319, A262686. - Antti Karttunen, Oct 12 2015
Luca (2005) proved that this seqeunce is infinite. - Amiram Eldar, Jul 26 2025

Antti Karttunen, Table of n, a(n) for n = 1..10000
Florian Luca, On numbers not of the form n-omega(n), Acta Mathematica Hungarica, Vol. 106, No. 1 (2005), pp. 117-135.

Top row of A262898.
Cf. A000005, A002183, A049820, A060990, A070319, A236562, A236565, A261100, A262511, A262686, A262901, A262902, A262903, A262909, A263081.
Cf. A263091 (primes in this sequence), A263095 (squares).
Cf. A259934 (gives the infinite trunk of the same tree, conjectured to be unique).

N:= 1000: # to get all terms <= N
sort(convert({$1..N} minus {seq(x - numtheory:-tau(x), x=1..2*(1+N))},list)); # Robert Israel, Oct 12 2015

lim = 10000; Take[Complement[Range@ lim, Sort@ DeleteDuplicates@ Table[n - DivisorSigma[0, n], {n, lim}]], 57] (* Michael De Vlieger, Oct 13 2015 *)

allocatemem((2^31)+(2^30));
uplim = 36756720 + 640; \\ = A002182(53) + A002183(53).
v060990 = vector(uplim);
for(n=3, uplim, v060990[n-numdiv(n)]++);
A060990 = n -> if(!n,2,v060990[n]);
uplim2 = 36756720;
n=0; k=1; while(n <= uplim2, if(0==A060990(n), write("b045765_big.txt", k, " ", n); k++); n++;);
\\ Antti Karttunen, Oct 09 2015

(define A045765 (ZERO-POS 1 1 A060990))
;; Using also IntSeq-library of Antti Karttunen, Oct 06 2015

A082284 a(n) = smallest number k such that k - tau(k) = n, or 0 if no such number exists, where tau(n) = the number of divisors of n (A000005).

Original entry on oeis.org

1, 3, 6, 5, 8, 7, 9, 0, 0, 11, 14, 13, 18, 0, 20, 17, 24, 19, 22, 0, 0, 23, 25, 27, 0, 0, 32, 29, 0, 31, 34, 35, 40, 0, 38, 37, 0, 0, 44, 41, 0, 43, 46, 0, 50, 47, 49, 51, 56, 0, 0, 53, 0, 57, 58, 0, 0, 59, 62, 61, 72, 65, 68, 0, 0, 67, 0, 0, 0, 71, 74, 73, 84, 77, 0, 0, 81, 79, 82, 0, 88

Offset: 0

Amarnath Murthy, Apr 14 2003

nonn

a(p-2) = p for odd primes p.

Antti Karttunen, Table of n, a(n) for n = 0..124340

Column 1 of A265751.
Cf. A000005, A002182, A002183, A049820, A060990, A261100.
Cf. A262686 (the largest such number), A262511 (positions where these are equal and nonzero).
Cf. A266114 (same sequence sorted into ascending order, with zeros removed).
Cf. A266115 (positive numbers missing from this sequence).
Cf. A266110 (number of iterations before zero is reached), A266116 (final nonzero value reached).
Cf. also tree A263267 and its illustration.

N:= 1000: # to get a(0) .. a(N)
V:= Array(0..N):
for k from 1 to 2*(N+1) do
  v:= k - numtheory:-tau(k);
  if v <= N and V[v] = 0 then V[v]:= k fi
od:
seq(V[n],n=0..N); # Robert Israel, Dec 21 2015

Table[k = 1; While[k - DivisorSigma[0, k] != n && k <= 2 (n + 1), k++]; If[k > 2 (n + 1), 0, k], {n, 0, 80}] (* Michael De Vlieger, Dec 22 2015 *)

allocatemem(123456789);
uplim1 = 2162160 + 320; \\ = A002182(41) + A002183(41).
uplim2 = 2162160;
v082284 = vector(uplim1);
A082284 = n -> if(!n,1,v082284[n]);
for(n=1, uplim1, k = n-numdiv(n); if((0 == A082284(k)), v082284[k] = n));
for(n=0, 124340, write("b082284.txt", n, " ", A082284(n)));
\\ Antti Karttunen, Dec 21 2015

(define (A082284 n) (if (zero? n) 1 (let ((u (+ n (A002183 (+ 2 (A261100 n)))))) (let loop ((k n)) (cond ((= (A049820 k) n) k) ((> k u) 0) (else (loop (+ 1 k))))))))
;; Antti Karttunen, Dec 21 2015

Other identities and observations. For all n >= 0:

a(n) <= A262686(n).

More terms from David Wasserman, Aug 31 2004

A262686 a(n) = largest number k such that k - d(k) = n, or 0 if no such number exists, where d(n) = the number of divisors of n (A000005).

Original entry on oeis.org

2, 4, 6, 5, 8, 7, 12, 0, 0, 11, 14, 16, 18, 0, 20, 17, 24, 21, 22, 0, 0, 23, 30, 27, 0, 0, 32, 36, 0, 33, 34, 35, 40, 0, 42, 39, 0, 0, 48, 45, 0, 43, 46, 0, 50, 47, 54, 51, 60, 0, 0, 55, 0, 57, 58, 0, 0, 64, 66, 61, 72, 65, 70, 0, 0, 69, 0, 0, 0, 75, 80, 73, 84, 77, 0, 0, 81, 79, 90, 0, 88, 85, 86, 87, 96, 0, 92, 91, 0, 93, 94, 100, 98, 99, 102, 97, 108, 105, 0, 101

Offset: 0

Antti Karttunen, Sep 28 2015

nonn

Antti Karttunen, Table of n, a(n) for n = 0..10080

Cf. A000005, A002183, A049820, A060990, A261100, A262512.

Cf. also A082284 (the smallest such number), A262511 (positions where these are equal and nonzero).

Table[k = 2 n + 3; While[Nor[k - DivisorSigma[0, k] == n, k == 0], k--]; k, {n, 0, 99}] (* Michael De Vlieger, Sep 29 2015 *)

(definec (A262686 n) (if (zero? n) 2 (let ((u (+ n (A002183 (+ 2 (A261100 n)))))) (let loop ((k u)) (cond ((= (A049820 k) n) k) ((< k n) 0) (else (loop (- k 1))))))))

A070319 a(n) = Max_{k=1..n} tau(k) where tau(x)=A000005(x) is the number of divisors of x.

Original entry on oeis.org

1, 2, 2, 3, 3, 4, 4, 4, 4, 4, 4, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12

Offset: 1

Benoit Cloitre, May 11 2002

Is this the same as A068509? - David Scambler, Sep 10 2012
They are different even asymptotically: A068509(n)=O(sqrt(n)), while a(n) does not have polynomial growth. One example where the sequences differ: a(625) = 24 < A068509(625). (The inequality is implied by the set {1,2,..,25} where each pair of the elements has lcm <= 625.) - Max Alekseyev, Sep 11 2012
The two sequences first differ when n = 336, due to the set of 21 elements {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 18, 21, 24, 30, 36, 42, 48} where each pair of elements has lcm <= 336, while no positive integer <= 336 has more than 20 divisors. Therefore A068509(336) = 21 and A070319(336) = 20. - William Rex Marshall, Sep 11 2012
Indices of records give A002182. - Omar E. Pol, Feb 18 2023

Sándor, J., Crstici, B., Mitrinović, Dragoslav S. Handbook of Number Theory I. Dordrecht: Kluwer Academic, 2006, p. 44.
S. Wigert, Sur l'ordre de grandeur du nombre des diviseurs d'un entier, Arkiv. for Math. 3 (1907), 1-9.

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
S. Ramanujan, Highly composite numbers, Proceedings of the London Mathematical Society, 2, XIV, 1915, 347 - 409.

Cf. A000005, A002182, A002183, A261100, A261104.

a070319 n = a070319_list !! (n-1)
a070319_list = scanl1 max $ map a000005 [1..]
-- Reinhard Zumkeller, Apr 01 2011

a = {0}; Do[AppendTo[a, Max[DivisorSigma[0, n], a[[n]]]], {n, 120}]; Rest@ a (* Michael De Vlieger, Sep 29 2015 *)

PARI
```
a(n)=vecmax(vector(n,k,numdiv(k)))
```

v=vector(100);v[1]=1;for(n=2,#v,v[n]=max(v[n-1],numdiv(n))); v \\ Charles R Greathouse IV, Sep 12 2012

A070319(n,m=1,s=2)={for(k=s,n,mM. F. Hasler, Sep 12 2012

{a=0;for(n=1,100,print1(a=A070319(n,a,n),","))} /* Using this pattern, computation of a(1..10^6) is faster than "normal" computation of a(1..3000). */

a(n) = exp(log(2) log(n) / log(log(n)) + O(log(n) log(log(log(n))) / (log(log(n)))^2)). (See Sándor reference for more formulas.) - Eric M. Schmidt, Jun 30 2013

a(n) = A002183(A261100(n)). - Antti Karttunen, Jun 06 2017

A181801 Number of divisors of n that are highly composite (A002182).

Original entry on oeis.org

1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 5, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 6, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 6, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 7, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 6, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 7, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 5, 1, 2, 1, 3, 1, 3, 1, 3, 1, 2, 1, 7, 1, 2, 1, 3, 1, 3, 1, 3, 1

Offset: 1

Matthew Vandermast, Nov 27 2010

nonn

A divisor d of integer n is highly composite iff more multiples of (n/d) divide n than divide any smaller positive integer. This is because the number of divisors of n that are multiples of (n/d) equals the number of divisors of d, or A000005(d). (Also see example.)

a(n) = a(n+12) if n is not a multiple of 12.

			6 is a multiple of 3 highly composite integers (1, 2 and 6); therefore a(6) = 3.
As the first comment implies, there are also a(6) = 3 values of m such that 6 sets a record for number of divisors that are multiples of m. These values of m are 1, 3 and 6. All four of 6's divisors are multiples of 1; two (namely, 3 and 6) are multiples of 3; and one (namely, 6) is a multiple of 6. Each of these totals exceeds the corresponding total for any positive integer smaller than 6.

Antti Karttunen, Table of n, a(n) for n = 1..20000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000
Eric Weisstein's World of Mathematics, Highly composite number.

Row n of A181802 gives highly composite divisors of n. Row n of A181803 gives values of m such that n sets a record for the number of its divisors that are multiples of m. Numbers that set records for a(n) are in A181806.
Cf. A002182, A181804, A181805, A322584, A352418.
Inverse Möbius transform of A322586.

v002182 = vector(128); v002182[1] = 1; \\ For memoization.
A002182(n) = { my(d,k); if(v002182[n],v002182[n], k = A002182(n-1); d = numdiv(k); while(numdiv(k) <= d, k=k+1); v002182[n] = k; k); };
A261100(n) = { my(k=1); while(A002182(k)<=n,k=k+1); (k-1); };
A322586(n) = if(1==n,1,(A261100(n)-A261100(n-1)));
A181801(n) = sumdiv(n,d,A322586(d)); \\ Antti Karttunen, Dec 20 2018

a(n) = Sum_{d|n} A322586(d). - Antti Karttunen, Dec 20 2018

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A352418 = 2.132872... . - Amiram Eldar, Jan 01 2024

cp's OEIS Frontend

A002182 Highly composite numbers: numbers n where d(n), the number of divisors of n (A000005), increases to a record.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Python

Formula

Extensions

A002183 Number of divisors of n-th highly composite number.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

Extensions

A060990 Number of solutions to x - d(x) = n, where d(n) is the number of divisors of n (A000005).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula

Extensions

A045765 k - d(k) never takes these values, where d(k) = A000005(k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Scheme

A082284 a(n) = smallest number k such that k - tau(k) = n, or 0 if no such number exists, where tau(n) = the number of divisors of n (A000005).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Scheme

Formula

Extensions

A262686 a(n) = largest number k such that k - d(k) = n, or 0 if no such number exists, where d(n) = the number of divisors of n (A000005).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Scheme