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.

Showing 1-10 of 47 results. Next

A000705 n-th superior highly composite number A002201(n) is product of first n terms of this sequence.

Original entry on oeis.org

2, 3, 2, 5, 2, 3, 7, 2, 11, 13, 2, 3, 5, 17, 19, 2, 23, 7, 29, 3, 31, 2, 37, 41, 43, 47, 5, 53, 59, 2, 11, 61, 3, 67, 71, 73, 79, 13, 83, 89, 2, 97, 101, 103, 107, 7, 109, 113, 17, 127, 131, 137, 139, 3, 5, 149, 151, 19, 2, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199
Offset: 1

Views

Author

Keywords

Comments

The Mathematica program uses the fact that the ratio of consecutive superior highly composite numbers is a prime, which was proved by Ramanujan. Ramanujan computed the first 50 terms of this sequence. Related sequences are A004490 and A073751, having to do with colossally abundant numbers.

References

  • S. Ramanujan, Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, p. 115.
  • S. Ramanujan, Ramanujan's Papers, pp. 147-9, Ed. B. J. Venkatachala et al., Prism Books, Bangalore 2000.
  • 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).

Crossrefs

Programs

  • Mathematica
    pFactor[f_List] := Module[{p = f[[1]], k = f[[2]]}, N[Log[(k + 2)/(k + 1)]/Log[p]]]; maxN = 100; f = {{2, 1}, {3, 0}}; primes = 1; lst = {2}; x = Table[pFactor[f[[i]]], {i, primes + 1}]; For[n = 2, n <= maxN, n++, i = Position[x, Max[x]][[1, 1]]; AppendTo[lst, f[[i, 1]]]; f[[i, 2]]++; If[i > primes, primes++; AppendTo[f, {Prime[i + 1], 0}]; AppendTo[x, pFactor[f[[ -1]]]]]; x[[i]] = pFactor[f[[i]]]]; lst (* T. D. Noe, Nov 01 2002 *)

Extensions

Edited by T. D. Noe, Nov 01 2002

A301416 Numbers k in A301413 such that k * A002110 (m) is in A002201.

Original entry on oeis.org

1, 2, 4, 12, 24, 48, 144, 720, 1440, 10080, 30240, 60480, 302400, 604800, 6652800, 19958400, 259459200, 518918400, 3632428800, 61751289600, 185253868800, 926269344000, 17599117536000, 35198235072000, 809559406656000, 1619118813312000, 4857356439936000
Offset: 1

Views

Author

Michael De Vlieger, Apr 09 2018

Keywords

Examples

			From _Michael De Vlieger_, May 14 2018: (Start)
Numbers m = A301416(x) * A002110(y) that are in A002201 are plotted below.
      1     2     3     4      5        6        7 ...
  +-----------------------------------------------
1 |   2
2 |   6    12
3 |        60   120   360
4 |                  2520   5040
5 |                        55440
6 |                       720720  1441440  4324320
...
(End)
		

Crossrefs

