A182862 -id:A182862 - cp's OEIS frontend

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

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

David W. Wilson

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

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

Franklin T. Adams-Watters, Table of n, a(n) for n = 1..10001 (first 291 terms from Will Nicholes)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972.
Kevin Broughan, Equivalents of the Riemann Hypothesis, Vol. 1: Arithmetic Equivalents, Cambridge University Press, 2017. See section 8.2, "Hardy-Ramanujan Numbers".
YoungJu Choie, Nicolas Lichiardopol, Pieter Moree and Patrick Solé, On Robin's criterion for the Riemann hypothesis, Journal de théorie des nombres de Bordeaux, Vol. 19, No. 2 (2007), pp. 357-372. See section 5, p. 367.
Asaf Cohen Antonir and Asaf Shapira, An Elementary Proof of a Theorem of Hardy and Ramanujan (2022). arXiv:2207.09410 [math.NT]
Michael De Vlieger, Relations of A025487 to A002110, A002182, and A002201.
Steven R. Finch, Errata and Addenda to Mathematical Constants, arXiv:2001.00578 [math.HO], 2020, pp. 9-10.
G. H. Hardy and S. Ramanujan, Asymptotic formulae for the distribution of integers of various types, Proc. London Math. Soc, Ser. 2, Vol. 16 (1917), pp. 112-132. Also published in the book Collected Papers of Srinivasa Ramanujan, Chelsea, 1962, pages 245-261.
Jeffery Kline, On the eigenstructure of sparse matrices related to the prime number theorem, Linear Algebra and its Applications (2020) Vol. 584, 409-430.
L. B. Richmond, Asymptotic results for partitions (I) and the distribution of certain integers, Journal of Number Theory, Vol. 8, No. 4 (1976), pp. 372-389. See page 388.

Subsequence of A055932, A191743, and A324583.
Cf. A025488, A051282, A036041, A051466, A061394, A124832, A161360, A166469, A181815, A181817, A283980, A306802, A322584, A322585 (characteristic function), A329897, A329898, A329899, A329900, A329904, A330683.
Cf. A085089, A101296 (left inverses).
Equals range of values taken by A046523.
Cf. A178799 (first differences), A247451 (squarefree kernel), A146288 (number of divisors).
Subsequences of this sequence include: A000079, A000142, A000400, A001013, A001813, A002110, A002182, A005179, A006939, A025527, A056836, A061742, A064350, A066120, A087980, A097212, A097213, A111059, A119840, A119845, A126098, A129912, A140999, A166338, A166470, A166472, A166473, A166475, A167448, A168262, A168263, A168264, A179215, A181555, A181804, A181806, A181809, A181818, A181822, A181823, A181824, A181825, A181826, A181827, A182763, A182862, A182863, A212170, A220264, A220423, A250269, A250270, A260633, A266047, A284456, A300357, A304938, A329894, A330687; also A037019 and A330681 (when sorted), possibly also A289132.
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.

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

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

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

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

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

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

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

\\ 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

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

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

Offset corrected by Matthew Vandermast, Oct 19 2008

Minor correction by Charles R Greathouse IV, Sep 03 2010

A212167 Numbers k such that the maximum exponent in its prime factorization is not greater than the number of positive exponents (A051903(k) <= A001221(k)).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83

Offset: 1

Matthew Vandermast, May 22 2012

nonn

Union of A212166 and A212168. Includes numerous subsequences that are subsequences of neither A212166 nor A212168.

			40 = 2^3*5^1 has 2 distinct prime factors, hence, 2 positive exponents in its prime factorization (although the 1 is often left implicit).  2 is less than the maximal exponent in 40's prime factorization, which is 3. Therefore, 40 does not belong to the sequence. But 10 = 2^1*5^1 and 20 = 2^2*5^1 belong, since the maximal exponents in their prime factorizations are 1 and 2 respectively.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
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].
Primefan, The First 2500 Integers Factored (first of 5 pages).

Complement of A212164. See also A212165.
Subsequences (none of which are subsequences of A212166 or A212168) include A002110, A051451, A129912, A179983, A181826, A181827, A182862, A182863. Includes all members of A003418.
Cf. A001221, A050326, A051903, A188654, A225230.

import Data.List (findIndices)
a212167 n = a212167_list !! (n-1)
a212167_list = map (+ 1) $ findIndices (>= 0) a225230_list
-- Reinhard Zumkeller, May 03 2013

isA212167 := proc(n)
    simplify(A051903(n) <= A001221(n)) ;
end proc:
for n from 1 to 1000 do
    if isA212167(n) then
        printf("%d,",n) ;
    end if;
