A057588 -id:A057588 - cp's OEIS frontend

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

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

Offset: 0

N. J. A. Sloane and J. H. Conway

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

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

A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 50.
G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, Cambridge, University Press, 1940, p. 49.
P. Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 4.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 114.
D. Wolfe and S. Hirshberg, Underspecified puzzles, in Tribute to A Mathemagician, Peters, 2005, pp. 73-74.

Alex Ermolaev, Table of n, a(n) for n = 0..350 (terms up to a(100) from T. D. Noe)
Iskander Aliev, Jesús De Loera, Fritz Eisenbrand, Timm Oertel, and Robert Weismantel, The Support of Integer Optimal Solutions, arXiv:1712.08923 [math.OC], 2017.
C. K. Caldwell, The Prime Glossary, Primorial.
Geoffrey Caveney, J.-L. Nicolas and J. Sondow, On SA, CA, and GA numbers, arXiv:1112.6010 [math.NT], 2011-2012; Ramanujan J., 29 (2012), 359-384.
Harvey Dubner, Factorial and primorial primes, J. Rec. Math., Vol. 19, No. 3 (1987), pp. 197-203. (Annotated scanned copy)
F. Ellermann, Illustration for A002110, A005867, A038110, A060753.
S. W. Golomb, The evidence for Fortune's conjecture, Math. Mag. 54 (1981), 209-210.
D. J. Greenhoe, MRA-Wavelet subspace architecture for logic, probability, and symbolic sequence processing, 2014.
Daniel J. Greenhoe, Frames and Bases: Structure and Design, Version 0.20, Signal Processing ABCs series (2019) Vol. 4, pp. 7, 81.
Daniel J. Greenhoe, A Book Concerning Transforms, Version 0.10, Signal Processing ABCs series (2019) Vol. 5, see page 7.
A. W. Lin and S. Zhou, A linear-time algorithm for the orbit problem over cyclic groups, preprint, CONCUR 2014 - Concurrency Theory, Volume 8704 of the series Lecture Notes in Computer Science pp. 327-341.
A. W. Lin and S. Zhou, A linear-time algorithm for the orbit problem over cyclic groups, CONCUR 2014 - Concurrency Theory, Lecture Notes in Computer Science, Volume 8704, 2014, pp. 327-341.
F. E. Masat, Letter to N. J. A. Sloane with attachment: "A note on prime number sequences" (unpublished manuscript), Apr. 1991.
R. Mestrovic, Euclid's theorem on the infinitude of primes: a historical survey of its proofs (300 BC--2012) and another new proof, arXiv:1202.3670 [math.HO], 2012.
Thomas Morrill, Further Development of "Non-Pythagorean" Musical Scales Based on Logarithms, arXiv:1804.08067 [math.HO], 2018.
J.-L. Nicolas, Petites valeurs de la fonction d'Euler, J. Number Theory 17, no.3 (1983), 375-388.
Patrick Sole and Michel Planat, The Robin inequality for 7-free integers, INTEGERS, 2011, #A65.
Andrew V. Sutherland, Order Computations in Generic Groups, Ph. D. Dissertation, Math. Dept., M.I.T., 2007.
G. Villemin's Almanach of Numbers, Primorielle.
Eric Weisstein's World of Mathematics, Primorial.
Robert G. Wilson v, Letter to N. J. A. Sloane, Jan. 1994.
Index to divisibility sequences
Index entries for "core" sequences
Index entries for sequences related to primorial base
Index entries for sequences related to primorial numbers

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

a002110 n = product $ take n a000040_list
a002110_list = scanl (*) 1 a000040_list
-- Reinhard Zumkeller, Feb 19 2012, May 03 2011

[1] cat [&*[NthPrime(i): i in [1..n]]: n in [1..20]]; // Bruno Berselli, Oct 24 2012

[1] cat [&*PrimesUpTo(p): p in PrimesUpTo(60)]; // Bruno Berselli, Feb 08 2015

A002110 := n -> mul(ithprime(i),i=1..n);