Programs

  • Mathematica
    t = Import["b002201.txt", "Data"][[All, -1]]; (* Uses b-file at A002201 Alternatively, use this conversion of terms at A000705 to a 10^4 term dataset for A002201. Processing 10^4 terms will take a long time: *) t = With[{s = Import["b000705.txt", "Data"][[All, -1]]}, FoldList[Times, s]]; f[n_] := With[{d = FactorInteger@ n}, If[n == 1, {0}, ReplacePart[Table[0, {PrimePi[d[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, d]]]; Union@ Array[Times @@ MapIndexed[Prime[First@#2]^#1 &, #3] & @@ {#1, Boole[First@ #2 > 0] Length@ #2, DeleteCases[-1 + #2, 0] /. -1 -> 0} & @@ {t[[#]], f@ t[[#]]} &, Length@ t]

A333655 Highly composite numbers (A002182) that are not superior highly composite numbers (A002201).

Original entry on oeis.org

1, 4, 24, 36, 48, 180, 240, 720, 840, 1260, 1680, 7560, 10080, 15120, 20160, 25200, 27720, 45360, 50400, 83160, 110880, 166320, 221760, 277200, 332640, 498960, 554400, 665280, 1081080, 2162160, 2882880, 3603600, 6486480, 7207200, 8648640, 10810800, 14414400
Offset: 1

Views

Author

Iain Fox, Aug 23 2020

Keywords

Comments

For a number n to be in this sequence, it must have the following conditions be true, where d(n) represents the number of divisors of n (A000005): d(n) > d(k), for all k < n, and there does not exist a number e > 0 such that d(n)/n^e >= d(k)/k^e for k < n and d(n)/n^e > d(k)/k^e for k > n.
This sequence is the same as A189228 until n=12, for which a(12) = 7560 and A189228(12) = 10080.

Examples

			4 is in the sequence because it has three factors, more than any preceding number, making it highly composite, but it is not a superior highly composite number.
		

Crossrefs

Highly composite numbers: A002182.
Superior highly composite numbers: A002201.

Programs

  • PARI
    lista(nn)=my(v, w=[1,2,4], r=1, p=primes(primepi(2^log(nn)))); v=setminus(Set(vector(nn, i, prod(n=1, primepi(2^log(i)), p[n]^floor(1/(p[n]^(1/log(i))-1))))), [1]); forstep(x=6, v[#v], 6, if(numdiv(x)>r, r=numdiv(x); w=setunion(w, [x]))); setminus(w, v)

A328549 1, together with the numbers that are simultaneously superior highly composite (A002201), colossally abundant (A004490), deeply composite (A095848), and miserable average divisor numbers (A263572).

Original entry on oeis.org

1, 2, 6, 12, 60, 120, 360, 2520, 5040, 55440, 720720, 1441440
Offset: 1

Views

Author

N. J. A. Sloane, Oct 20 2019

Keywords

Comments

Presumably there are no further terms.
From Hal M. Switkay, Nov 04 2019: (Start)
1. a(n+1) is the product of the first n terms of A328852.
2. This sequence is most rapidly constructed as the intersection of A095849 and A224078. It is designed to list all potential solutions to a question. Let n be a natural number, k real <= 0, e real > 0. Let P(n,k,e) state: on the domain of natural numbers, sigma_k(x)/x^e reaches a maximum at x = n. This implies Q(n,k): sigma_k(n) > sigma_k(m) for m < n a natural number. We ask: for which natural numbers n is it true for all real k <= 0 that there is a real e > 0 such that P(n,k,e)?
If any such n exist, they must belong to the present sequence. A095849 consists of all natural numbers n such that for all real k <= 0, Q(n,k) holds. A224078 consists of all natural numbers n such that for some real e0 and e1 both > 0, P(n,0,e0) and P(n,-1,e1) hold. It would be interesting to see the list of n for which there is an e2 > 0 such that P(n,-2,e2) holds.
Conjecture: the solutions to this problem, if any, form an initial sequence of the present sequence. (End)
Every term of this sequence is also in A065385: a record for the cototient function. - Hal M. Switkay, Feb 27 2021
Every term of this sequence, except the first, is also in A210594: factor-dense numbers. - Hal M. Switkay, Mar 29 2021

References

Crossrefs

1 together with the intersection of A002201, A004490, A095848, A263572.

A002110 Primorial numbers (first definition): product of first n primes. Sometimes written prime(n)#.

Original entry on oeis.org

1, 2, 6, 30, 210, 2310, 30030, 510510, 9699690, 223092870, 6469693230, 200560490130, 7420738134810, 304250263527210, 13082761331670030, 614889782588491410, 32589158477190044730, 1922760350154212639070, 117288381359406970983270, 7858321551080267055879090
Offset: 0

Views

Author

Keywords

Comments

See A034386 for the second definition of primorial numbers: product of primes in the range 2 to n.
a(n) is the least number N with n distinct prime factors (i.e., omega(N) = n, cf. A001221). - Lekraj Beedassy, Feb 15 2002
Phi(n)/n is a new minimum for each primorial. - Robert G. Wilson v, Jan 10 2004
Smallest number stroked off n times after the n-th sifting process in an Eratosthenes sieve. - Lekraj Beedassy, Mar 31 2005
Apparently each term is a new minimum for phi(x)*sigma(x)/x^2. 6/Pi^2 < sigma(x)*phi(x)/x^2 < 1 for n > 1. - Jud McCranie, Jun 11 2005
Let f be a multiplicative function with f(p) > f(p^k) > 1 (p prime, k > 1), f(p) > f(q) > 1 (p, q prime, p < q). Then the record maxima of f occur at n# for n >= 1. Similarly, if 0 < f(p) < f(p^k) < 1 (p prime, k > 1), 0 < f(p) < f(q) < 1 (p, q prime, p < q), then the record minima of f occur at n# for n >= 1. - David W. Wilson, Oct 23 2006
Wolfe and Hirshberg give ?, ?, ?, ?, ?, 30030, ?, ... as a puzzle.
Records in number of distinct prime divisors. - Artur Jasinski, Apr 06 2008
For n >= 2, the digital roots of a(n) are multiples of 3. - Parthasarathy Nambi, Aug 19 2009 [with corrections by Zak Seidov, Aug 30 2015]
Denominators of the sum of the ratios of consecutive primes (see A094661). - Vladimir Joseph Stephan Orlovsky, Oct 24 2009
Where record values occur in A001221. - Melinda Trang (mewithlinda(AT)yahoo.com), Apr 15 2010
It can be proved that there are at least T prime numbers less than N, where the recursive function T is: T = N - N*Sum_{i = 0..T(sqrt(N))} A005867(i)/A002110(i). This can show for example that at least 0.16*N numbers are primes less than N for 29^2 > N > 23^2. - Ben Paul Thurston, Aug 23 2010
The above comment from Parthasarathy Nambi follows from the observation that digit summing produces a congruent number mod 9, so the digital root of any multiple of 3 is a multiple of 3. prime(n)# is divisible by 3 for n >= 2. - Christian Schulz, Oct 30 2013
The peaks (i.e., local maximums) in a graph of the number of repetitions (i.e., the tally of values) vs. value, as generated by taking the differences of all distinct pairs of odd prime numbers within a contiguous range occur at regular periodic intervals given by the primorial numbers 6 and greater. Larger primorials yield larger (relative) peaks, however the range must be >50% larger than the primorial to be easily observed. Secondary peaks occur at intervals of those "near-primorials" divisible by 6 (e.g., 42). See A259629. Also, periodicity at intervals of 6 and 30 can be observed in the local peaks of all possible sums of two, three or more distinct odd primes within modest contiguous ranges starting from p(2) = 3. - Richard R. Forberg, Jul 01 2015
If a number k and a(n) are coprime and k < (prime(n+1))^b < a(n), where b is an integer, then k has fewer than b prime factors, counting multiplicity (i.e., bigomega(k) < b, cf. A001222). - Isaac Saffold, Dec 03 2017
If n > 0, then a(n) has 2^n unitary divisors (A034444), and a(n) is a record; i.e., if k < a(n) then k has fewer unitary divisors than a(n) has. - Clark Kimberling, Jun 26 2018
Unitary superabundant numbers: numbers k with a record value of the unitary abundancy index, A034448(k)/k > A034448(m)/m for all m < k. - Amiram Eldar, Apr 20 2019
Psi(n)/n is a new maximum for each primorial (psi = A001615) [proof in link: Patrick Sole and Michel Planat, proposition 1 page 2]; compare with comment 2004: Phi(n)/n is a new minimum for each primorial. - Bernard Schott, May 21 2020
The term "primorial" was coined by Harvey Dubner (1987). - Amiram Eldar, Apr 16 2021
a(n)^(1/n) is approximately (n log n)/e. - Charles R Greathouse IV, Jan 03 2023
Subsequence of A267124. - Frank M Jackson, Apr 14 2023

Examples

			a(9) = 23# = 2*3*5*7*11*13*17*19*23 = 223092870 divides the difference 5283234035979900 in the arithmetic progression of 26 primes A204189. - _Jonathan Sondow_, Jan 15 2012
		

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 50.
  • G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, Cambridge, University Press, 1940, p. 49.
  • P. Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 4.
  • 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 114.
  • D. Wolfe and S. Hirshberg, Underspecified puzzles, in Tribute to A Mathemagician, Peters, 2005, pp. 73-74.

Crossrefs

A034386 gives the second version of the primorial numbers.
Subsequence of A005117 and of A064807. Apart from the first term, a subsequence of A083207.
Cf. A001615, A002182, A002201, A003418, A005235, A006862, A034444 (unitary divisors), A034448, A034387, A033188, A035345, A035346, A036691 (compositorial numbers), A049345 (primorial base representation), A057588, A060735 (and integer multiples), A061742 (squares), A072938, A079266, A087315, A094348, A106037, A121572, A053589, A064648, A132120, A260188.
Cf. A061720 (first differences), A143293 (partial sums).
Cf. also A276085, A276086.
The following fractions are all related to each other: Sum 1/n: A001008/A002805, Sum 1/prime(n): A024451/A002110 and A106830/A034386, Sum 1/nonprime(n): A282511/A282512, Sum 1/composite(n): A250133/A296358.

Programs

  • Haskell
    a002110 n = product $ take n a000040_list
    a002110_list = scanl (*) 1 a000040_list
    -- Reinhard Zumkeller, Feb 19 2012, May 03 2011
    
  • Magma
    [1] cat [&*[NthPrime(i): i in [1..n]]: n in [1..20]]; // Bruno Berselli, Oct 24 2012
    
  • Magma
    [1] cat [&*PrimesUpTo(p): p in PrimesUpTo(60)]; // Bruno Berselli, Feb 08 2015
    
  • Maple
    A002110 := n -> mul(ithprime(i),i=1..n);
  • Mathematica
    FoldList[Times, 1, Prime[Range[20]]]
    primorial[n_] := Product[Prime[i], {i, n}]; Array[primorial,20] (* José María Grau Ribas, Feb 15 2010 *)
    Join[{1}, Denominator[Accumulate[1/Prime[Range[20]]]]] (* Harvey P. Dale, Apr 11 2012 *)
  • PARI
    a(n)=prod(i=1,n, prime(i)) \\ Washington Bomfim, Sep 23 2008
    
  • PARI
    p=1; for (n=0, 100, if (n, p*=prime(n)); write("b002110.txt", n, " ", p) )  \\ Harry J. Smith, Nov 13 2009
    
  • PARI
    a(n) = factorback(primes(n)) \\ David A. Corneth, May 06 2018
    
  • Python
    from sympy import primorial
    def a(n): return 1 if n < 1 else primorial(n)
    [a(n) for n in range(51)]  # Indranil Ghosh, Mar 29 2017
    
  • Sage
    [sloane.A002110(n) for n in (1..20)] # Giuseppe Coppoletta, Dec 05 2014
    
  • Scheme
    ; with memoization-macro definec
    (definec (A002110 n) (if (zero? n) 1 (* (A000040 n) (A002110 (- n 1))))) ;; Antti Karttunen, Aug 30 2016

Formula

Asymptotic expression for a(n): exp((1 + o(1)) * n * log(n)) where o(1) is the "little o" notation. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 08 2001
a(n) = A054842(A002275(n)).
Binomial transform = A136104: (1, 3, 11, 55, 375, 3731, ...). Equals binomial transform of A121572: (1, 1, 3, 17, 119, 1509, ...). - Gary W. Adamson, Dec 14 2007
a(0) = 1, a(n+1) = prime(n)*a(n). - Juri-Stepan Gerasimov, Oct 15 2010
a(n) = Product_{i=1..n} A000040(i). - Jonathan Vos Post, Jul 17 2008
a(A051838(n)) = A116536(n) * A007504(A051838(n)). - Reinhard Zumkeller, Oct 03 2011
A000005(a(n)) = 2^n. - Carlos Eduardo Olivieri, Jun 16 2015
a(n) = A035345(n) - A005235(n) for n > 0. - Jonathan Sondow, Dec 02 2015
For all n >= 0, a(n) = A276085(A000040(n+1)), a(n+1) = A276086(A143293(n)). - Antti Karttunen, Aug 30 2016
A054841(a(n)) = A002275(n). - Michael De Vlieger, Aug 31 2016
a(n) = A270592(2*n+2) - A270592(2*n+1) if 0 <= n <= 4 (conjectured for all n by Alon Kellner). - Jonathan Sondow, Mar 25 2018
Sum_{n>=1} 1/a(n) = A064648. - Amiram Eldar, Oct 16 2020
Sum_{n>=1} (-1)^(n+1)/a(n) = A132120. - Amiram Eldar, Apr 12 2021
Theta being Chebyshev's theta function, a(0) = exp(theta(1)), and for n > 0, a(n) = exp(theta(m)) for A000040(n) <= m < A000040(n+1) where m is an integer. - Miles Englezou, Nov 26 2024

A025487 Least integer of each prime signature A124832; also products of primorial numbers A002110.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 72, 96, 120, 128, 144, 180, 192, 210, 216, 240, 256, 288, 360, 384, 420, 432, 480, 512, 576, 720, 768, 840, 864, 900, 960, 1024, 1080, 1152, 1260, 1296, 1440, 1536, 1680, 1728, 1800, 1920, 2048, 2160, 2304, 2310
Offset: 1

Views

Author

Keywords

Comments

All numbers of the form 2^k1*3^k2*...*p_n^k_n, where k1 >= k2 >= ... >= k_n, sorted.
A111059 is a subsequence. - Reinhard Zumkeller, Jul 05 2010
Choie et al. (2007) call these "Hardy-Ramanujan integers". - Jean-François Alcover, Aug 14 2014
The exponents k1, k2, ... can be read off Abramowitz & Stegun p. 831, column labeled "pi".
For all such sequences b for which it holds that b(n) = b(A046523(n)), the sequence which gives the indices of records in b is a subsequence of this sequence. For example, A002182 which gives the indices of records for A000005, A002110 which gives them for A001221 and A000079 which gives them for A001222. - Antti Karttunen, Jan 18 2019
The prime signature corresponding to a(n) is given in row n of A124832. - M. F. Hasler, Jul 17 2019

Examples

			The first few terms are 1, 2, 2^2, 2*3, 2^3, 2^2*3, 2^4, 2^3*3, 2*3*5, ...
		

Crossrefs

Subsequence of A055932, A191743, and A324583.
Cf. A085089, A101296 (left inverses).
Equals range of values taken by A046523.
Cf. A178799 (first differences), A247451 (squarefree kernel), A146288 (number of divisors).
Rearrangements of this sequence include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181821, A181822, A322827, A329886, A329887.
Cf. also array A124832 (row n = prime signature of a(n)) and A304886, A307056.

Programs

  • Haskell
    import Data.Set (singleton, fromList, deleteFindMin, union)
    a025487 n = a025487_list !! (n-1)
    a025487_list = 1 : h [b] (singleton b) bs where
       (_ : b : bs) = a002110_list
       h cs s xs'@(x:xs)
         | m <= x    = m : h (m:cs) (s' `union` fromList (map (* m) cs)) xs'
         | otherwise = x : h (x:cs) (s  `union` fromList (map (* x) (x:cs))) xs
         where (m, s') = deleteFindMin s
    -- Reinhard Zumkeller, Apr 06 2013
    
  • Maple
    isA025487 := proc(n)
        local pset,omega ;
        pset := sort(convert(numtheory[factorset](n),list)) ;
        omega := nops(pset) ;
        if op(-1,pset) <> ithprime(omega) then
            return false;
        end if;
        for i from 1 to omega-1 do
            if padic[ordp](n,ithprime(i)) < padic[ordp](n,ithprime(i+1)) then
                return false;
            end if;
        end do:
        true ;
    end proc:
    A025487 := proc(n)
        option remember ;
        local a;
        if n = 1 then
            1 ;
        else
            for a from procname(n-1)+1 do
                if isA025487(a) then
                    return a;
                end if;
            end do:
        end if;
    end proc:
    seq(A025487(n),n=1..100) ; # R. J. Mathar, May 25 2017
  • Mathematica
    PrimeExponents[n_] := Last /@ FactorInteger[n]; lpe = {}; ln = {1}; Do[pe = Sort@PrimeExponents@n; If[ FreeQ[lpe, pe], AppendTo[lpe, pe]; AppendTo[ln, n]], {n, 2, 2350}]; ln (* Robert G. Wilson v, Aug 14 2004 *)
    (* Second program: generate all terms m <= A002110(n): *)
    f[n_] := {{1}}~Join~
      Block[{lim = Product[Prime@ i, {i, n}],
       ww = NestList[Append[#, 1] &, {1}, n - 1], dec},
       dec[x_] := Apply[Times, MapIndexed[Prime[First@ #2]^#1 &, x]];
       Map[Block[{w = #, k = 1},
          Sort@ Prepend[If[Length@ # == 0, #, #[[1]]],
            Product[Prime@ i, {i, Length@ w}] ] &@ Reap[
             Do[
              If[# < lim,
                 Sow[#]; k = 1,
                 If[k >= Length@ w, Break[], k++]] &@ dec@ Set[w,
                 If[k == 1,
                   MapAt[# + 1 &, w, k],
                   PadLeft[#, Length@ w, First@ #] &@
                     Drop[MapAt[# + Boole[i > 1] &, w, k], k - 1] ]],
               {i, Infinity}] ][[-1]]
    ] &, ww]]; Sort[Join @@ f@ 13] (* Michael De Vlieger, May 19 2018 *)
  • PARI
    isA025487(n)=my(k=valuation(n,2),t);n>>=k;forprime(p=3,default(primelimit),t=valuation(n,p);if(t>k,return(0),k=t);if(k,n/=p^k,return(n==1))) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    factfollow(n)={local(fm, np, n2);
      fm=factor(n); np=matsize(fm)[1];
      if(np==0,return([2]));
      n2=n*nextprime(fm[np,1]+1);
      if(np==1||fm[np,2]Franklin T. Adams-Watters, Dec 01 2011 */
    
  • PARI
    is(n) = {if(n==1, return(1)); my(f = factor(n));  f[#f~, 1] == prime(#f~) && vecsort(f[, 2],,4) == f[, 2]} \\ David A. Corneth, Feb 14 2019
    
  • PARI
    upto(Nmax)=vecsort(concat(vector(logint(Nmax,2),n,select(t->t<=Nmax,if(n>1,[factorback(primes(#p),Vecrev(p)) || p<-partitions(n)],[1,2]))))) \\ M. F. Hasler, Jul 17 2019
    
  • PARI
    \\ For fast generation of large number of terms, use this program:
    A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)}; \\ From A283980
    A025487list(e) = { my(lista = List([1, 2]), i=2, u = 2^e, t); while(lista[i] != u, if(2*lista[i] <= u, listput(lista,2*lista[i]); t = A283980(lista[i]); if(t <= u, listput(lista,t))); i++); vecsort(Vec(lista)); }; \\ Returns a list of terms up to the term 2^e.
    v025487 = A025487list(101);
    A025487(n) = v025487[n];
    for(n=1,#v025487,print1(A025487(n), ", ")); \\ Antti Karttunen, Dec 24 2019
    
  • Sage
    def sharp_primorial(n): return sloane.A002110(prime_pi(n))
    N = 2310
    nmax = 2^floor(log(N,2))
    sorted([j for j in (prod(sharp_primorial(t[0])^t[1] for k, t in enumerate(factor(n))) for n in (1..nmax)) if j <= N])
    # Giuseppe Coppoletta, Jan 26 2015

Formula

What can be said about the asymptotic behavior of this sequence? - Franklin T. Adams-Watters, Jan 06 2010
Hardy & Ramanujan prove that there are exp((2 Pi + o(1))/sqrt(3) * sqrt(log x/log log x)) members of this sequence up to x. - Charles R Greathouse IV, Dec 05 2012
From Antti Karttunen, Jan 18 & Dec 24 2019: (Start)
A085089(a(n)) = n.
A101296(a(n)) = n [which is the first occurrence of n in A101296, and thus also a record.]
A001221(a(n)) = A061395(a(n)) = A061394(n).
A007814(a(n)) = A051903(a(n)) = A051282(n).
a(A101296(n)) = A046523(n).
a(A306802(n)) = A002182(n).
a(n) = A108951(A181815(n)) = A329900(A181817(n)).
If A181815(n) is odd, a(n) = A283980(a(A329904(n))), otherwise a(n) = 2*a(A329904(n)).
(End)
Sum_{n>=1} 1/a(n) = Product_{n>=1} 1/(1 - 1/A002110(n)) = A161360. - Amiram Eldar, Oct 20 2020

Extensions

Offset corrected by Matthew Vandermast, Oct 19 2008
Minor correction by Charles R Greathouse IV, Sep 03 2010

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

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

References

  • 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.

Crossrefs

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).

Programs

  • Mathematica
    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 *)
  • PARI
    print1(r=1); forstep(n=2,1e5,2, if(numdiv(n)>r, r=numdiv(n); print1(", "n))) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    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
    
  • PARI
    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
    
  • Python
    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

Formula

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

Extensions

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

A003418 Least common multiple (or LCM) of {1, 2, ..., n} for n >= 1, a(0) = 1.

Original entry on oeis.org

1, 1, 2, 6, 12, 60, 60, 420, 840, 2520, 2520, 27720, 27720, 360360, 360360, 360360, 720720, 12252240, 12252240, 232792560, 232792560, 232792560, 232792560, 5354228880, 5354228880, 26771144400, 26771144400, 80313433200, 80313433200, 2329089562800, 2329089562800
Offset: 0

Views

Author

Roland Anderson (roland.anderson(AT)swipnet.se)

Keywords

Comments

The minimal exponent of the symmetric group S_n, i.e., the least positive integer for which x^a(n)=1 for all x in S_n. - Franz Vrabec, Dec 28 2008
Product over all primes of highest power of prime less than or equal to n. a(0) = 1 by convention.
Also smallest number whose set of divisors contains an n-term arithmetic progression. - Reinhard Zumkeller, Dec 09 2002
An assertion equivalent to the Riemann hypothesis is: | log(a(n)) - n | < sqrt(n) * log(n)^2. - Lekraj Beedassy, Aug 27 2006. (This is wrong for n = 1 and n = 2. Should "for n large enough" be added? - Georgi Guninski, Oct 22 2011)
Corollary 3 of Farhi gives a proof that a(n) >= 2^(n-1). - Jonathan Vos Post, Jun 15 2009
Appears to be row products of the triangle T(n,k) = b(A010766) where b = A130087/A130086. - Mats Granvik, Jul 08 2009
Greg Martin (see link) proved that "the product of the Gamma function sampled over the set of all rational numbers in the open interval (0,1) whose denominator in lowest terms is at most n" equals (2*Pi)^(1/2)*a(n)^(-1/2). - Jonathan Vos Post, Jul 28 2009
a(n) = lcm(A188666(n), A188666(n)+1, ..., n). - Reinhard Zumkeller, Apr 25 2011
a(n+1) is the smallest integer such that all polynomials a(n+1)*(1^i + 2^i + ... + m^i) in m, for i=0,1,...,n, are polynomials with integer coefficients. - Vladimir Shevelev, Dec 23 2011
It appears that A020500(n) = a(n)/a(n-1). - Asher Auel, corrected by Bill McEachen, Apr 05 2024
n-th distinct value = A051451(n). - Matthew Vandermast, Nov 27 2009
a(n+1) = least common multiple of n-th row in A213999. - Reinhard Zumkeller, Jul 03 2012
For n > 2, (n-1) = Sum_{k=2..n} exp(a(n)*2*i*Pi/k). - Eric Desbiaux, Sep 13 2012
First column minus second column of A027446. - Eric Desbiaux, Mar 29 2013
For n > 0, a(n) is the smallest number k such that n is the n-th divisor of k. - Michel Lagneau, Apr 24 2014
Slowest growing integer > 0 in Z converging to 0 in Z^ when considered as profinite integer. - Herbert Eberle, May 01 2016
What is the largest number of consecutive terms that are all equal? I found 112 equal terms from a(370261) to a(370372). - Dmitry Kamenetsky, May 05 2019
Answer: there exist arbitrarily long sequences of consecutive terms with the same value; also, the maximal run of consecutive terms with different values is 5 from a(1) to a(5) (see link Roger B. Eggleton). - Bernard Schott, Aug 07 2019
Related to the inequality (54) in Ramanujan's paper about highly composite numbers A002182, also used in A199337: a(A329570(m))^2 is a (not minimal) bound above which all highly composite numbers are divisible by m, according to the right part of that inequality. - M. F. Hasler, Jan 04 2020
For n > 2, a(n) is of the form 2^e_1 * p_2^e_2 * ... * p_m^e_m, where e_m = 1 and e = floor(log_2(p_m)) <= e_1. Therefore, 2^e * p_m^e_m is a primitive Zumkeler number (A180332). Therefore, 2^e_1 * p_m^e_m is a Zumkeller number (A083207). Therefore, for n > 2, a(n) = 2^e_1 * p_m^e_m * r, where r is relatively prime to 2*p_m, is a Zumkeller number (see my proof at A002182 for details). - Ivan N. Ianakiev, May 10 2020
For n > 1, 2|(a(n)+2) ... n|(a(n)+n), so a(n)+2 .. a(n)+n are all composite and (part of) a prime gap of at least n. (Compare n!+2 .. n!+n). - Stephen E. Witham, Oct 09 2021

Examples

			LCM of {1,2,3,4,5,6} = 60. The primes up to 6 are 2, 3 and 5. floor(log(6)/log(2)) = 2 so the exponent of 2 is 2.
floor(log(6)/log(3)) = 1 so the exponent of 3 is 1.
floor(log(6)/log(5)) = 1 so the exponent of 5 is 1. Therefore, a(6) = 2^2 * 3^1 * 5^1 = 60. - _David A. Corneth_, Jun 02 2017
		

References

  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 365.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Row products of A133233.
Cf. A025528 (number of prime factors of a(n) with multiplicity).
Cf. A275120 (lengths of runs of consecutive equal terms), A276781 (ordinal transform from term a(1)=1 onward).

Programs

  • Haskell
    a003418 = foldl lcm 1 . enumFromTo 2
    -- Reinhard Zumkeller, Apr 04 2012, Apr 25 2011
    
  • Magma
    [1] cat [Exponent(SymmetricGroup(n)) : n in [1..28]]; // Arkadiusz Wesolowski, Sep 10 2013
    
  • Magma
    [Lcm([1..n]): n in [0..30]]; // Bruno Berselli, Feb 06 2015
    
  • Maple
    A003418 := n-> lcm(seq(i,i=1..n));
    HalfFarey := proc(n) local a,b,c,d,k,s; a := 0; b := 1; c := 1; d := n; s := NULL; do k := iquo(n + b, d); a, b, c, d := c, d, k*c - a, k*d - b; if 2*a > b then break fi; s := s,(a/b); od: [s] end: LCM := proc(n) local i; (1/2)*mul(2*sin(Pi*i),i=HalfFarey(n))^2 end: # Peter Luschny
    # next Maple program:
    a:= proc(n) option remember; `if`(n=0, 1, ilcm(n, a(n-1))) end:
    seq(a(n), n=0..33);  # Alois P. Heinz, Jun 10 2021
  • Mathematica
    Table[LCM @@ Range[n], {n, 1, 40}] (* Stefan Steinerberger, Apr 01 2006 *)
    FoldList[ LCM, 1, Range@ 28]
    A003418[0] := 1; A003418[1] := 1; A003418[n_] := A003418[n] = LCM[n,A003418[n-1]]; (* Enrique Pérez Herrero, Jan 08 2011 *)
    Table[Product[Prime[i]^Floor[Log[Prime[i], n]], {i, PrimePi[n]}], {n, 0, 28}] (* Wei Zhou, Jun 25 2011 *)
    Table[Product[Cyclotomic[n, 1], {n, 2, m}], {m, 0, 28}] (* Fred Daniel Kline, May 22 2014 *)
    a1[n_] := 1/12 (Pi^2+3(-1)^n (PolyGamma[1,1+n/2] - PolyGamma[1,(1+n)/2])) // Simplify
    a[n_] := Denominator[Sqrt[a1[n]]];
    Table[If[IntegerQ[a[n]], a[n], a[n]*(a[n])[[2]]], {n, 0, 28}] (* Gerry Martens, Apr 07 2018 [Corrected by Vaclav Kotesovec, Jul 16 2021] *)
  • PARI
    a(n)=local(t); t=n>=0; forprime(p=2,n,t*=p^(log(n)\log(p))); t
    
  • PARI
    a(n)=if(n<1,n==0,1/content(vector(n,k,1/k)))
    
  • PARI
    a(n)=my(v=primes(primepi(n)),k=sqrtint(n),L=log(n+.5));prod(i=1,#v,if(v[i]>k,v[i],v[i]^(L\log(v[i])))) \\ Charles R Greathouse IV, Dec 21 2011
    
  • PARI
    a(n)=lcm(vector(n,i,i)) \\ Bill Allombert, Apr 18 2012 [via Charles R Greathouse IV]
    
  • PARI
    n=1; lim=100; i=1; j=1; until(n==lim, a=lcm(j,i+1); i++; j=a; n++; print(n" "a);); \\ Mike Winkler, Sep 07 2013
    
  • Python
    from functools import reduce
    from operator import mul
    from sympy import sieve
    def integerlog(n,b): # find largest integer k>=0 such that b^k <= n
        kmin, kmax = 0,1
        while b**kmax <= n:
            kmax *= 2
        while True:
            kmid = (kmax+kmin)//2
            if b**kmid > n:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmin
    def A003418(n):
        return reduce(mul,(p**integerlog(n,p) for p in sieve.primerange(1,n+1)),1) # Chai Wah Wu, Mar 13 2021
    
  • Python
    # generates initial segment of sequence
    from math import gcd
    from itertools import accumulate
    def lcm(a, b): return a * b // gcd(a, b)
    def aupton(nn): return [1] + list(accumulate(range(1, nn+1), lcm))
    print(aupton(30)) # Michael S. Branicky, Jun 10 2021
  • Sage
    [lcm(range(1,n)) for n in range(1, 30)] # Zerinvary Lajos, Jun 06 2009
    
  • Scheme
    (define (A003418 n) (let loop ((n n) (m 1)) (if (zero? n) m (loop (- n 1) (lcm m n))))) ;; Antti Karttunen, Jan 03 2018
    

Formula

The prime number theorem implies that lcm(1,2,...,n) = exp(n(1+o(1))) as n -> infinity. In other words, log(lcm(1,2,...,n))/n -> 1 as n -> infinity. - Jonathan Sondow, Jan 17 2005
a(n) = Product (p^(floor(log n/log p))), where p runs through primes not exceeding n (i.e., primes 2 through A007917(n)). - Lekraj Beedassy, Jul 27 2004
Greg Martin showed that a(n) = lcm(1,2,3,...,n) = Product_{i = Farey(n), 0 < i < 1} 2*Pi/Gamma(i)^2. This can be rewritten (for n > 1) as a(n) = (1/2)*(Product_{i = Farey(n), 0 < i <= 1/2} 2*sin(i*Pi))^2. - Peter Luschny, Aug 08 2009
Recursive formula useful for computations: a(0)=1; a(1)=1; a(n)=lcm(n,a(n-1)). - Enrique Pérez Herrero, Jan 08 2011
From Enrique Pérez Herrero, Jun 01 2011: (Start)
a(n)/a(n-1) = A014963(n).
if n is a prime power p^k then a(n)=a(p^k)=p*a(n-1), otherwise a(n)=a(n-1).
a(n) = Product_{k=2..n} (1 + (A007947(k)-1)*floor(1/A001221(k))), for n > 1. (End)
a(n) = A079542(n+1, 2) for n > 1.
a(n) = exp(Sum_{k=1..n} Sum_{d|k} moebius(d)*log(k/d)). - Peter Luschny, Sep 01 2012
a(n) = A025529(n) - A027457(n). - Eric Desbiaux, Mar 14 2013
a(n) = exp(Psi(n)) = 2 * Product_{k=2..A002088(n)} (1 - exp(2*Pi*i * A038566(k+1) / A038567(k))), where i is the imaginary unit, and Psi the second Chebyshev's function. - Eric Desbiaux, Aug 13 2014
a(n) = A064446(n)*A038610(n). - Anthony Browne, Jun 16 2016
a(n) = A000142(n) / A025527(n) = A000793(n) * A225558(n). - Antti Karttunen, Jun 02 2017
log(a(n)) = Sum_{k>=1} (A309229(n, k)/k - 1/k). - Mats Granvik, Aug 10 2019
From Petros Hadjicostas, Jul 24 2020: (Start)
Nair (1982) proved that 2^n <= a(n) <= 4^n for n >= 9. See also Farhi (2009). Nair also proved that
a(n) = lcm(m*binomial(n,m): 1 <= m <= n) and
a(n) = gcd(a(m)*binomial(n,m): n/2 <= m <= n). (End)
Sum_{n>=1} 1/a(n) = A064859. - Bernard Schott, Aug 24 2020

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

Views

Author

Keywords

Comments

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

References

  • 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).

Crossrefs

Programs

  • Haskell
    import Data.List (nub)
    a002183 n = a002183_list !! (n-1)
    a002183_list = nub $ map (a000005 . a061799) [1..]
    -- Reinhard Zumkeller, Apr 01 2011
  • Mathematica
    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 *)

Formula

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

Extensions

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

A004490 Colossally abundant numbers: m for which there is a positive exponent epsilon such that sigma(m)/m^{1 + epsilon} >= sigma(k)/k^{1 + epsilon} for all k > 1, so that m attains the maximum value of sigma(m)/m^{1 + epsilon}.

Original entry on oeis.org

2, 6, 12, 60, 120, 360, 2520, 5040, 55440, 720720, 1441440, 4324320, 21621600, 367567200, 6983776800, 160626866400, 321253732800, 9316358251200, 288807105787200, 2021649740510400, 6064949221531200, 224403121196654400
Offset: 1

Views

Author

N. J. A. Sloane, Jan 22 2001

Keywords

References

  • S. Ramanujan, Highly composite numbers, Proc. London Math. Soc., 14 (1915), 347-407. Reprinted in Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, pp. 78-129. See esp. pp. 87, 115.

Crossrefs

A subsequence of A004394 (superabundant numbers).
Cf. A002093 (highly abundant numbers), A002182, A005101 (abundant numbers), A006038, A189228 (superabundant numbers that are not colossally abundant).

Formula

a(n) = Product_{k=1..n} A073751(k). - Jeppe Stig Nielsen, Nov 28 2021
Showing 1-10 of 47 results. Next