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 44 results. Next

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

N. J. A. Sloane

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.

Links

Crossrefs

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

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

A050376 "Fermi-Dirac primes": numbers of the form p^(2^k) where p is prime and k >= 0.

Original entry on oeis.org

2, 3, 4, 5, 7, 9, 11, 13, 16, 17, 19, 23, 25, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241
Offset: 1

Views

Author

Christian G. Bower, Nov 15 1999

Keywords

Comments

Every number n is a product of a unique subset of these numbers. This is sometimes called the Fermi-Dirac factorization of n (see A182979). Proof: In the prime factorization n = Product_{j>=1} p(j)^e(j) expand every exponent e(j) as binary number and pick the terms of this sequence corresponding to the positions of the ones in binary (it is clear that both n and n^2 have the same number of factors in this sequence, and that each factor appears with exponent 1 or 0).
Or, a(1) = 2; for n>1, a(n) = smallest number which cannot be obtained as the product of previous terms. This is evident from the unique factorization theorem and the fact that every number can be expressed as the sum of powers of 2. - Amarnath Murthy, Jan 09 2002
Except for the first term, same as A084400. - David Wasserman, Dec 22 2004
The least number having 2^n divisors (=A037992(n)) is the product of the first n terms of this sequence according to Ramanujan.
According to the Bose-Einstein distribution of particles, an unlimited number of particles may occupy the same state. On the other hand, according to the Fermi-Dirac distribution, no two particles can occupy the same state (by the Pauli exclusion principle). Unique factorizations of the positive integers by primes (A000040) and over terms of A050376 one can compare with two these distributions in physics of particles. In the correspondence with this, the factorizations over primes one can call "Bose-Einstein factorizations", while the factorizations over distinct terms of A050376 one can call "Fermi-Dirac factorizations". - Vladimir Shevelev, Apr 16 2010
The numbers of the form p^(2^k), where p is prime and k >= 0, might thus be called the "Fermi-Dirac primes", while the classic primes might be called the "Bose-Einstein primes". - Daniel Forgues, Feb 11 2011
In the theory of infinitary divisors, the most natural name of the terms is "infinitary primes" or "i-primes". Indeed, n is in the sequence, if and only if it has only two infinitary divisors. Since 1 and n are always infinitary divisors of n>1, an i-prime has no other infinitary divisors. - Vladimir Shevelev, Feb 28 2011
{a(n)} is the minimal set including all primes and closed with respect to squaring. In connection with this, note that n and n^2 have the same number of factors in their Fermi-Dirac representations. - Vladimir Shevelev, Mar 16 2012
In connection with this sequence, call an integer compact if the factors in its Fermi-Dirac factorization are pairwise coprime. The density of such integers equals (6/Pi^2)*Product_{prime p} (1+(Sum_{i>=1} p^(-(2^i-1))/(p+1))) = 0.872497... It is interesting that there exist only 7 compact factorials listed in A169661. - Vladimir Shevelev, Mar 17 2012
The first k terms of the sequence solve the following optimization problem:
Let x_1, x_2,..., x_k be integers with the restrictions: 2<=x_1A064547(Product{i=1..k} x_i) >= k. Let the goal function be Product_{i=1..k} x_i. Then the minimal value of the goal function is Product_{i=1..k} a(i). - Vladimir Shevelev, Apr 01 2012
From Joerg Arndt, Mar 11 2013: (Start)
Similarly to the first comment, for the sequence "Numbers of the form p^(3^k) or p^(2*3^k) where p is prime and k >= 0" one obtains a factorization into distinct factors by using the ternary expansion of the exponents (here n and n^3 have the same number of such factors).
The generalization to base r would use "Numbers of the form p^(r^k), p^(2*r^k), p^(3*r^k), ..., p^((r-1)*r^k) where p is prime and k >= 0" (here n and n^r have the same number of (distinct) factors). (End)
The first appearance of this sequence as a multiplicative basis in number theory with some new notions, formulas and theorems may have been in my 1981 paper (see the Abramovich reference). - Vladimir Shevelev, Apr 27 2014
Numbers n for which A064547(n) = 1. - Antti Karttunen, Feb 10 2016
Lexicographically earliest sequence of distinct nonnegative integers such that no term is a product of 2 or more distinct terms. Removing the distinctness requirement, the sequence becomes A000040 (the prime numbers); and the equivalent sequence where the product is of 2 distinct terms is A026416 (without its initial term, 1). - Peter Munn, Mar 05 2019
The sequence was independently developed as a multiplicative number system in 1985-1986 (and first published in 1995, see the Uhlmann reference) using a proof method involving representations of positive integers as sums of powers of 2. This approach offers an arguably simpler and more flexible means for analyzing the sequence. - Jeffrey K. Uhlmann, Nov 09 2022

Examples

			Prime powers which are not terms of this sequence:
  8 = 2^3 = 2^(1+2), 27 = 3^3 = 3^(1+2), 32 = 2^5 = 2^(1+4),
  64 = 2^6 = 2^(2+4), 125 = 5^3 = 5^(1+2), 128 = 2^7 = 2^(1+2+4)