FoldList[Times, 1, Prime[Range[20]]]
primorial[n_] := Product[Prime[i], {i, n}]; Array[primorial,20] (* José María Grau Ribas, Feb 15 2010 *)
Join[{1}, Denominator[Accumulate[1/Prime[Range[20]]]]] (* Harvey P. Dale, Apr 11 2012 *)

a(n)=prod(i=1,n, prime(i)) \\ Washington Bomfim, Sep 23 2008

p=1; for (n=0, 100, if (n, p*=prime(n)); write("b002110.txt", n, " ", p) )  \\ Harry J. Smith, Nov 13 2009

a(n) = factorback(primes(n)) \\ David A. Corneth, May 06 2018

from sympy import primorial
def a(n): return 1 if n < 1 else primorial(n)
[a(n) for n in range(51)]  # Indranil Ghosh, Mar 29 2017

[sloane.A002110(n) for n in (1..20)] # Giuseppe Coppoletta, Dec 05 2014

; with memoization-macro definec
(definec (A002110 n) (if (zero? n) 1 (* (A000040 n) (A002110 (- n 1))))) ;; Antti Karttunen, Aug 30 2016

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

A276086 Primorial base exp-function: digits in primorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Original entry on oeis.org

1, 2, 3, 6, 9, 18, 5, 10, 15, 30, 45, 90, 25, 50, 75, 150, 225, 450, 125, 250, 375, 750, 1125, 2250, 625, 1250, 1875, 3750, 5625, 11250, 7, 14, 21, 42, 63, 126, 35, 70, 105, 210, 315, 630, 175, 350, 525, 1050, 1575, 3150, 875, 1750, 2625, 5250, 7875, 15750, 4375, 8750, 13125, 26250, 39375, 78750, 49, 98, 147, 294, 441, 882, 245, 490, 735, 1470, 2205, 4410, 1225, 2450

Offset: 0

Antti Karttunen, Aug 21 2016

Prime product form of primorial base expansion of n.
Sequence is a permutation of A048103. It maps the smallest prime not dividing n to the smallest prime dividing n, that is, A020639(a(n)) = A053669(n) holds for all n >= 1.
The sequence satisfies the exponential function identity, a(x + y) = a(x) * a(y), whenever A329041(x,y) = 1, that is, when adding x and y together will not generate any carries in the primorial base. Examples of such pairs of x and y are A328841(n) & A328842(n), and also A328770(n) (when added with itself). - Antti Karttunen, Oct 31 2019
From Antti Karttunen, Feb 18 2022: (Start)
The conjecture given in A327969 asks whether applying this function together with the arithmetic derivative (A003415) in some combination or another can eventually transform every positive integer into zero.
Another related open question asks whether there are any other numbers than n=6 such that when starting from that n and by iterating with A003415, one eventually reaches a(n). See comments in A351088.
This sequence is used in A351255 to list the terms of A099308 in a different order, by the increasing exponents of the successive primes in their prime factorization. (End)
From Bill McEachen, Oct 15 2022: (Start)
From inspection, the least significant decimal digits of a(n) terms form continuous chains of 30 as follows. For n == i (mod 30), i=0..5, there are 6 ordered elements of these 8 {1,2,3,6,9,8,7,4}. Then for n == i (mod 30), i=6..29, there are 12 repeated pairs = {5,0}.
Moreover, when the individual elements of any of the possible groups of 6 are transformed via (7*digit) (mod 10), the result matches one of the other 7 groupings (not all 7 may be seen). As example, {1,2,3,6,9,8} transforms to {7,4,1,2,3,6}. (End)
The least significant digit of a(n) in base 4 is given by A353486, and in base 6 by A358840. - Antti Karttunen, Oct 25 2022, Feb 17 2024

			For n = 24, which has primorial base representation (see A049345) "400" as 24 = 4*A002110(2) + 0*A002110(1) + 0*A002110(0) = 4*6 + 0*2 + 0*1, thus a(24) = prime(3)^4 * prime(2)^0 * prime(1)^0 = 5^4 = 625.
For n = 35 = "1021" as 35 = 1*A002110(3) + 0*A002110(2) + 2*A002110(1) + 1*A002110(0) = 1*30 + 0*6 + 2*2 + 1*1, thus a(35) = prime(4)^1 * prime(2)^2 * prime(1) = 7 * 3*3 * 2 = 126.

