A111059 -id:A111059 - 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

A209862 Permutation of nonnegative integers which maps A209642 into ascending order (A209641).

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 11, 13, 14, 15, 16, 17, 18, 20, 24, 19, 21, 25, 22, 26, 28, 23, 27, 29, 30, 31, 32, 33, 34, 36, 40, 48, 35, 37, 41, 49, 38, 42, 50, 44, 52, 56, 39, 43, 51, 45, 53, 57, 46, 54, 58, 60, 47, 55, 59, 61, 62, 63, 64, 65, 66, 68, 72, 80, 96, 67, 69, 73, 81, 97, 70, 74, 82, 98, 76, 84, 100, 88, 104, 112, 71, 75, 83

Offset: 0

Antti Karttunen, Mar 24 2012

nonn

Conjecture: For all n, a(A054429(n)) = A054429(a(n)), i.e. A054429 acts as a homomorphism (automorphism) of the cyclic group generated by this permutation. This implies also a weaker conjecture given in A209860.
From Gus Wiseman, Aug 24 2021: (Start)
As a triangle with row lengths 2^n, T(n,k) for n > 0 appears (verified up to n = 2^15) to be the unique nonnegative integer whose binary indices are the k-th subset of {1..n} containing n. Here, a binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion, and sets are sorted first by length, then lexicographically. For example, the triangle begins:
   1
   2  3
   4  5  6  7
   8  9 10 12 11 13 14 15
  16 17 18 20 24 19 21 25 22 26 28 23 27 29 30 31