end do: # R. J. Mathar, Jan 06 2021

okQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Max[f] <= Length[f]]; Select[Range[1000], okQ] (* T. D. Noe, May 24 2012 *)

is(k) = {my(e = factor(k)[, 2]); !(#e) || vecmax(e) <= #e; } \\ Amiram Eldar, Sep 09 2024

A225230(a(n)) >= 0; A050326(a(n)) > 0. - Reinhard Zumkeller, May 03 2013

A182860 Number of distinct prime signatures represented among the unitary divisors of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 3, 2, 2, 2, 3, 2, 4, 2, 3, 3, 2, 2, 4, 2, 4, 3, 3, 2, 4, 2, 3, 2, 4, 2, 4, 2, 2, 3, 3, 3, 3, 2, 3, 3, 4, 2, 4, 2, 4, 4, 3, 2, 4, 2, 4, 3, 4, 2, 4, 3, 4, 3, 3, 2, 6, 2, 3, 4, 2, 3, 4, 2, 4, 3, 4, 2, 4, 2, 3, 4, 4, 3, 4, 2, 4, 2, 3, 2, 6, 3, 3, 3, 4, 2, 6, 3, 4, 3, 3, 3, 4, 2, 4, 4, 3, 2, 4, 2, 4, 4

Offset: 1

Matthew Vandermast, Jan 14 2011

nonn

a(n) = number of members m of A025487 such that d(m^k) divides d(n^k) for all values of k. (Here d(n) represents the number of divisors of n, or A000005(n).)

a(n) depends only on prime signature of n (cf. A025487).

			60 has 8 unitary divisors (1, 3, 4, 5, 12, 15, 20 and 60). Primes 3 and 5 have the same prime signature, as do 12 (2^2*3) and 20 (2^2*5); each of the other four numbers listed is the only unitary divisor of 60 with its particular prime signature.  Since a total of 6 distinct prime signatures appear among the unitary divisors of 60, a(60) = 6.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Unitary Divisor
Index entries for sequences computed from exponents in factorization of n

Cf. A000005, A034444, A064553, A085082, A146289, A146290, A182861, A182862, A212180, A328830.

Table[Length@ Union@ Map[Sort[FactorInteger[#] /. {p_, e_} /; p > 0 :> If[p == 1, 0, e]] &, Select[Divisors@ n, CoprimeQ[#, n/#] &]], {n, 105}] (* Michael De Vlieger, Jul 19 2017 *)

A181819(n) = {my(f=factor(n)); prod(k=1, #f~, prime(f[k, 2])); }; \\ From A181819
A182860(n) = numdiv(A181819(n)); \\ Antti Karttunen, Jul 19 2017

a(n) = A000005(A181819(n)).
If the canonical factorization of n into prime powers is Product p^e(p), then the formula for d(n^k) is Product_p (ek + 1). (See also A146289, A146290.)
a(n) = A064553(A328830(n)). - Antti Karttunen, Apr 29 2022

A182863 Members m of A025487 such that, if k appears in m's prime signature, k-1 appears at least as often as k (for any integer k > 1).

Original entry on oeis.org

1, 2, 6, 12, 30, 60, 210, 360, 420, 1260, 2310, 2520, 4620, 13860, 27720, 30030, 60060, 75600, 138600, 180180, 360360, 510510, 831600, 900900, 1021020, 1801800, 3063060, 6126120, 9699690, 10810800, 15315300, 19399380, 30630600, 37837800

Offset: 1

Matthew Vandermast, Jan 14 2011

nonn

Members m of A025487 such that A181819(m) is also a member of A025487.
If prime signatures are considered as partitions, these are the members of A025487 whose prime signature is conjugate to the prime signature of a member of A181818.
Also the least number with each sorted prime metasignature, where a number's metasignature is the sequence of multiplicities of exponents in its prime factorization. For example, 2520 has prime indices {1,1,1,2,2,3,4}, sorted prime signature {1,1,2,3}, and sorted prime metasignature {1,1,2}. - Gus Wiseman, May 21 2022

			The prime signature of 360360 = 2^3*3^2*5*7*11*13 is (3,2,1,1,1,1). 2 appears as many times as 3 in 360360's prime signature, and 1 appears more times than 2. Since 360360 is also a member of A025487, it is a member of this sequence.
From _Gus Wiseman_, May 21 2022: (Start)
The terms together with their sorted prime signatures and sorted prime metasignatures begin:
      1: {}                -> {}            -> {}
      2: {1}               -> {1}           -> {1}
      6: {1,2}             -> {1,1}         -> {2}
     12: {1,1,2}           -> {1,2}         -> {1,1}
     30: {1,2,3}           -> {1,1,1}       -> {3}
     60: {1,1,2,3}         -> {1,1,2}       -> {1,2}
    210: {1,2,3,4}         -> {1,1,1,1}     -> {4}
    360: {1,1,1,2,2,3}     -> {1,2,3}       -> {1,1,1}
    420: {1,1,2,3,4}       -> {1,1,1,2}     -> {1,3}
   1260: {1,1,2,2,3,4}     -> {1,1,2,2}     -> {2,2}
   2310: {1,2,3,4,5}       -> {1,1,1,1,1}   -> {5}
   2520: {1,1,1,2,2,3,4}   -> {1,1,2,3}     -> {1,1,2}
   4620: {1,1,2,3,4,5}     -> {1,1,1,1,2}   -> {1,4}
  13860: {1,1,2,2,3,4,5}   -> {1,1,1,2,2}   -> {2,3}
  27720: {1,1,1,2,2,3,4,5} -> {1,1,1,2,3}   -> {1,1,3}
  30030: {1,2,3,4,5,6}     -> {1,1,1,1,1,1} -> {6}
  60060: {1,1,2,3,4,5,6}   -> {1,1,1,1,1,2} -> {1,5}
(End)

David A. Corneth, Table of n, a(n) for n = 1..10000 (first 1444 terms from Amiram Eldar)
Eric Weisstein's World of Mathematics, Conjugate Partition.

Intersection of A025487 and A179983.
Subsequence of A129912 and A181826.
Includes all members of A182862.
Positions of first appearances in A353742, unordered version A238747.
A001222 counts prime factors with multiplicity, distinct A001221.
A003963 gives product of prime indices.
A005361 gives product of prime signature, firsts A353500 (sorted A085629).
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
A182850 gives frequency depth of prime indices, counted by A225485.
A323014 gives adjusted frequency depth of prime indices, counted by A325280.
Cf. A000040, A055932, A070175, A097318, A304678, A325238, A353507, A353745.

nn=1000;
r=Table[Sort[Length/@Split[Sort[Last/@If[n==1,{},FactorInteger[n]]]]],{n,nn}];
Select[Range[nn],!MemberQ[Take[r,#-1],r[[#]]]&] (* Gus Wiseman, May 21 2022 *)

A179983 Positive integers m such that, if k appears in m's prime signature, k-1 appears at least as often as k (for any integer k > 1).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83, 84, 85, 86, 87, 89, 90

Offset: 1

Matthew Vandermast, Jan 15 2011

nonn

Numbers m such that A181819(m) is a term of A025487.

			The prime signature of 20 = 2^2*5 is (2,1). Since the largest number appearing in 20's prime signature is 2, and 1 appears as many times as 2, 20 is a member of this sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Index to sequences related to prime signature.

Includes all squarefree numbers (A005117); also includes all members of A054753, A085987, A163569, A182862, A182863.

Cf. A025487, A181819.

isA179983 := proc(n)
    local es,me,k ;
    # list of exponents in prime signature
    es := [seq(op(2,pe), pe =ifactors(n)[2])] ;
    # maximum exponent
    me := max(op(es)) ;
    for k from me to 2 by -1 do
        if numboccur(es,k-1) < numboccur(es,k) then
            return false;
        end if;
    end do:
    true ;
end proc:
for n from 1 to 100 do
    if isA179983(n) then
        printf("%d,",n) ;
    end if;
end do: # R. J. Mathar, Mar 21 2023

q[n_] := Module[{t = SortBy[Tally[FactorInteger[n][[;; , 2]]], First], t1, t2}, t1 = t[[;; , 1]]; t2 = t[[;; , 2]]; Sort[t1] == Range[Length[t1]] && Max[Differences[t2]] < 1]; Select[Range[100], q] (* Amiram Eldar, Aug 04 2024 *)

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Sage

Formula

Extensions

A212167 Numbers k such that the maximum exponent in its prime factorization is not greater than the number of positive exponents (A051903(k) <= A001221(k)).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A182860 Number of distinct prime signatures represented among the unitary divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A182863 Members m of A025487 such that, if k appears in m's prime signature, k-1 appears at least as often as k (for any integer k > 1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A179983 Positive integers m such that, if k appears in m's prime signature, k-1 appears at least as often as k (for any integer k > 1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A182861 Number of distinct prime signatures represented among the unitary divisors of A025487(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A306491 Numbers that set a record for the number of distinct prime signatures represented among their divisors (A085082).