A096932 -id:A096932 - cp's OEIS frontend

A001358 Semiprimes (or biprimes): products of two primes.

4, 6, 9, 10, 14, 15, 21, 22, 25, 26, 33, 34, 35, 38, 39, 46, 49, 51, 55, 57, 58, 62, 65, 69, 74, 77, 82, 85, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 121, 122, 123, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166, 169, 177, 178, 183, 185, 187

Offset: 1

N. J. A. Sloane, R. K. Guy

Numbers of the form p*q where p and q are primes, not necessarily distinct.
These numbers are sometimes called semiprimes or 2-almost primes.
Numbers n such that Omega(n) = 2 where Omega(n) = A001222(n) is the sum of the exponents in the prime decomposition of n.
Complement of A100959; A064911(a(n)) = 1. - Reinhard Zumkeller, Nov 22 2004
The graph of this sequence appears to be a straight line with slope 4. However, the asymptotic formula shows that the linearity is an illusion and in fact a(n)/n ~ log(n)/log(log(n)) goes to infinity. See also the graph of A066265 = number of semiprimes < 10^n.
For numbers between 33 and 15495, semiprimes are more plentiful than any other k-almost prime. See A125149.
Numbers that are divisible by exactly 2 prime powers (not including 1). - Jason Kimberley, Oct 02 2011
The (disjoint) union of A006881 and A001248. - Jason Kimberley, Nov 11 2015
An equivalent definition of this sequence is a'(n) = smallest composite number which is not divided by any smaller composite number a'(1),...,a'(n-1). - Meir-Simchah Panzer, Jun 22 2016
The above characterization can be simplified to "Composite numbers not divisible by a smaller term." This shows that this is the equivalent of primes computed via Eratosthenes's sieve, but starting with the set of composite numbers (i.e., complement of 1 union primes) instead of all positive integers > 1. It's easy to see that iterating the method (using Eratosthenes's sieve each time on the remaining numbers, complement of the previously computed set) yields numbers with bigomega = k for k = 0, 1, 2, 3, ..., i.e., {1}, A000040, this, A014612, etc. - M. F. Hasler, Apr 24 2019
For all n except n = 2, a(n) is a deficient number. - Amrit Awasthi, Sep 10 2024
It is reasonable to assume that the "comforting numbers" which John T. Williams found in Chapter 3 of Milne's book "The House at Pooh Corner" are these semiprimes. Winnie-the-Pooh wonders whether he has 14 or 15 honey pots and concludes: "It's sort of comforting." To arrange a semiprime number of honey pots in a rectangular way, let's say on a shelf, with the larger divisor parallel to the wall, there is only one solution and this is for a simple mind like Winnie-the-Pooh comforting. - Ruediger Jehn, Dec 12 2024

			From _Gus Wiseman_, May 27 2021: (Start)
The sequence of terms together with their prime factors begins:
   4 = 2*2     46 = 2*23     91 = 7*13    141 = 3*47
   6 = 2*3     49 = 7*7      93 = 3*31    142 = 2*71
   9 = 3*3     51 = 3*17     94 = 2*47    143 = 11*13
  10 = 2*5     55 = 5*11     95 = 5*19    145 = 5*29
  14 = 2*7     57 = 3*19    106 = 2*53    146 = 2*73
  15 = 3*5     58 = 2*29    111 = 3*37    155 = 5*31
  21 = 3*7     62 = 2*31    115 = 5*23    158 = 2*79
  22 = 2*11    65 = 5*13    118 = 2*59    159 = 3*53
  25 = 5*5     69 = 3*23    119 = 7*17    161 = 7*23
  26 = 2*13    74 = 2*37    121 = 11*11   166 = 2*83
  33 = 3*11    77 = 7*11    122 = 2*61    169 = 13*13
  34 = 2*17    82 = 2*41    123 = 3*41    177 = 3*59
  35 = 5*7     85 = 5*17    129 = 3*43    178 = 2*89
  38 = 2*19    86 = 2*43    133 = 7*19    183 = 3*61
  39 = 3*13    87 = 3*29    134 = 2*67    185 = 5*37
(End)

Archimedeans Problems Drive, Eureka, 17 (1954), 8.
Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; Chapter II, Problem 60.
Edmund Landau, Handbuch der Lehre von der Verteilung der Primzahlen, Vol. 1, Teubner, Leipzig; third edition: Chelsea, New York (1974). See p. 211.
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).
John T. Williams, Pooh and the Philosophers, Dutton Books, 1995.