"Fermi-Dirac factorizations":
  6 = 2*3, 8 = 2*4, 24 = 2*3*4, 27 = 3*9, 32 = 2*16, 64 = 4*16,
  108 = 3*4*9, 120 = 2*3*4*5, 121 = 121, 125 = 5*25, 128 = 2*4*16.

References

  • V. S. Abramovich, On an analog of the Euler function, Proceeding of the North-Caucasus Center of the Academy of Sciences of the USSR (Rostov na Donu) (1981) No. 2, 13-17 (Russian; MR0632989(83a:10003)).
  • S. Ramanujan, Highly Composite Numbers, Collected Papers of Srinivasa Ramanujan, p. 125, Ed. G. H. Hardy et al., AMS Chelsea 2000.
  • V. S. Shevelev, Multiplicative functions in the Fermi-Dirac arithmetic, Izvestia Vuzov of the North-Caucasus region, Nature sciences 4 (1996), 28-43 (in Russian; MR 2000f: 11097, pp. 3912-3913).
  • J. K. Uhlmann, Dynamic map building and localization: new theoretical foundations, Doctoral Dissertation, University of Oxford, Appendix 16, 1995.

Links

Crossrefs

Cf. A000040 (primes, is a subsequence), A026416, A026477, A037992 (partial products), A050377-A050380, A052330, A064547, A066724, A084400, A176699, A182979.
Cf. A268388 (complement without 1).
Cf. A124010, subsequence of A000028, A000961, A213925, A223490.
Cf. A228520, A186945 (Fermi-Dirac analog of Ramanujan primes, A104272, and Labos primes, A080359).
Cf. also A268385, A268391, A268392.

Programs

  • Haskell
    a050376 n = a050376_list !! (n-1)
a050376_list = filter ((== 1) . a209229 . a100995) [1..]
-- Reinhard Zumkeller, Mar 19 2013
  • Maple
    isA050376 := proc(n)
    local f,e;
    f := ifactors(n)[2] ;
    if nops(f) = 1 then
        e := op(2,op(1,f)) ;
        if isA000079(e) then
            true;
        else
            false;
        end if;
    else
        false;
    end if;