Mathematica: Table[Total[2^(Append[#,n]-1)]&/@Subsets[Range[n-1]],{n,5}]
Row lengths are A000079 (shifted right). Also Column k = 1.
Row sums are A010036.
Using reverse-lexicographic order gives A059893.
Using lexicographic order gives A059894.
Taking binary indices to prime indices gives A339195 (or A019565).
The ordering of sets is A344084.
A version using Heinz numbers is A344085.
(End)

			From _Gus Wiseman_, Aug 24 2021: (Start)
The terms, their binary expansions, and their binary indices begin:
   0:      ~ {}
   1:    1 ~ {1}
   2:   10 ~ {2}
   3:   11 ~ {1,2}
   4:  100 ~ {3}
   5:  101 ~ {1,3}
   6:  110 ~ {2,3}
   7:  111 ~ {1,2,3}
   8: 1000 ~ {4}
   9: 1001 ~ {1,4}
  10: 1010 ~ {2,4}
  12: 1100 ~ {3,4}
  11: 1011 ~ {1,2,4}
  13: 1101 ~ {1,3,4}
  14: 1110 ~ {2,3,4}
  15: 1111 ~ {1,2,3,4}
(End)

Antti Karttunen, Table of n, a(n) for n = 0..32767
Index entries for sequences that are permutations of the natural numbers

Inverse permutation: A209861. Cf. A209860, A209863, A209864, A209865, A209866, A209867, A209868.

Cf. A010036, A026793, A048793, A111059, A147655, A246688, A246867, A261144, A272020, A294648, A339360.

a(n) = A209859(A036044(A209641(n))) = A209859(A056539(A209641(n))).

A112141 Product of the first n semiprimes.

Original entry on oeis.org

4, 24, 216, 2160, 30240, 453600, 9525600, 209563200, 5239080000, 136216080000, 4495130640000, 152834441760000, 5349205461600000, 203269807540800000, 7927522494091200000, 364666034728195200000, 17868635701681564800000, 911300420785759804800000

Offset: 1

Jonathan Vos Post, Nov 28 2005

Semiprime analog of primorial (A002110). Equivalent for product of what A062198 is for sum.

			a(10) = 4*6*9*10*14*15*21*22*25*26 = 136216080000, the product of the first 10 semiprimes.
From _Gus Wiseman_, Dec 06 2020: (Start)
The sequence of terms together with their prime signatures begins:
                        4: (2)
                       24: (3,1)
                      216: (3,3)
                     2160: (4,3,1)
                    30240: (5,3,1,1)
                   453600: (5,4,2,1)
                  9525600: (5,5,2,2)
                209563200: (6,5,2,2,1)
               5239080000: (6,5,4,2,1)
             136216080000: (7,5,4,2,1,1)
            4495130640000: (7,6,4,2,2,1)
          152834441760000: (8,6,4,2,2,1,1)
         5349205461600000: (8,6,5,3,2,1,1)
       203269807540800000: (9,6,5,3,2,1,1,1)
      7927522494091200000: (9,7,5,3,2,2,1,1)
    364666034728195200000: (10,7,5,3,2,2,1,1,1)
  17868635701681564800000: (10,7,5,5,2,2,1,1,1)
(End)

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

Partial sums of semiprimes are A062198.
First differences of semiprimes are A065516.
A000040 lists primes, with partial products A002110 (primorials).
A000142 lists factorials, with partial products A000178 (superfactorials).
A001358 lists semiprimes, with partial products A112141 (this sequence).
A005117 lists squarefree numbers, with partial products A111059.
A006881 lists squarefree semiprimes, with partial products A339191.
A101048 counts partitions into semiprimes (restricted: A338902).
A320655 counts factorizations into semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes, with product/sum/difference A087794/A176504/A176506.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.
Cf. A001222, A084126, A084127, A115392, A168472, A320732, A320892.

A112141 := proc(n)
    mul(A001358(i),i=1..n) ;
end proc:
seq(A112141(n),n=1..10) ; # R. J. Mathar, Jun 30 2020

NextSemiPrime[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sp = n + sgn; While[c < Abs[k], While[ PrimeOmega[sp] != 2, If[sgn < 0, sp--, sp++]]; If[sgn < 0, sp--, sp++]; c++]; sp + If[sgn < 0, 1, -1]]; f[n_] := Times @@ NestList[ NextSemiPrime@# &, 2^2, n - 1]; Array[f, 18] (* Robert G. Wilson v, Jun 13 2013 *)
FoldList[Times,Select[Range[30],PrimeOmega[#]==2&]] (* Gus Wiseman, Dec 06 2020 *)

a(n)=my(v=vector(n),i,k=3);while(iCharles R Greathouse IV, Apr 04 2013

from sympy import factorint
def aupton(terms):
    alst, k, p = [], 1, 1
    while len(alst) < terms:
        if sum(factorint(k).values()) == 2:
            p *= k
            alst.append(p)
        k += 1
    return alst
print(aupton(18)) # Michael S. Branicky, Aug 31 2021

a(n) = Product_{i=1..n} A001358(i).

A001222(a(n)) = 2*n.

A322608 Values of k such that (product of squarefree numbers <= k) / (sum of squarefree numbers <= k) is an integer.

Original entry on oeis.org

1, 3, 11, 14, 17, 21, 23, 33, 34, 37, 46, 47, 55, 58, 59, 61, 62, 67, 69, 73, 82, 83, 87, 94, 95, 97, 101, 106, 107, 109, 114, 115, 119, 123, 127, 133, 134, 141, 146, 151, 157, 158, 159, 161, 165, 166, 173, 181, 187, 197, 202, 203, 210, 218, 219, 223, 226, 230

Offset: 1

Paolo P. Lava, Dec 20 2018

			3 is in the sequence because (1*2*3)/(1+2+3) = 1.
11 is in the sequence because (1*2*3*5*6*7*10*11)/(1+2+3+5+6+7+10+11) = 138600/45 = 3080.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A005117, A051838, A111059, A116536, A141092, A173143, A322607 (values of the quotient), A347690 (numbers of terms in the numerators).

with(numtheory): P:=proc(q) local a,b,c,n; a:=1; b:=0; c:=[];
for n from 1 to q do if issqrfree(n) then a:=a*n; b:=b+n;
if frac(a/b)=0 then c:=[op(c),n];
fi; fi; od; op(c); end: P(60);

seq = {}; sum = 0; prod = 1; Do[If[SquareFreeQ[n], sum += n; prod *= n; If[Divisible[prod, sum], AppendTo[seq, n]]], {n, 1, 230}]; seq (* Amiram Eldar, Mar 05 2021 *)

Definition corrected by N. J. A. Sloane, Sep 19 2021 at the suggestion of Harvey P. Dale.

A322607 Numbers that can be expressed as the ratio between the product and the sum of consecutive squarefree numbers starting from 1.

Original entry on oeis.org

1, 3080, 350350, 61850250, 17823180375, 6871260396000, 88909822914869880000, 2746644314348614680000, 2980109081068246927800000, 9638057975990853416623724908800000, 424217819372970387341691005411520000, 51912228216508515627667235880347808000000, 152157812632066726080765311397008321568000000

Offset: 1

Paolo P. Lava, Dec 20 2018

			1 is a term because 1/1 = (1*2*3)/(1+2+3) = 1.
3080 is a term because (1*2*3*5*6*7*10*11)/(1+2+3+5+6+7+10+11) = 138600/45 = 3080.

Cf. A005117, A111059, A116536, A141092, A173143, A322608.

with(numtheory): P:=proc(q) local a,b,c,n; a:=1; b:=0; c:=[];
for n from 1 to q do if issqrfree(n) then a:=a*n; b:=b+n;
if frac(a/b)=0 then if n>1 then c:=[op(c),a/b];
fi; fi; fi; od; op(c); end: P(60);

a(n) = A111059(A322608(n+1))/A173143(A322608(n+1)).

A339191 Partial products of squarefree semiprimes (A006881).

Original entry on oeis.org

6, 60, 840, 12600, 264600, 5821200, 151351200, 4994589600, 169816046400, 5943561624000, 225855341712000, 8808358326768000, 405184483031328000, 20664408634597728000, 1136542474902875040000, 64782921069463877280000, 3757409422028904882240000

Offset: 1

Gus Wiseman, Nov 30 2020

nonn

A squarefree semiprime is a product of any two distinct prime numbers.

Do all terms belong to A242031 (weakly decreasing prime signature)?

			The sequence of terms together with their prime indices begins:
          6: {1,2}
         60: {1,1,2,3}
        840: {1,1,1,2,3,4}
      12600: {1,1,1,2,2,3,3,4}
     264600: {1,1,1,2,2,2,3,3,4,4}
    5821200: {1,1,1,1,2,2,2,3,3,4,4,5}
  151351200: {1,1,1,1,1,2,2,2,3,3,4,4,5,6}
The sequence of terms together with their prime signatures begins:
                   6: (1,1)
                  60: (2,1,1)
                 840: (3,1,1,1)
               12600: (3,2,2,1)
              264600: (3,3,2,2)
             5821200: (4,3,2,2,1)
           151351200: (5,3,2,2,1,1)
          4994589600: (5,4,2,2,2,1)
        169816046400: (6,4,2,2,2,1,1)
       5943561624000: (6,4,3,3,2,1,1)
     225855341712000: (7,4,3,3,2,1,1,1)
    8808358326768000: (7,5,3,3,2,2,1,1)
  405184483031328000: (8,5,3,3,2,2,1,1,1)

A000040 lists the primes, with partial products A002110 (primorials).
A001358 lists semiprimes, with partial products A112141.
A002100 counts partitions into squarefree semiprimes (restricted: A338903)
A000142 lists factorial numbers, with partial products A000178.
A005117 lists squarefree numbers, with partial products A111059.
A006881 lists squarefree semiprimes, with partial sums A168472.
A166237 gives first differences of squarefree semiprimes.
A320655 counts factorizations into semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A338898/A338912/A338913 give prime indices of semiprimes.
A338899/A270650/A270652 give prime indices of squarefree semiprimes.
A338901 gives first appearances in the list of squarefree semiprimes.
A339113 gives products of primes of squarefree semiprime index.
Cf. A001221, A112798, A167171, A320732, A320891, A320892, A320894, A320911, A338900, A338902, A339003, A339004.

FoldList[Times,Select[Range[20],SquareFreeQ[#]&&PrimeOmega[#]==2&]]

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

A209862 Permutation of nonnegative integers which maps A209642 into ascending order (A209641).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A112141 Product of the first n semiprimes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A322608 Values of k such that (product of squarefree numbers <= k) / (sum of squarefree numbers <= k) is an integer.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

A322607 Numbers that can be expressed as the ratio between the product and the sum of consecutive squarefree numbers starting from 1.

Views

Author

Keywords

Examples

Crossrefs

Programs

Maple

Formula

A339191 Partial products of squarefree semiprimes (A006881).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A347690 Values of k such that (product of first k squarefree numbers) / (sum of first k squarefree numbers) is an integer.

Views

Author

Keywords

Comments

Links

Crossrefs