Antti Karttunen, Table of n, a(n) for n = 0..2310
Antti Karttunen, Program in LODA-assembly
Antti Karttunen, Program in LODA-assembly [Cached copy]
Index entries for sequences related to primorial base

Cf. A276085 (a left inverse) and also A276087, A328403.
Cf. A000040, A001221, A001222, A002110, A020639, A049345, A053669, A055396, A057588, A071178, A143293, A257993, A267263, A276084, A276088, A276092, A276093, A276147, A276150, A276151, A276153, A276156, A283477, A324198 (= gcd(n, a(n))), A328584 (= lcm(n, a(n))), A324646, A324289, A328386, A328403, A328475, A328571, A328572, A328578, A328612, A328613, A328620, A328624, A328627, A328763, A328766, A328828, A328835, A328841, A328842, A328843, A328844, A329041, A324580 [= n*a(n)], A324895 (largest proper divisor of a(n)), A351252, A353486 (reduced modulo 4), A358840 (modulo 6), A353489, A353516.
Cf. A048103 (terms sorted into ascending order), A100716 (natural numbers not present in this sequence).
Cf. A278226 (associated filter-sequence), A286626 (and its rgs-version), A328477.
Cf. A328316 (iterates started from zero).
Cf. A327858, A327859, A327860, A327963, A328097, A328098, A328099, A328110, A328112, A328382 for various combinations with arithmetic derivative (A003415).
Cf. also A327167, A329037.
Cf. A019565 and A054842 for base-2 and base-10 analogs and A276076 for the analogous "factorial base exp-function", from which this differs for the first time at n=24, where a(24)=625 while A276076(24)=7.
Cf. A327969, A351088, A351458 for sequences with conjectures involving this sequence.