N. J. A. Sloane, Table of n, a(n) for n = 1..20000 (first 10000 terms from T. D. Noe)
Dragos Crisan and Radek Erban, On the counting function of semiprimes, INTEGERS, Vol. 21 (2021), #A122.
Daniel A. Goldston, Sidney W. Graham, János Pintz and Cem Y. Yildirim, Small gaps between primes or almost primes, Transactions of the American Mathematical Society, Vol. 361, No. 10 (2009), pp. 5285-5330, arXiv preprint, arXiv:math/0506067 [math.NT], 2005.
Richard K. Guy, Letters to N. J. A. Sloane, June-August 1968
Sh. T. Ishmukhametov and F. F. Sharifullina, On distribution of semiprime numbers, Izvestiya Vysshikh Uchebnykh Zavedenii. Matematika, 2014, No. 8, pp. 53-59. English translation, Russian Mathematics, Vol. 58, No. 8 (2014), pp. 43-48, alternative link.
Donovan Johnson, Jonathan Vos Post, and Robert G. Wilson v, Selected n and a(n). (2.5 MB)
Dixon Jones, Quickie 593, Mathematics Magazine, Vol. 47, No. 3, May 1974, p. 167.
Edmund Landau, Handbuch der Lehre von der Verteilung der Primzahlen, vol. 1 and vol. 2, Leipzig, Berlin, B. G. Teubner, 1909. See Vol. 1, p. 211.
Xianmeng Meng, On sums of three integers with a fixed number of prime factors, Journal of Number Theory, Vol. 114, No. 1 (2005), pp. 37-65.
Michael Penn, What makes a number "good"?, YouTube video, 2022.
Carlos Rivera, Conjecture 108: On the counting function of semiprimes, The Prime Puzzles & Problems Connection.
Eric Weisstein's World of Mathematics, Semiprime.
Eric Weisstein's World of Mathematics, Almost Prime.
Wikipedia, Almost prime.
Robert G. Wilson v, Subsequences at various powers of 10.
Index to sequences related to sums of cubes
Index entries for "core" sequences

Cf. A064911 (characteristic function).
Cf. A048623, A048639, A000040 (primes), A014612 (products of 3 primes), A014613, A014614, A072000 ("pi" for semiprimes), A065516 (first differences).
Cf. A077554, A077555, A002024, A072966, A100592, A014673, A068318, A061299, A087718, A089994, A089995, A096916, A096932, A106550, A106554, A108541, A108542, A126663, A131284, A138510, A138511, A072931, A088183, A171963, A237040 (semiprimes of form n^3 + 1).
Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r=1), this sequence (r=2), A014612 (r=3), A014613 (r=4), A014614 (r=5), A046306 (r=6), A046308 (r=7), A046310 (r=8), A046312 (r=9), A046314 (r=10), A069272 (r=11), A069273 (r=12), A069274 (r=13), A069275 (r=14), A069276 (r=15), A069277 (r=16), A069278 (r=17), A069279 (r=18), A069280 (r=19), A069281 (r=20).
These are the Heinz numbers of length-2 partitions, counted by A004526.
The squarefree case is A006881 with odd/even terms A046388/A100484 (except 4).
Including primes gives A037143.
The odd/even terms are A046315/A100484.
Partial sums are A062198.
The prime factors are A084126/A084127.
Grouping by greater factor gives A087112.
The product/sum/difference of prime indices is A087794/A176504/A176506.
Positions of even/odd terms are A115392/A289182.
The terms with relatively prime/divisible prime indices are A300912/A318990.
Factorizations using these terms are counted by A320655.
The prime indices are A338898/A338912/A338913.
Grouping by weight (sum of prime indices) gives A338904, with row sums A024697.
The terms with even/odd weight are A338906/A338907.
The terms with odd/even prime indices are A338910/A338911.
The least/greatest term of weight n is A339114/A339115.
Cf. A014342, A098350, A112141, A320732, A332765, A339003/A339004, A339116.