end proc:
A050376 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        2 ;
    else
        for a from procname(n-1)+1 do
            if isA050376(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, May 26 2017
  • Mathematica
    nn = 300; t = {}; k = 1; While[lim = nn^(1/k); lim > 2,  t = Join[t, Prime[Range[PrimePi[lim]]]^k]; k = 2 k]; t = Union[t] (* T. D. Noe, Apr 05 2012 *)
  • PARI
    {a(n)= local(m, c, k, p); if(n<=1, 2*(n==1), n--; c=0; m=2; while( cMichael Somos, Apr 15 2005; edited by Michel Marcus, Aug 07 2021
  • PARI
    lst(lim)=my(v=primes(primepi(lim)),t); forprime(p=2,sqrt(lim),t=p; while((t=t^2)<=lim,v=concat(v,t))); vecsort(v) \\ Charles R Greathouse IV, Apr 10 2012
  • PARI
    is_A050376(n)=2^#binary(n=isprimepower(n))==n*2 \\ M. F. Hasler, Apr 08 2015
  • PARI
    ispow2(n)=n && n>>valuation(n,2)==1
is(n)=ispow2(isprimepower(n)) \\ Charles R Greathouse IV, Sep 18 2015
  • PARI
    isok(n)={my(e=isprimepower(n)); e && !bitand(e,e-1)} \\ Andrew Howroyd, Oct 16 2024
  • Python
    from sympy import isprime, perfect_power
def ok(n):
  if isprime(n): return True
  answer = perfect_power(n)
  if not answer: return False
  b, e = answer
  if not isprime(b): return False
  while e%2 == 0: e //= 2
  return e == 1
def aupto(limit):
  alst, m = [], 1
  for m in range(1, limit+1):
    if ok(m): alst.append(m)
  return alst
print(aupto(241)) # Michael S. Branicky, Feb 03 2021
  • Python
    from sympy import primepi, integer_nthroot
def A050376(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(integer_nthroot(x,1<Chai Wah Wu, Feb 18-19 2025
  • Scheme
    (define A050376 (MATCHING-POS 1 1 (lambda (n) (= 1 (A064547 n)))))
;; Requires also my IntSeq-library. - Antti Karttunen, Feb 09 2016

Formula

From Vladimir Shevelev, Mar 16 2012: (Start)
Product_{i>=1} a(i)^k_i = n!, where k_i = floor(n/a(i)) - floor(n/a(i)^2) + floor(n/a(i)^3) - floor(n/a(i)^4) + ...
Denote by A(x) the number of terms not exceeding x.
Then A(x) = pi(x) + pi(x^(1/2)) + pi(x^(1/4)) + pi(x^(1/8)) + ...
Conversely, pi(x) = A(x) - A(sqrt(x)). For example, pi(37) = A(37) - A(6) = 16-4 = 12. (End)
A209229(A100995(a(n))) = 1. - Reinhard Zumkeller, Mar 19 2013
From Vladimir Shevelev, Aug 31 2013: (Start)
A Fermi-Dirac analog of Euler product: Zeta(s) = Product_{k>=1} (1+a(k)^(-s)), for s > 1.
In particular, Product_{k>=1} (1+a(k)^(-2)) = Pi^2/6. (End)
a(n) = A268385(A268392(n)). [By their definitions.] - Antti Karttunen, Feb 10 2016
A000040 union A001248 union A030514 union A179645 union A030635 union .... - R. J. Mathar, May 26 2017

Extensions

Edited by Charles R Greathouse IV, Mar 17 2010
More examples from Daniel Forgues, Feb 09 2011

A005179 Smallest number with exactly n divisors.

Original entry on oeis.org

1, 2, 4, 6, 16, 12, 64, 24, 36, 48, 1024, 60, 4096, 192, 144, 120, 65536, 180, 262144, 240, 576, 3072, 4194304, 360, 1296, 12288, 900, 960, 268435456, 720, 1073741824, 840, 9216, 196608, 5184, 1260, 68719476736, 786432, 36864, 1680, 1099511627776, 2880
Offset: 1

Views

Author

N. J. A. Sloane, David Singmaster

Keywords

Comments

A number n is called ordinary iff a(n)=A037019(n). Brown shows that the ordinary numbers have density 1 and all squarefree numbers are ordinary. See A072066 for the extraordinary or exceptional numbers. - M. F. Hasler, Oct 14 2014
All terms are in A025487. Therefore, a(n) is even for n > 1. - David A. Corneth, Jun 23 2017 [corrected by Charles R Greathouse IV, Jul 05 2023]

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 1, p. 52.
  • Joe Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 86.
  • 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 89.

Links

Crossrefs

Cf. A000005, A007416, A099316, A003586, A025487, A099311, A099313, A050376, A037992, A061799, A162247, A262981, A262983.

Programs

  • Haskell
    import Data.List (elemIndex)
import Data.Maybe (fromJust)
a005179 n = succ $ fromJust $ elemIndex n $ map a000005 [1..]
-- Reinhard Zumkeller, Apr 01 2011
  • Maple
    A005179_list := proc(SearchLimit, ListLength)
local L, m, i, d; m := 1;
L := array(1..ListLength,[seq(0,i=1..ListLength)]);
for i from 1 to SearchLimit while m <= ListLength do
  d := numtheory[tau](i);
  if d <= ListLength and 0 = L[d] then L[d] := i;
  m := m + 1; fi
od:
print(L) end: A005179_list(65537,18);
# If a '0' appears in the list the search limit has to be increased. - Peter Luschny, Mar 09 2011
# alternative
# Construct list of ordered lists of factorizations of n with
# minimum divisors mind.
# Returns a list with A001055(n) entries if called with mind=2.
# Example: print(ofact(10^3,2))
ofact := proc(n,mind)
    local fcts,d,rec,r ;
    fcts := [] ;
    for d in numtheory[divisors](n) do
        if d >= mind then
            if d = n then
                fcts := [op(fcts),[n]] ;
            else
                # recursive call supposed one more factor fixed now
                rec := procname(n/d,max(d,mind)) ;
                for r in rec do
                    fcts := [op(fcts),[d,op(r)]] ;
                end do:
            end if;
        end if;
    end do:
    return fcts ;
end proc:
A005179 := proc(n)
    local Lexp,a,eList,cand,maxxrt ;
    if n = 1 then
        return 1;
    end if;
    Lexp := ofact(n,2) ;
    a := 0 ;
    for eList in Lexp do
        maxxrt := ListTools[Reverse](eList) ;
        cand := mul( ithprime(i)^ ( op(i,maxxrt)-1),i=1..nops(maxxrt)) ;
        if a =0 or cand < a then
            a := cand ;
        end if;
    end do:
    a ;
end proc:
seq(A005179(n),n=1..40) ; # R. J. Mathar, Jun 06 2024
  • Mathematica
    a = Table[ 0, {43} ]; Do[ d = Length[ Divisors[ n ]]; If[ d < 44 && a[[ d ]] == 0, a[[ d]] = n], {n, 1, 1099511627776} ]; a
(* Second program: *)
Function[s, Map[Lookup[s, #] &, Range[First@ Complement[Range@ Max@ #, #] - 1]] &@ Keys@ s]@ Map[First, KeySort@ PositionIndex@ Table[DivisorSigma[0, n], {n, 10^7}]] (* Michael De Vlieger, Dec 11 2016, Version 10 *)
mp[1, m_] := {{}}; mp[n_, 1] := {{}}; mp[n_?PrimeQ, m_] := If[m < n, {}, {{n}}]; mp[n_, m_] := Join @@ Table[Map[Prepend[#, d] &, mp[n/d, d]], {d, Select[Rest[Divisors[n]], # <= m &]}]; mp[n_] := mp[n, n]; Table[mulpar = mp[n] - 1; Min[Table[Product[Prime[s]^mulpar[[j, s]], {s, 1, Length[mulpar[[j]]]}], {j, 1, Length[mulpar]}]], {n, 1, 100}] (* Vaclav Kotesovec, Apr 04 2021 *)
a[n_] := Module[{e = f[n] - 1}, Min[Times @@@ ((Prime[Range[Length[#], 1, -1]]^#) & /@ e)]]; Array[a, 100] (* Amiram Eldar, Jul 26 2025 using the function f by T. D. Noe at A162247 *)
  • PARI
    (prodR(n,maxf)=my(dfs=divisors(n),a=[],r); for(i=2,#dfs, if( dfs[i]<=maxf, if(dfs[i]==n, a=concat(a,[[n]]), r=prodR(n/dfs[i],min(dfs[i],maxf)); for(j=1,#r, a=concat(a,[concat(dfs[i],r[j])]))))); a); A005179(n)=my(pf=prodR(n,n),a=1,b); for(i=1,#pf, b=prod(j=1,length(pf[i]),prime(j)^(pf[i][j]-1)); if(bA005179(n)", ")) \\ R. J. Mathar, May 26 2008, edited by M. F. Hasler, Oct 11 2014
  • Python
    from math import prod
from sympy import isprime, divisors, prime
def A005179(n):
    def mult_factors(n):
        if isprime(n):
            return [(n,)]
        c = []
        for d in divisors(n,generator=True):
            if 1Chai Wah Wu, Aug 17 2024

Formula

a(p) = 2^(p-1) for primes p: a(A000040(n)) = A061286(n); a(p^2) = 6^(p-1) for primes p: a(A001248(n)) = A061234(n); a(p*q) = 2^(q-1)*3^(p-1) for primes p<=q: a(A001358(n)) = A096932(n); a(p*m*q) = 2^(q-1) * 3^(m-1) * 5^(p-1) for primes pA005179(A007304(n)) = A061299(n). - Reinhard Zumkeller, Jul 15 2004 [This can be continued to arbitrarily many distinct prime factors since no numbers in A072066 (called "exceptional" or "extraordinary") are squarefree. - Jianing Song, Jul 18 2025]
a(p^n) = (2*3...*p_n)^(p-1) for p > log p_n / log 2. Unpublished proof from Andrzej Schinzel. - Thomas Ordowski, Jul 22 2005
If p is a prime and n=p^k then a(p^k)=(2*3*...*s_k)^(p-1) where (s_k) is the numbers of the form q^(p^j) for every q and j>=0, according to Grost (1968), Theorem 4. For example, if p=2 then a(2^k) is the product of the first k members of the A050376 sequence: number of the form q^(2^j) for j>=0, according to Ramanujan (1915). - Thomas Ordowski, Aug 30 2005
a(2^k) = A037992(k). - Thomas Ordowski, Aug 30 2005
a(n) <= A037019(n) with equality except for n in A072066. - M. F. Hasler, Jun 15 2022

Extensions

More terms from David W. Wilson

A037445 Number of infinitary divisors (or i-divisors) of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 4, 2, 2, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 8, 4, 8, 4, 4, 2, 8, 2, 4, 4, 4, 4, 8, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 4, 8, 2, 4, 2, 4, 2, 8, 4, 4, 4, 8, 2, 8, 4, 4, 4, 4, 4, 8, 2, 4, 4, 4, 2, 8, 2, 8, 8
Offset: 1

Views

Author

Yasutoshi Kohmoto

Keywords

Comments

A divisor of n is called infinitary if it is a product of divisors of the form p^{y_a 2^a}, where p^y is a prime power dividing n and sum_a y_a 2^a is the binary representation of y.
The smallest number m with exactly 2^n infinitary divisors is A037992(n); for these values m, a(m) increases also to a new record. - Bernard Schott, Mar 09 2023

Examples

			For n = 8, n = 2^3 = 2^"11" (writing 3 in binary) so the infinitary divisors are 2^"00" = 1, 2^"01" = 2, 2^"10" = 4 and 2^"11" = 8, so a(8) = 4.
For n = 90, n = 2*5*9 where 2,5,9 are in A050376, so a(90) = 2^3 = 8.

Links

Crossrefs

Cf. A000120, A001316, A004607, A007358, A007357, A037992, A038148, A049417, A064547, A074848, A077609, A124010, A156552, A286575.

Programs

  • Haskell
    a037445 = product . map (a000079 . a000120) . a124010_row
-- Reinhard Zumkeller, Mar 19 2013
  • Maple
    A037445 := proc(n)
    local a,p;
    a := 1 ;
    for p in ifactors(n)[2] do
        a := a*2^wt(p[2]) ;
    end do:
    a ;
end proc: # R. J. Mathar, May 16 2016
  • Mathematica
    Table[Length@((Times @@ (First[it]^(#1 /. z -> List)) & ) /@
Flatten[Outer[z, Sequence @@ bitty /@
Last[it = Transpose[FactorInteger[k]]], 1]]), {k, 2, 240}]
bitty[k_] := Union[Flatten[Outer[Plus, Sequence @@ ({0, #1} & ) /@ Union[2^Range[0, Floor[Log[2, k]]]*Reverse[IntegerDigits[k, 2]]]]]]
y[n_] := Select[Range[0, n], BitOr[n, # ] == n & ] divisors[Infinity][1] := {1}
divisors[Infinity][n_] := Sort[Flatten[Outer[Times, Sequence @@ (FactorInteger[n] /. {p_, m_Integer} :> p^y[m])]]] Length /@ divisors[Infinity] /@ Range[105] (* Paul Abbott (paul(AT)physics.uwa.edu.au), Apr 29 2005 *)
a[1] = 1; a[n_] := Times @@ Flatten[ 2^DigitCount[#, 2, 1]&  /@ FactorInteger[n][[All, 2]] ]; Table[a[n], {n, 1, 105}] (* Jean-François Alcover, Aug 19 2013, after Reinhard Zumkeller *)
  • PARI
    A037445(n) = factorback(apply(a -> 2^hammingweight(a), factorint(n)[,2])) \\ Andrew Lelechenko, May 10 2014
  • Python
    from sympy import factorint
def wt(n): return bin(n).count("1")
def a(n):
    f=factorint(n)
    return 2**sum([wt(f[i]) for i in f]) # Indranil Ghosh, May 30 2017
  • Scheme
    (define (A037445 n) (if (= 1 n) n (* (A001316 (A067029 n)) (A037445 (A028234 n))))) ;; Antti Karttunen, May 28 2017

Formula

Multiplicative with a(p^e) = 2^A000120(e). - David W. Wilson, Sep 01 2001
Let n = q_1*...*q_k, where q_1,...,q_k are different terms of A050376. Then a(n) = 2^k (the number of subsets of a set with k elements is 2^k). - Vladimir Shevelev, Feb 19 2011.
a(n) = Product_{k=1..A001221(n)} A000079(A000120(A124010(n,k))). - Reinhard Zumkeller, Mar 19 2013
From Antti Karttunen, May 28 2017: (Start)
a(n) = A286575(A156552(n)). [Because multiplicative with a(p^e) = A001316(e).]
a(n) = 2^A064547(n). (End)
a(A037992(n)) = 2^n. - Bernard Schott, Mar 10 2023

Extensions

Corrected and extended by Naohiro Nomoto, Jun 21 2001

A064547 Sum of binary digits (or count of 1-bits) in the exponents of the prime factorization of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 2, 2, 1, 2, 2, 3, 1, 3, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 3, 2, 3, 2, 2, 1, 3, 1, 2, 2, 2, 2, 3, 1, 2, 2, 3, 1, 3, 1, 2, 2, 2, 2, 3, 1, 2, 1, 2, 1, 3, 2, 2, 2, 3, 1, 3, 2, 2, 2, 2, 2, 3, 1, 2, 2, 2, 1, 3, 1, 3, 3
Offset: 1

Views

Author

Wouter Meeussen, Oct 09 2001

Keywords

Comments

This sequence is different from A058061 for n containing 6th, 8th, ..., k-th powers in its prime decomposition, where k runs through the integers missing from A064548.
For n > 1, n is a product of a(n) distinct members of A050376. - Matthew Vandermast, Jul 13 2004
For n > 1: a(n) = length of n-th row in A213925. - Reinhard Zumkeller, Mar 20 2013
Number of Fermi-Dirac factors of n. - Peter Munn, Dec 27 2019

Examples

			For n = 54, n = 2^1 * 3^3 with exponents (1) and (11) in binary, so a(54) = A000120(1) + A000120(3) = 1 + 2 = 3.

Links

Crossrefs

Cf. A000028 (positions of odd terms), A000379 (of even terms).
Cf. A050376 (positions of ones), A268388 (terms larger than ones).
Row lengths of A213925.
A000120, A007814, A028234, A037445, A052331, A064989, A067029, A156552, A223491, A286574 are used in formulas defining this sequence.
Cf. A005117, A058061 (to which A064548 relates), A138302.
Cf. other sequences counting factors of n: A001221, A001222.
Cf. other sequences where a(n) depends only on the prime signature of n: A181819, A267116, A268387.
A003961, A007913, A008833, A059895, A059896, A059897, A225546 are used to express relationship between terms of this sequence.
Cf. A077761, A176699, A037992, A382294.

Programs

  • Haskell
    a064547 1 = 0
a064547 n = length $ a213925_row n  -- Reinhard Zumkeller, Mar 20 2013
  • Maple
    expts:=proc(n) local t1,t2,t3,t4,i; if n=1 then RETURN([0]); fi; if isprime(n) then RETURN([1]); fi; t1:=ifactor(n); if nops(factorset(n))=1 then RETURN([op(2,t1)]); fi; t2:=nops(t1); t3:=[]; for i from 1 to t2 do t4:=op(i,t1); if nops(t4) = 1 then t3:=[op(t3),1]; else t3:=[op(t3),op(2,t4)]; fi; od; RETURN(t3); end;
A000120 := proc(n) local w,m,i; w := 0; m := n; while m > 0 do i := m mod 2; w := w+i; m := (m-i)/2; od; w; end:
LamMos:= proc(n) local t1,t2,t3,i; t1:=expts(n); add( A000120(t1[i]),i=1..nops(t1)); end; # N. J. A. Sloane, Dec 20 2007
# alternative Maple program:
A064547:= proc(n) local F;
F:= ifactors(n)[2];
add(convert(convert(f[2],base,2),`+`),f=F)
end proc:
map(A064547,[$1..100]); # Robert Israel, May 17 2016
  • Mathematica
    Table[Plus@@(DigitCount[Last/@FactorInteger[k], 2, 1]), {k, 105}]
  • PARI
    a(n) = {my(f = factor(n)[,2]); sum(k=1, #f, hammingweight(f[k]));} \\ Michel Marcus, Feb 10 2016
  • Python
    from sympy import factorint
def wt(n): return bin(n).count("1")
def a(n):
    f=factorint(n)
    return sum([wt(f[i]) for i in f]) # Indranil Ghosh, May 30 2017
  • Scheme
    ;; uses memoizing-macro definec
(definec (A064547 n) (cond ((= 1 n) 0) (else (+ (A000120 (A067029 n)) (A064547 (A028234 n))))))
;; Antti Karttunen, Feb 09 2016
  • Scheme
    ;; uses memoizing-macro definec
(definec (A064547 n) (if (= 1 n) 0 (+ (A000120 (A007814 n)) (A064547 (A064989 n)))))
;; Antti Karttunen, Feb 09 2016

Formula

a(m*n) <= a(m)*a(n). - Reinhard Zumkeller, Mar 20 2013
From Antti Karttunen, Feb 09 2016: (Start)
a(1) = 0, and for n > 1, a(n) = A000120(A067029(n)) + a(A028234(n)).
a(1) = 0, and for n > 1, a(n) = A000120(A007814(n)) + a(A064989(n)).
(End)
a(n) = log_2(A037445(n)). - Vladimir Shevelev, May 13 2016
a(n) = A286574(A156552(n)). - Antti Karttunen, May 28 2017
Additive with a(p^e) = A000120(e). - Jianing Song, Jul 28 2018
a(n) = A000120(A052331(n)). - Peter Munn, Aug 26 2019
From Peter Munn, Dec 18 2019: (Start)
a(A000379(n)) mod 2 = 0.
a(A000028(n)) mod 2 = 1.
A001221(n) <= a(n) <= A001222(n).
A001221(n) < a(n) => a(n) < A001222(n).
a(n) = A001222(n) if and only if n is in A005117.
a(n) = A001221(n) if and only if n is in A138302.
a(n^2) = a(n).
a(A003961(n)) = a(n).
a(A225546(n)) = a(n).
a(n) = a(A007913(n)) + a(A008833(n)).
a(A050376(n)) = 1.
a(A059897(n,k)) + 2 * a(A059895(n,k)) = a(n) + a(k).
a(A059896(n,k)) + a(A059895(n,k)) = a(n) + a(k).
Alternative definition: a(1) = 0; a(n * m) = a(n) + 1 for m = A050376(k) > A223491(n).
(End)
Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761) and C = Sum_{p prime} f(1/p) = 0.13605447049622836522... (A382294), where f(x) = -x + Sum_{k>=0} x^(2^k)/(1+x^(2^k)). - Amiram Eldar, Sep 28 2023
a(n) << log n/log log n. - Charles R Greathouse IV, Nov 29 2024

A061283 Smallest number with exactly 2n-1 divisors.

Original entry on oeis.org

1, 4, 16, 64, 36, 1024, 4096, 144, 65536, 262144, 576, 4194304, 1296, 900, 268435456, 1073741824, 9216, 5184, 68719476736, 36864, 1099511627776, 4398046511104, 3600, 70368744177664, 46656, 589824, 4503599627370496, 82944
Offset: 1

Views

Author

Labos Elemer, May 22 2001

Keywords

Comments

The terms are always squares (because the divisors of a nonsquare N come in pairs, d and N/d, and so their number is always even - N. J. A. Sloane, Dec 26 2018).

Examples

			For n=15, a(15)=144 with 15 divisors: 1,2,3,4,6,8,9,12,16,18,24,36,48,72 and 144.

Links

Crossrefs

Cf. A000005, A000290, A005408, A003680, A016017, A037992, A055079, A048691.
Second bisection of A005179.

Programs

  • Mathematica
    mp[1, m_] := {{}}; mp[n_, 1] := {{}}; mp[n_?PrimeQ, m_] := If[m < n, {}, {{n}}]; mp[n_, m_] := Join @@ Table[Map[Prepend[#, d] &, mp[n/d, d]], {d, Select[Rest[Divisors[n]], # <= m &]}]; mp[n_] := mp[n, n]; Table[mulpar = mp[2*n-1] - 1; Min[Table[Product[Prime[s]^mulpar[[j, s]], {s, 1, Length[mulpar[[j]]]}], {j, 1, Length[mulpar]}]], {n, 1, 100}] (* Vaclav Kotesovec, Apr 04 2021 *)

Formula

a(n) = Min{k | A000005(k)=2n-1}.
a((p+1)/2) = 2^(p-1) for odd prime p. [Corrected by Jianing Song, Aug 30 2021]
From Jianing Song, Aug 30 2021: (Start)
a(n) = A016017(n)^2.
a(n) <= 2^(2n-2), where the equality holds if and only if n=1 or 2n-1 is prime. (End)

A073904 Smallest multiple k*n of n having n divisors.

Original entry on oeis.org

1, 2, 9, 8, 625, 12, 117649, 24, 36, 80, 25937424601, 60, 23298085122481, 448, 2025, 384, 48661191875666868481, 180, 104127350297911241532841, 240, 35721, 11264, 907846434775996175406740561329, 360, 10000, 53248, 26244, 1344
Offset: 1

Views

Author

Amarnath Murthy, Aug 18 2002

Keywords

Comments

Smallest refactorable number, m, such that m=k*n has n divisors. - Robert G. Wilson v, Oct 31 2005

Examples

			Smallest multiple a(n)=k*n; a(1)=1*1, a(2)=1*2, a(3)=3*3, a(4)=2*4, a(5)=125*5, a(6)=2*6, ... having d(k*n)=n divisors; d(1)=1, d(2)=2, d(3^2)=3, d(2^3)=4, d(5^4)=5, d(2^2*3)=3*2=6, ...

Links

Crossrefs

Cf. A076931, A050376, A005179, A037992, A050376, A111172, A299795.
Cf. A033950 (refactorable numbers, also known as tau numbers).
Cf. A110821 (SuperRefactorable numbers).

Programs

  • Mathematica
    f[n_] := Block[{k = 1}, If[ PrimeQ[n], n^(n - 1), While[d = DivisorSigma[0, k*n]; d != n, k++ ]; k*n]]; Table[ f[n], {n, 28}] (* Robert G. Wilson v *)

Formula

If p is a prime then a(p) = p^(p-1). If n = p^2 then a(n) = 2^(p-1)*p^(p-1).
a(p^r) = (2*3*5*...*p_r)^(p-1) for r < p <= p_r. a(p^r) = (2*3*...*p_(r-1))^(p-1)*p^(p-1) for p > p_r. Else a(p^r) = ...? for r >= p. Problem a(2^r) = ...? Cf. A005179(p^n)=(2*3*...*p_n)^(p-1) for p_n < 2^p. - Thomas Ordowski, Aug 20 2005
a(p^r) = (2*3...*p_(r-1)*p)^(p-1) for p > p_r; else a(p^r) = (2*3...*p...*p_m)^(p-1)*p^(p^k-p) for p <= p_r and p_m < 2^p, where m=r-k+1 for smallest k such that p^k > r, so k=floor(log(r)/log(p))+1 and p > log(p_m)/log(2). Examples: If k=1 then a(p^r) = (2*3*...*p_r)^(p-1) for r < p <= p_r. If p=2 then a(2^r) = (2*3*...*p_m)*2^(2^k-2) for r < 5. For instance, let r=4 so k=3, m=2 and a(2^4)=384. - Thomas Ordowski, Aug 22 2005
If p is a prime and n=p^r then a(p^r) = (s_1*s_2*...*s_r)^(p-1) where (s_r) is a permutation of the (ascending sequence) numbers of the form q^(p^j) for every prime q and j>=0; permutation such that s_(p^j)=p^(p^j) and shifted remainder. For example, if p=3 then (s_r): 3, 2, 3^3, 5, 7, 2^3, 11, 13, 3^9, 17, 19, ... so a(3^r) = (3*2*27*5*...*s_r)^2. - Thomas Ordowski, Aug 29 2005
If n=2^r then a(2^r) is the product of the first r members of the A109429 sequence. - Thomas Ordowski, Aug 29 2005
a(n) = n * A076931(n). - Thomas Ordowski, Oct 07 2005
a(4) = 8; a(2*prime(n)) = A299795(n), for n>1. - Bernard Schott, Nov 06 2022

Extensions

a(12) corrected by Thomas Ordowski, Aug 18 2005
Further corrections from Thomas Ordowski, Oct 07 2005
a(21), a(27) & a(28) from Robert G. Wilson v, Oct 31 2005

A309141 Nonunitary highly composite numbers: numbers with a record number of nonunitary divisors.

Original entry on oeis.org

1, 4, 8, 16, 24, 36, 48, 72, 144, 216, 288, 360, 576, 720, 1080, 1440, 2160, 2880, 3600, 4320, 5040, 7200, 7560, 10080, 15120, 20160, 25200, 30240, 45360, 50400, 60480, 75600, 100800, 110880, 151200, 221760, 277200, 302400, 332640, 453600, 498960, 554400, 665280
Offset: 1

Views

Author

Amiram Eldar, Jul 14 2019

Keywords

Comments

Numbers k with A048105(k) > A048105(j) for all j < k.
The corresponding values of records are 0, 1, 2, 3, 4, 5, 6, 8, 11, 12, 14, 16, 17, 22, 24, 28, 32, 34, 37, 40, 44, 46, 48, 56, 64, 68, 74, 80, 84, 92, 96, ... (see the link for more values)

Links

Crossrefs

Cf. A048105, A002182 (highly composite), A002110 (unitary), A037992 (infinitary), A293185 (bi-unitary), A318278 (exponential).

Programs

  • Mathematica
    f[n_] := DivisorSigma[0, n] - 2^PrimeNu[n]; fm=-1; s={}; Do[f1 = f[n]; If[f1 > fm, fm = f1; AppendTo[s, n]], {n, 1, 10^5}]; s

A016017 Smallest k such that 1/k can be written as a sum of exactly 2 unit fractions in n ways.

Original entry on oeis.org

1, 2, 4, 8, 6, 32, 64, 12, 256, 512, 24, 2048, 36, 30, 16384, 32768, 96, 72, 262144, 192, 1048576, 2097152, 60, 8388608, 216, 768, 67108864, 288, 1536, 536870912, 1073741824, 120, 576, 8589934592, 6144, 34359738368, 68719476736, 180, 864
Offset: 1

Views

Author

Robert G. Wilson v

Keywords

Comments

From Jianing Song, Aug 30 2021: (Start)
a(n) is the smallest number whose square has exactly 2n-1 divisors.
a(n) is the earliest occurrence of 2n-1 in A048691. (End)

Examples

			a(1)=1 and a(2)=2 because 1/2 = 1/3 + 1/6 = 1/4 + 1/4.
a(3)=4 because 1/4 = 1/5 + 1/20 = 1/6 + 1/12 = 1/8 + 1/8.
a(4)=8 because 1/8 = 1/9 + 1/72 = 1/10 + 1/40 = 1/12 + 1/24 = 1/16 + 1/16.
a(5)=6 because 1/6 = 1/7 + 1/42 = 1/8 + 1/24 = 1/9 + 1/18 = 1/10 + 1/15 = 1/12 + 1/12.

Links

Crossrefs

Identical to A071571 shifted right.
Cf. A000005, A000290, A005408, A005179, A003680, A037992, A055079, A048691.

Programs

  • Mathematica
    f[j_, n_] := (Times @@ (j(Last /@ FactorInteger[n]) + 1) + j - 1)/j; t = Table[0, {50}]; Do[a = f[2, n]; If[a < 51 && t[[a]] == 0, t[[a]] = n; Print[{a, n}]], {n, 2^30}] (* Robert G. Wilson v, Aug 03 2005 *)
  • PARI
    a(n) = {k = 1; while (numdiv(k^2) != (2*n-1), k++); return (k); }; \\ Amiram Eldar, Jan 07 2019 after Michel Marcus at A071571

Formula

a(n+1) <= 2^n.
From Labos Elemer, May 22 2001: (Start)
a(n) = sqrt(A061283(n)).
a(n) = sqrt(Min{k| A000005(k)=2n-1}).
a((p+1)/2) = 2^((p-1)/2) = 2^A005097(i) if p is the i-th odd prime. [Corrected by Jianing Song, Aug 30 2021] (End)
a(n) is the least k such that (tau(k^2) + 1)/2 = n. - Vladeta Jovovic, Aug 01 2001

Extensions

Entry revised by N. J. A. Sloane, Aug 14 2005
Offset corrected by David W. Wilson, Dec 27 2018

A306736 Exponential infinitary highly composite numbers: where the number of exponential infinitary divisors (A307848) increases to record.

Original entry on oeis.org

1, 4, 36, 576, 14400, 705600, 57153600, 6915585600, 1168733966400, 337764116289600, 121932845980545600, 64502475523708622400, 40314047202317889000000, 33904113697149344649000000, 32581853262960520207689000000, 44604557116992952164326241000000, 74980260513665152588232411121000000
Offset: 1

Views

Author

Amiram Eldar, May 01 2019

Keywords

Comments

Subsequence of A025487.
All the terms have prime factors with multiplicities which are infinitary highly composite number (A037992) > 1, similarly to exponential highly composite numbers (A318278) whose prime factors have multiplicities which are highly composite numbers (A002182). Thus all the terms are squares. Their square roots are 1, 2, 6, 24, 120, 840, 7560, 83160, 1081080, 18378360, 349188840, 8031343320, 200783583000, 5822723907000, 180504441117000, ...
Differs from A307845 (exponential unitary highly composite numbers) from n >= 107. a(107) = 2^24 * (3 * 5 * ... * 19)^6 * (23 * 29 * ... * 509)^2 ~ 2.370804... * 10^456, while A307845(107) = (2 * 3 * 5 * ... * 19)^6 * (23 * 29 * ... * 521)^2 ~ 2.454885... * 10^456.

Links

Crossrefs

Cf. A002182, A025487, A037445, A037992, A307845, A307848, A318278.

Programs

  • Mathematica
    di[1] = 1; di[n_] := Times @@ Flatten[2^DigitCount[#, 2, 1] & /@ FactorInteger[n][[All, 2]]]; fun[p_, e_] := di[e]; a[1] = 1; a[n_] := Times @@ (fun @@@ FactorInteger[n]); s = {}; am = 0; Do[a1 = a[n]; If[a1 > am, am = a1; AppendTo[s, n]], {n, 1, 10^6}]; s (* after Jean-François Alcover at A037445 *)

Formula

A307848(a(n)) = 2^(n-1).
Showing 1-10 of 44 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.