b = MixedRadix[Reverse@ Prime@ Range@ 12]; Table[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[n, b], {n, 0, 51}] (* Michael De Vlieger, Aug 23 2016, Version 10.2 *)
f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; Table[Times @@ Flatten@ MapIndexed[Prime[#2]^#1 &, Reverse@ f@ n], {n, 0, 73}] (* Michael De Vlieger, Aug 30 2016, Pre-Version 10 *)
a[n0_] := Module[{m = 1, i = 1, n = n0, p}, While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Dec 01 2021, after Antti Karttunen's Sage code *)

A276086(n) = { my(i=0,m=1,pr=1,nextpr); while((n>0),i=i+1; nextpr = prime(i)*pr; if((n%nextpr),m*=(prime(i)^((n%nextpr)/pr));n-=(n%nextpr));pr=nextpr); m; }; \\ Antti Karttunen, May 12 2017

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; \\ (Better than above one, avoids unnecessary construction of primorials). - Antti Karttunen, Oct 14 2019

from sympy import prime
def a(n):
    i=0
    m=pr=1
    while n>0:
        i+=1
        N=prime(i)*pr
        if n%N!=0:
            m*=(prime(i)**((n%N)/pr))
            n-=n%N
        pr=N
    return m # Indranil Ghosh, May 12 2017, after Antti Karttunen's PARI code

from sympy import nextprime
def a(n):
    m, p = 1, 2
    while n > 0:
        n, r = divmod(n, p)
        m *= p**r
        p = nextprime(p)
    return m
print([a(n) for n in range(74)])  # Peter Luschny, Apr 20 2024

def A276086(n):
    m=1
    i=1
    while n>0:
        p = sloane.A000040(i)
        m *= (p**(n%p))
        n = floor(n/p)
        i += 1
    return (m)
# Antti Karttunen, Oct 14 2019, after Indranil Ghosh's Python code above, and my own leaner PARI code from Oct 14 2019. This avoids unnecessary construction of primorials.

(define (A276086 n) (let loop ((n n) (t 1) (i 1)) (if (zero? n) t (let* ((p (A000040 i)) (d (modulo n p))) (loop (/ (- n d) p) (* t (expt p d)) (+ 1 i))))))

(definec (A276086 n) (if (zero? n) 1 (* (expt (A053669 n) (A276088 n)) (A276086 (A276093 n))))) ;; Needs macro definec from http://oeis.org/wiki/Memoization#Scheme

(definec (A276086 n) (if (zero? n) 1 (* (A053669 n) (A276086 (- n (A002110 (A276084 n))))))) ;; Needs macro definec from http://oeis.org/wiki/Memoization#Scheme

a(0) = 1; for n >= 1, a(n) = A053669(n) * a(A276151(n)) = A053669(n) * a(n-A002110(A276084(n))).
a(0) = 1; for n >= 1, a(n) = A053669(n)^A276088(n) * a(A276093(n)).
a(n) = A328841(a(n)) + A328842(a(n)) = A328843(n) + A328844(n).
a(n) = a(A328841(n)) * a(A328842(n)) = A328571(n) * A328572(n).
a(n) = A328475(n) * A328580(n) = A328476(n) + A328580(n).
a(A002110(n)) = A000040(n+1). [Maps primorials to primes]
a(A143293(n)) = A002110(n+1). [Maps partial sums of primorials to primorials]
a(A057588(n)) = A276092(n).
a(A276156(n)) = A019565(n).
a(A283477(n)) = A324289(n).
a(A003415(n)) = A327859(n).
Here the text in brackets shows how the right hand side sequence is a function of the primorial base expansion of n:
A001221(a(n)) = A267263(n). [Number of nonzero digits]
A001222(a(n)) = A276150(n). [Sum of digits]
A067029(a(n)) = A276088(n). [The least significant nonzero digit]
A071178(a(n)) = A276153(n). [The most significant digit]
A061395(a(n)) = A235224(n). [Number of significant digits]
A051903(a(n)) = A328114(n). [Largest digit]
A055396(a(n)) = A257993(n). [Number of trailing zeros + 1]
A257993(a(n)) = A328570(n). [Index of the least significant zero digit]
A079067(a(n)) = A328620(n). [Number of nonleading zeros]
A056169(a(n)) = A328614(n). [Number of 1-digits]
A056170(a(n)) = A328615(n). [Number of digits larger than 1]
A277885(a(n)) = A328828(n). [Index of the least significant digit > 1]
A134193(a(n)) = A329028(n). [The least missing nonzero digit]
A005361(a(n)) = A328581(n). [Product of nonzero digits]
A072411(a(n)) = A328582(n). [LCM of nonzero digits]
A001055(a(n)) = A317836(n). [Number of carry-free partitions of n in primorial base]
Various number theoretical functions applied:
A000005(a(n)) = A324655(n). [Number of divisors of a(n)]
A000203(a(n)) = A324653(n). [Sum of divisors of a(n)]
A000010(a(n)) = A324650(n). [Euler phi applied to a(n)]
A023900(a(n)) = A328583(n). [Dirichlet inverse of Euler phi applied to a(n)]
A069359(a(n)) = A329029(n). [Sum a(n)/p over primes p dividing a(n)]
A003415(a(n)) = A327860(n). [Arithmetic derivative of a(n)]
Other identities:
A276085(a(n)) = n.          [A276085 is a left inverse]
A020639(a(n)) = A053669(n). [The smallest prime not dividing n -> the smallest prime dividing n]
A046523(a(n)) = A278226(n). [Least number with the same prime signature as a(n)]
A246277(a(n)) = A329038(n).
A181819(a(n)) = A328835(n).
A053669(a(n)) = A326810(n), A326810(a(n)) = A328579(n).
A257993(a(n)) = A328570(n), A328570(a(n)) = A328578(n).
A328613(a(n)) = A328763(n), A328620(a(n)) = A328766(n).
A328828(a(n)) = A328829(n).
A053589(a(n)) = A328580(n). [Greatest primorial number which divides a(n)]
A276151(a(n)) = A328476(n). [... and that primorial subtracted from a(n)]
A111701(a(n)) = A328475(n).
A328114(a(n)) = A328389(n). [Greatest digit of primorial base expansion of a(n)]
A328389(a(n)) = A328394(n), A328394(a(n)) = A328398(n).
A235224(a(n)) = A328404(n), A328405(a(n)) = A328406(n).
a(A328625(n)) = A328624(n), a(A328626(n)) = A328627(n). ["Twisted" variants]
a(A108951(n)) = A324886(n).
a(n) mod n = A328386(n).
a(a(n)) = A276087(n), a(a(a(n))) = A328403(n). [2- and 3-fold applications]
a(2n+1) = 2 * a(2n). - Antti Karttunen, Feb 17 2022

Name edited and new link-formulas added by Antti Karttunen, Oct 29 2019

Name changed again by Antti Karttunen, Feb 05 2022

A048103 Numbers not divisible by p^p for any prime p.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98

Offset: 1

N. J. A. Sloane

If a(n) = Product p_i^e_i then p_i > e_i for all i.
Complement of A100716; A129251(a(n)) = 0. - Reinhard Zumkeller, Apr 07 2007
Density is 0.72199023441955... = Product_{p>=2} (1 - p^-p) where p runs over the primes. - Charles R Greathouse IV, Jan 25 2012
A027748(a(n),k) <= A124010(a(n),k), 1<=k<=A001221(a(n)). - Reinhard Zumkeller, Apr 28 2012
Range of A276086. Also numbers not divisible by m^m for any natural number m > 1. - Antti Karttunen, Nov 18 2024

			6 = 2^1 * 3^1 is OK but 12 = 2^2 * 3^1 is not.
625 = 5^4 is present because it is not divisible by 5^5.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Complement: A100716.
Positions of 0's in A129251, A342023, A376418, positions of 1's in A327936, A342007, A359550 (characteristic function).
Cf. A048102, A048104, A051674 (p^p), A054743, A054744, A377982 (a left inverse, partial sums of char. fun, see also A328402).
Cf. A276086 (permutation of this sequence, see also A376411, A376413).
Subsequences: A002110, A005117, A006862, A024451 (after its initial 0), A057588, A099308 (after its initial 0), A276092, A328387, A328832, A359547, A370114, A371083, A373848, A377871, A377992.
Disjoint union of {1}, A327934 and A358215.
Also A276078 is a subsequence, from which this differs for the first time at n=451 where a(451)=625, while that value is missing from A276078.

a048103 n = a048103_list !! (n-1)
a048103_list = filter (\x -> and $
   zipWith (>) (a027748_row x) (map toInteger $ a124010_row x)) [1..]
-- Reinhard Zumkeller, Apr 28 2012

{1}~Join~Select[Range@ 120, Times @@ Boole@ Map[First@ # > Last@ # &, FactorInteger@ #] > 0 &] (* Michael De Vlieger, Aug 19 2016 *)

isok(n) = my(f=factor(n)); for (i=1, #f~, if (f[i,1] <= f[i,2], return(0))); return(1); \\ Michel Marcus, Nov 13 2020

A359550(n) = { my(pp); forprime(p=2, , pp = p^p; if(!(n%pp), return(0)); if(pp > n, return(1))); }; \\ (A359550 is the characteristic function for A048103) - Antti Karttunen, Nov 18 2024

from itertools import count, islice
from sympy import factorint
def A048103_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:all(map(lambda d:d[1]A048103_list = list(islice(A048103_gen(),30)) # Chai Wah Wu, Jan 05 2023

;; With Antti Karttunen's IntSeq-library.
(define A048103 (ZERO-POS 1 1 A129251))
;; Antti Karttunen, Aug 18 2016

a(n) ~ kn with k = 1/Product_{p>=2}(1 - p^-p) = Product_{p>=2}(1 + 1/(p^p - 1)) = 1.3850602852..., where the product is over all primes p. - Charles R Greathouse IV, Jan 25 2012

For n >= 1, A377982(a(n)) = n. - Antti Karttunen, Nov 18 2024

More terms from James Sellers, Apr 22 2000

A006862 Euclid numbers: 1 + product of the first n primes.

Original entry on oeis.org

2, 3, 7, 31, 211, 2311, 30031, 510511, 9699691, 223092871, 6469693231, 200560490131, 7420738134811, 304250263527211, 13082761331670031, 614889782588491411, 32589158477190044731, 1922760350154212639071, 117288381359406970983271, 7858321551080267055879091

Offset: 0

Simon Plouffe and N. J. A. Sloane

It is an open question whether all terms of this sequence are squarefree.
a(n) is the smallest x > 1 such that x^prime(n) == 1 (mod prime(i)) i=1,2,3,...,n-1. - Benoit Cloitre, May 30 2002
Numbers n such that n/phi(n-1) is a record. - Arkadiusz Wesolowski, Nov 22 2012
Nyblom (theorem 2.3) proves that this sequence contains no proper powers, e.g., is a subsequence of A007916. - Charles R Greathouse IV, Mar 02 2016
It is an open question if there are an infinite number of prime Euclid numbers. - Mike Winkler, Feb 05 2017
These numbers are not pairwise relatively prime; the first example is gcd(a(7), a(17)) = 277. Also gcd(a(47), a(131)) = 1051, which is probably the second example (wrt. greater index which is here 131). It is easy to find other primes like 277 and 1051. - Jeppe Stig Nielsen, Mar 24 2017
Subsequence of A048103. Proof: For all primes p, when i >= A000720(p), neither p itself nor p^p divides a(i), but neither does p^p divide a(i) when i < A000720(p), as p^p > 1 + A034386(p). - Antti Karttunen, Nov 17 2024

			It is a universal convention that an empty product is 1 (just as an empty sum is 0), and since this sequence has offset 0, the first term is 1+1 = 2. - _N. J. A. Sloane_, Dec 02 2015

John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 134.
J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 211, p. 61, Ellipses, Paris 2008.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
F. Smarandache, Properties of numbers, Arizona State University Special Collections, 1973.
I. Vardi, Computational Recreations in Mathematica, Addison-Wesley, 1991, sections 5.1 and 5.2.
S. Wagon, Mathematica in Action, Freeman, NY, 1991, p. 35.

Derek Maciel, Table of n, a(n) for n = 0..348 (first 101 terms from T. D. Noe)
S. W. Golomb, The evidence for Fortune's conjecture, Math. Mag. 54 (1981), 209-210.
H. Ibstedt, A Few Smarandache Sequences, Smarandache Notions Journal, Vol. 8, No. 1-2-3, 1997, 170-183.
R. Mestrovic, Euclid's theorem on the infinitude of primes: a historical survey of its proofs (300 BC--2012) and another new proof, arXiv preprint arXiv:1202.3670 [math.HO], 2012.
Hisanori Mishima, Factorizations of many number sequences
Michael A. Nyblom, On the construction of a family of almost power free sequences, Fibonacci Quart. 46/47 (2008/2009), no. 4, 366-368.
Shubhankar Paul, Ten Problems of Number Theory, International Journal of Engineering and Technical Research (IJETR), ISSN: 2321-0869, Volume-1, Issue-9, November 2013.
Eric Weisstein's World of Mathematics, Euclid Number
Eric Weisstein's World of Mathematics, Fortunate Prime
R. G. Wilson v, Explicit factorizations

Cf. A005867, A007916, A014545, A018239 (primes in sequence), A034386, A057588, A377871.

Subsequence of A048103.

[2] cat [&*PrimesUpTo(p)+1: p in PrimesUpTo(70)]; // Vincenzo Librandi, Dec 03 2015

with(numtheory): A006862 := proc(n) local i; if n=0 then 2 else 1+product('ithprime(i)','i'=1..n); fi; end;
# second Maple program:
a:= proc(n) option remember; `if`(n=0, 2,
      1+ithprime(n)*(a(n-1)-1))
    end:
seq(a(n), n=0..20);  # Alois P. Heinz, Feb 06 2021

Table[Product[Prime[k], {k, 1, n}] + 1, {n, 1, 18}]
1 + FoldList[Times, 1, Prime@ Range@ 19] (* Harvey P. Dale, Dec 02 2015 and modified by Robert G. Wilson v, Mar 25 2017 *)

a(n)=my(v=primes(n)); prod(i=1,#v,v[i])+1 \\ Charles R Greathouse IV, Nov 20 2012

from sympy import primorial
def A006862(n):
    if n == 0: return 2
    else: return 1 + primorial(n) # Karl-Heinz Hofmann, Aug 21 2024

a(n) = A002110(n) + 1.

For n >= 1, a(n) = A057588(n) + 2. - Antti Karttunen, Nov 17 2024

A369054 Number of representations of n as a sum (pq + pr + q*r) with three odd primes p <= q <= r.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1

Offset: 0

Antti Karttunen, Jan 20 2024

nonn

Number of solutions to n = x', where x' is the arithmetic derivative of x (A003415), and x is a product of three odd primes (not all necessarily distinct, A046316).

See the conjecture in A369055.

			a(27) = 1 as 27 can be expressed in exactly one way in the form (p*q + p*r + q*r), with p, q, r all being 3 in this case, as 27 = (3*3 + 3*3 + 3*3).
a(311) = 5 as 311 = (3*5 + 3*37 + 5*37) = (3*7 + 3*29 + 7*29) = (3*13 + 3*17 + 13*17) = (5*7 + 5*23 + 7*23) = (7*11 + 7*13 + 11*13). Expressed in the terms of arithmetic derivatives, of the A099302(311) = 8 antiderivatives of 311 [366, 430, 494, 555, 609, 663, 805, 1001], only the last five are products of three odd primes: 555 = 3*5*37, 609 = 3*7*29, 663 = 3*13*17, 805 = 5*7*23, 1001 = 7 * 11 * 13.

Antti Karttunen, Table of n, a(n) for n = 0..50000
Victor Ufnarovski and Bo Ahlander, How to Differentiate a Number, J. Integer Seqs., Vol. 6, 2003.

Cf. A369055 [quadrisection, a(4n-1)], and its trisections A369460 [= a((12*n)-9)], A369461 [= a((12*n)-5)], A369462 [= a((12*n)-1)].
Cf. A002620, A003415, A046316, A099302, A369056, A369058, A369252.
Cf. A369251 (positions of terms > 0), A369464 (positions of 0's).
Cf. A369063 (positions of records), A369064 (values of records).
Cf. A369241 [= a(2^n - 1)], A369242 [= a(n!-1)], A369245 [= a(A006862(n))], A369247 [= a(3*A057588(n))].

\\ Use this for building up a list up to a certain n. We iterate over weakly increasing triplets of odd primes:
A369054list(up_to) = { my(v = [3,3,3], ip = #v, d, u = vector(up_to)); while(1, d = ((v[1]*v[2]) + (v[1]*v[3]) + (v[2]*v[3])); if(d > up_to, ip--, ip = #v; u[d]++); if(!ip, return(u)); v[ip] = nextprime(1+v[ip]); for(i=1+ip,#v,v[i]=v[i-1])); };
v369054 = A369054list(100001);
A369054(n) = if(!n,n,v369054[n]);

\\ Use this for computing the value of arbitrary n. We iterate over weakly increasing pairs of odd primes:
A369054(n) = if(3!=(n%4),0, my(v = [3,3], ip = #v, r, c=0); while(1, r = (n-(v[1]*v[2])) / (v[1]+v[2]); if(r < v[2], ip--, ip = #v; if(1==denominator(r) && isprime(r),c++)); if(!ip, return(c)); v[ip] = nextprime(1+v[ip]); for(i=1+ip,#v,v[i]=v[i-1])));

a(n) = Sum_{i=1..A002620(n)} A369058(i)*[A003415(i)==n], where [ ] is the Iverson bracket.

For n >= 2, a(n) <= A099302(n).

A057705 Primorial primes: primes p such that p+1 is a primorial number (A002110).

Original entry on oeis.org

5, 29, 2309, 30029, 304250263527209, 23768741896345550770650537601358309, 19361386640700823163471425054312320082662897612571563761906962414215012369856637179096947335243680669607531475629148240284399976569

Offset: 1

Labos Elemer, Oct 24 2000

Chris K. Caldwell's The Top Twenty, Primorial.
Romeo Mestrovic, Euclid's theorem on the infinitude of primes: a historical survey of its proofs (300 BC--2012) and another new proof, arXiv preprint arXiv:1202.3670 [math.HO], 2012. - _N. J. A. Sloane_, Jun 13 2012

See A006794 and A057704 (the main entries for this sequence) for more terms.
Cf. A014545, A018239.
Cf. A002110, A010051.
Subsequence of A057588.

a057705 n = a057705_list !! (n-1)
a057705_list = filter ((== 1) . a010051) a057588_list
-- Reinhard Zumkeller, Mar 27 2013

Select[FoldList[Times, 1, Prime[Range[70]]], PrimeQ[# - 1] &] - 1 (* Harvey P. Dale, Jan 27 2014 *)

a(n) = A002110(A057704(n)) - 1.

A100778 Integer powers of primorial numbers.

Original entry on oeis.org

1, 2, 4, 6, 8, 16, 30, 32, 36, 64, 128, 210, 216, 256, 512, 900, 1024, 1296, 2048, 2310, 4096, 7776, 8192, 16384, 27000, 30030, 32768, 44100, 46656, 65536, 131072, 262144, 279936, 510510, 524288, 810000, 1048576, 1679616, 2097152, 4194304, 5336100

Offset: 1

Amarnath Murthy, Nov 28 2004

Smallest squarefree numbers or their powers with distinct prime signatures. Or least numbers with prime signatures (p*q*r*...)^k, where p,q,r,... are primes and k is a whole number.
Also Heinz numbers of uniform integer partitions whose union is an initial interval of positive integers. An integer partition is uniform if all parts appear with the same multiplicity. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k). The sequence of all uniform integer partitions whose Heinz numbers belong to the sequence begins: (1), (11), (12), (111), (1111), (123), (11111), (1122), (111111), (1111111), (1234), (111222), (11111111), (111111111), (112233), (1111111111). - Gus Wiseman, Dec 26 2018
From Amiram Eldar, Sep 26 2023: (Start)
Intersection of A025487 and A072774.
The distinct terms of A046523(A072774(n)) in ascending orders.
The k-th power of the n-th primorial number, A002110(n)^k, has (k+1)^n divisors which are the set of the (k+1)-free prime(n)-smooth numbers. (End)

			10 is not a term as 6 is a member with the same prime signature 10 > 6.
216 is a term as 216 = (2*3)^3. 243 is not a term as 32 represents that prime signature.

David A. Corneth, Table of n, a(n) for n = 1..8606 (terms <= 10^1000)

Cf. A000961, A001597, A002110, A007947, A025487, A046523, A055932, A056239, A057588, A072774, A072777, A112798, A304250, A319169, A322792, A322793.

unintQ[n_]:=And[SameQ@@Last/@FactorInteger[n],Length[FactorInteger[n]]==PrimePi[FactorInteger[n][[-1,1]]]];
Select[Range[1000],unintQ] (* Gus Wiseman, Dec 26 2018 *)

Sum_{n>=1} 1/a(n) = 1 + Sum_{n>=1} 1/A057588(n) = 2.2397359032... - Amiram Eldar, Oct 20 2020; corrected by Hal M. Switkay and Amiram Eldar, Apr 12 2021

More terms and simpler definition from Ray Chandler, Nov 29 2004

cp's OEIS Frontend

A369247 Number of representations of 3*A057588(n) as a sum (p*q + p*r + q*r) with three odd primes p <= q <= r, where A057588 is Kummer numbers, one less than primorials.

Views

Author

Keywords

Crossrefs

Programs

PARI

Formula

A068489 m for which prime(m) is the least prime dividing #prime(n) - 1, i.e., one less than primorial n-th prime (A057588).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A125549 Composite Kummer numbers: -1 + product of first n consecutive primes (A057588).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Extensions

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Sage

Scheme

Formula

A276086 Primorial base exp-function: digits in primorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Python

Sage

Scheme

Scheme

Scheme

Formula

Extensions

A048103 Numbers not divisible by p^p for any prime p.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

A369247 Number of representations of 3A057588(n) as a sum (pq + pr + qr) with three odd primes p <= q <= r, where A057588 is Kummer numbers, one less than primorials.

A369054 Number of representations of n as a sum (pq + pr + q*r) with three odd primes p <= q <= r.