a001358 n = a001358_list !! (n-1)
a001358_list = filter ((== 2) . a001222) [1..]

[n: n in [2..200] | &+[d[2]: d in Factorization(n)] eq 2]; // Bruno Berselli, Sep 09 2015

A001358 := proc(n) option remember; local a; if n = 1 then 4; else for a from procname(n-1)+1 do if numtheory[bigomega](a) = 2 then return a; end if; end do: end if; end proc:
seq(A001358(n), n=1..120) ; # R. J. Mathar, Aug 12 2010

Select[Range[200], Plus@@Last/@FactorInteger[#] == 2 &] (* Zak Seidov, Jun 14 2005 *)
Select[Range[200], PrimeOmega[#]==2&] (* Harvey P. Dale, Jul 17 2011 *)

select( isA001358(n)={bigomega(n)==2}, [1..199]) \\ M. F. Hasler, Apr 09 2008; added select() Apr 24 2019

list(lim)=my(v=List(),t);forprime(p=2, sqrt(lim), t=p;forprime(q=p, lim\t, listput(v,t*q))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Sep 11 2011

A1358=List(4); A001358(n)={while(#A1358M. F. Hasler, Apr 24 2019

from sympy import factorint
def ok(n): return sum(factorint(n).values()) == 2
print([k for k in range(1, 190) if ok(k)]) # Michael S. Branicky, Apr 30 2022

from math import isqrt
from sympy import primepi, prime
def A001358(n):
    def f(x): return int(n+x-sum(primepi(x//prime(k))-k+1 for k in range(1, primepi(isqrt(x))+1)))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m # Chai Wah Wu, Jul 23 2024

a(n) ~ n*log(n)/log(log(n)) as n -> infinity [Landau, p. 211], [Ayoub].
Recurrence: a(1) = 4; for n > 1, a(n) = smallest composite number which is not a multiple of any of the previous terms. - Amarnath Murthy, Nov 10 2002
A174956(a(n)) = n. - Reinhard Zumkeller, Apr 03 2010
a(n) = A088707(n) - 1. - Reinhard Zumkeller, Feb 20 2012
Sum_{n>=1} 1/a(n)^s = (1/2)*(P(s)^2 + P(2*s)), where P is the prime zeta function. - Enrique Pérez Herrero, Jun 24 2012
sigma(a(n)) + phi(a(n)) - mu(a(n)) = 2*a(n) + 1. mu(a(n)) = ceiling(sqrt(a(n))) - floor(sqrt(a(n))). - Wesley Ivan Hurt, May 21 2013
mu(a(n)) = -Omega(a(n)) + omega(a(n)) + 1, where mu is the Moebius function (A008683), Omega is the count of prime factors with repetition, and omega is the count of distinct prime factors. - Alonso del Arte, May 09 2014
a(n) = A078840(2,n). - R. J. Mathar, Jan 30 2019
A100484 UNION A046315. - R. J. Mathar, Apr 19 2023
Conjecture: a(n)/n ~ (log(n)/log(log(n)))*(1-(M/log(log(n)))) as n -> oo, where M is the Mertens's constant (A077761). - Alain Rocchelli, Feb 02 2025

More terms from James Sellers, Aug 22 2000

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

N. J. A. Sloane, David Singmaster

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]

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.

Matthew House, Table of n, a(n) for n = 1..3322 (terms 1..2000 from Don Reble)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Ron Brown, The minimal number with a given number of divisors, Journal of Number Theory 116 (2006) 150-158.
M. E. Grost, The smallest number with a given number of divisors, Amer. Math. Monthly, 75 (1968), 725-729.
Joe Roberts, Lure of the Integers, Annotated scanned copy of pp. 81, 86 with notes.
Anna K. Savvopoulou and Christopher M. Wedrychowicz, On the smallest number with a given number of divisors, The Ramanujan Journal, 2015, Vol. 37, pp. 51-64.
David Singmaster, Letter to N. J. A. Sloane, Oct 03 1982.
T. Verhoeff, Rectangular and Trapezoidal Arrangements, J. Integer Sequences, Vol. 2 (1999), Article 99.1.6.
Eric Weisstein's World of Mathematics, Divisor.
Robert G. Wilson v, Letter to N. J. A. Sloane, Dec 17 1991.

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

import Data.List (elemIndex)
import Data.Maybe (fromJust)
a005179 n = succ $ fromJust $ elemIndex n $ map a000005 [1..]
-- Reinhard Zumkeller, Apr 01 2011

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

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

(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

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

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

More terms from David W. Wilson

A084127 Prime factor >= other prime factor of n-th semiprime.

Original entry on oeis.org

2, 3, 3, 5, 7, 5, 7, 11, 5, 13, 11, 17, 7, 19, 13, 23, 7, 17, 11, 19, 29, 31, 13, 23, 37, 11, 41, 17, 43, 29, 13, 31, 47, 19, 53, 37, 23, 59, 17, 11, 61, 41, 43, 19, 67, 47, 71, 13, 29, 73, 31, 79, 53, 23, 83, 13, 59, 89, 61, 37, 17, 97, 67, 101, 29, 41, 103, 19, 71, 107, 43, 31

Offset: 1

Reinhard Zumkeller, May 15 2003

nonn

Largest nontrivial divisor of n-th semiprime. [Juri-Stepan Gerasimov, Apr 18 2010]

Greater of the prime factors of A001358(n). - Jianing Song, Aug 05 2022

Zak Seidov, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Semiprime

Cf. A001358 (the semiprimes), A084126 (lesser of the prime factors of the semiprimes).

Cf. A014673, A061299, A068318, A087718, A087794, A089994, A089995, A096932, A106550, A106554, A108542, A126663, A131284, A138510, A138511.

a084127 = a006530 . a001358  -- Reinhard Zumkeller, Nov 25 2012

FactorInteger[#][[-1, 1]]& /@ Select[Range[1000], PrimeOmega[#] == 2&] (* Jean-François Alcover, Nov 17 2021 *)

lista(nn) = {for (n=2, nn, if (bigomega(n)==2, f = factor(n); print1(f[length(f~),1], ", ")););} \\ Michel Marcus, Jun 05 2013

from math import isqrt
from sympy import primepi, primerange, primefactors
def A084127(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: 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 int(n+x+((t:=primepi(s:=isqrt(x)))*(t-1)>>1)-sum(primepi(x//p) for p in primerange(s+1)))
    return max(primefactors(bisection(f,n,n))) # Chai Wah Wu, Oct 23 2024

a(n) = A006530(A001358(n)).
a(n) = A001358(n)/A020639(A001358(n)). [corrected by Michel Marcus, Jul 18 2020]
a(n) = A001358(n)/A084126(n).

Corrected by T. D. Noe, Nov 15 2006

A137492 Numbers with 29 divisors.

Original entry on oeis.org

268435456, 22876792454961, 37252902984619140625, 459986536544739960976801, 144209936106499234037676064081, 15502932802662396215269535105521, 28351092476867700887730107366063041

Offset: 1

R. J. Mathar, Apr 22 2008

Maple implementation: see A030513.

28th powers of primes. The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. - Omar E. Pol, May 06 2008

T. D. Noe, Table of n, a(n) for n = 1..1000
OEIS Wiki, Index entries for number of divisors

Cf. A030513-A030516, A030626, A030627, A030634-A030638, A005179, A003680, A096932, A061286, A061283.

Cf. A000040.

Prime[Range[13]]^28 (* Vladimir Joseph Stephan Orlovsky, May 05 2011 *)

a(n)=prime(n)^28 \\ Charles R Greathouse IV, Jun 19 2016

A000005(a(n))=29.

a(n)=A000040(n)^(29-1)=A000040(n)^(28). - Omar E. Pol, May 06 2008

A061299 Least number whose number of divisors is A007304(n) (the n-th number that is the product of 3 distinct primes).

Original entry on oeis.org

720, 2880, 46080, 25920, 184320, 2949120, 129600, 414720, 11796480, 1658880, 188743680, 3732480, 2073600, 26542080, 12079595520, 14929920, 48318382080, 106168320, 8294400, 3092376453120, 1698693120, 18662400, 238878720

Offset: 1

Labos Elemer, Jun 05 2001

nonn

All terms are divisible by a(1) = 720, the first entry.

All terms [=a(j)], not only arguments [=j] have 3 distinct prime factors at exponents determined by the p,q,r factors of their arguments: a(pqr) = RPQ.

			For n = 5: A007304(5) = 78 = 2*3*13, A005179(78) = 184320 = (2^12)*(3^2)*(5^1) = a(5).

Amiram Eldar, Table of n, a(n) for n = 1..3706 (terms 1..835 from David A. Corneth)

Cf. A000005, A005179, A007304, A061148, A061149.

Cf. A096932, A061234, A061286.

a(n) = A005179(A007304(n)).
Min{x; A000005(x) = pqr} p, q, r are distinct primes. If k = pqr and p > q > r then A005179(k) = 2^(p-1)*3^(q-1)*5^(r-1).
From Reinhard Zumkeller, Jul 15 2004: (Start)
A000005(a(n)) = A007304(n).
A000005(m) != A007304(n) for m < a(n).
a(n) = A005179(A007304(n)).
a(p*m*q) = 2^(q-1) * 3^(m-1) * 5^(p-1) for primes p < m < q.
a(A000040(i)*A000040(j)*A000040(k)) = 2^(A084127(k)-1) * 3^(A084127(j)-1) * 5^(A084127(i)-1) for i < j < k. (End)

Edited by N. J. A. Sloane, Apr 20 2007

A137488 Numbers with 25 divisors.

Original entry on oeis.org

1296, 10000, 38416, 50625, 194481, 234256, 456976, 1185921, 1336336, 1500625, 2085136, 2313441, 4477456, 6765201, 9150625, 10556001, 11316496, 14776336, 16777216, 17850625, 22667121, 29986576, 35153041, 45212176, 52200625

Offset: 1

R. J. Mathar, Apr 22 2008

nonn

Maple implementation: see A030513.

Numbers of the form p^24 (24th powers of A000040, subset of A010812) or p^4*q^4 (A189991), where p and q are distinct primes. - R. J. Mathar, Mar 01 2010

T. D. Noe, Table of n, a(n) for n = 1..1000

Cf. A000005, A010812, A030513-A030516, A030626, A030627, A030634-A030638, A005179, A003680, A096932, A061286, A061283, A135581, A175755.

a137488 n = a137488_list !! (n-1)
a137488_list = m (map (^ 24) a000040_list) (map (^ 4) a006881_list) where
   m xs'@(x:xs) ys'@(y:ys) | x < y = x : m xs ys'
                           | otherwise = y : m xs' ys
-- Reinhard Zumkeller, Nov 29 2011

lst = {}; Do[If[DivisorSigma[0, n] == 25, Print[n]; AppendTo[lst, n]], {n, 55000000}]; lst (* Vladimir Joseph Stephan Orlovsky, May 03 2011 *)
Select[Range[5221*10^4],DivisorSigma[0,#]==25&] (* Harvey P. Dale, Mar 11 2019 *)

is(n)=numdiv(n)==25 \\ Charles R Greathouse IV, Jun 19 2016

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A137488(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 int(n+x+(t:=primepi(s:=isqrt(y:=integer_nthroot(x,4)[0])))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)))-primepi(integer_nthroot(x,24)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 22 2025

A000005(a(n)) = 25.

Sum_{n>=1} 1/a(n) = (P(4)^2 - P(8))/2 + P(24) = 0.000933328..., where P is the prime zeta function. - Amiram Eldar, Jul 03 2022

cp's OEIS Frontend

A001358 Semiprimes (or biprimes): products of two primes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Python

Formula

Extensions

A005179 Smallest number with exactly n divisors.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A084127 Prime factor >= other prime factor of n-th semiprime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A137492 Numbers with 29 divisors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A061299 Least number whose number of divisors is A007304(n) (the n-th number that is the product of 3 distinct primes).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A137488 Numbers with 25 divisors.

Views

Author

Keywords

Comments

Links

Crossrefs