A116536 -id:A116536 - cp's OEIS frontend

A159578 Dividend associated with A116536.

2, 30, 9699690, 304250263527210, 267064515689275851355624017992790, 5766152219975951659023630035336134306565384015606066319856068810, 962947420735983927056946215901134429196419130606213075415963491270, 29819592777931214269172453467810429868925511217482600306406141434158090

Offset: 1

Enoch Haga, Apr 16 2009

a(2)-a(4) are mentioned by Alladi and Erdős (1977). They conjectured that this sequence is infinite. - Amiram Eldar, Nov 02 2020

			a(2) = 30 because 2*3*5 = 30, 2+3+5 = 10, and 30/10 = 3 in A116536.

Amiram Eldar, Table of n, a(n) for n = 1..80
Krishnaswami Alladi and Paul Erdős, On an additive arithmetic function, Pacific Journal of Mathematics, Vol. 71, No. 2 (1977), pp. 275-294, alternative link. See p. 290.

Intersection of A002110 and A036844.

Cf. A051838, A116536.

# First define t1, the sequence A051838.
t1:=[1,3,8,13,23,38,39,41,43,48,50,53,56,57,58,66,68,
70,73,77,84,90,94,98,126,128,134,140,143,145,149,
151,153,157,160,164,167,168,172,174,176,182,191,
194,196,200,210,212,215,217,218,219,222,225,228,
229];
p:=ithprime;
num:=n->mul(p(i),i=1..t1[n]);
s:=[num(i),i=1..11)];

seq = {}; sum = 0; prod = 1; p = 1; Do[p = NextPrime[p]; prod *= p; sum += p; If[Divisible[prod, sum], AppendTo[seq, prod]], {50}]; seq (* Amiram Eldar, Nov 02 2020 *)
Module[{nn=50,s,p},s=Accumulate[Prime[Range[nn]]];p=FoldList[Times,Prime[Range[ nn]]]; Select[Thread[{p,s}],Divisible[#[[1]],#[[2]]]&]][[All,1]] (* Harvey P. Dale, Jun 07 2022 *)

a(n) = A002110(A051838(n)). - Amiram Eldar, Nov 02 2020

Corrected by N. J. A. Sloane, Oct 02 2011 (all the terms were wrong).

A159580 Integral quotients occur consecutively this many times (in sequence associated with A051838 and A116536).

Original entry on oeis.org

2, 3, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 3, 2, 2, 2, 3, 2, 2, 3, 3, 2, 3, 2, 4, 4, 2, 2, 2, 2, 2, 2, 3, 2, 2, 3, 2, 2, 3, 3, 2, 2, 2, 2, 2, 3, 3, 2, 3, 2, 2, 3, 2, 3, 3, 2, 2, 2, 3, 2, 2, 3, 2, 2, 2, 5, 2, 2, 2, 2, 2, 2, 2, 4, 2

Offset: 1

Enoch Haga, Apr 16 2009

			a(1)=2 because the 1st subsequence (38,39) of consecutive integers has run length 2.
a(2)=3 because the 2nd subsequence is (56,57,58) which has run length 3.
a(3)=2 because the 3rd subsequence is (167,168) which has run length 2.

Cf. A051838, A116536, A140763, A159578, A159581.

Recomputed based on b-file of A051838. - R. J. Mathar, Aug 28 2018

A159581 First divisors at which integral quotients occur consecutively in A116536 (and associated with A159580).

Original entry on oeis.org

2747, 6601, 75130, 133386, 148827, 208385, 255445, 799846, 814006, 890299, 993730, 1037571, 1049698, 1382738, 1723170, 1869711, 2168747, 2530165, 2569285, 2615298, 2838135, 2963531, 3020151, 3434151, 3510885, 3627674, 3706710, 3941521, 3999326, 4434411, 4700809, 5887533, 6133305

Offset: 1

Enoch Haga, Apr 16 2009

			The first run of consecutive integers in A051838 is A051838(6)=38 and A051838(7)=39, therefore a(1) = A140763(6)= 2747.
The second run of consecutive integers in A051838 is A051838(13)=56, A051838(14)= 57, A051838(15)=58, therefore a(2) = A140763(13) = 6601.

Cf. A051838, A116536, A140763, A159578, A159580.

Recomputed by R. J. Mathar, Aug 28 2018

A159577 Erroneous version of A116536.

Original entry on oeis.org

3, 125970, 44618574, 1050274875, 835668708875, 87740042053973197350, 249875247243953981660010, 958593992854426686094086750, 28077829144670491508241956182896750

Offset: 1

dead

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

Original entry on oeis.org

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

A007504 Sum of the first n primes.

Original entry on oeis.org

0, 2, 5, 10, 17, 28, 41, 58, 77, 100, 129, 160, 197, 238, 281, 328, 381, 440, 501, 568, 639, 712, 791, 874, 963, 1060, 1161, 1264, 1371, 1480, 1593, 1720, 1851, 1988, 2127, 2276, 2427, 2584, 2747, 2914, 3087, 3266, 3447, 3638, 3831, 4028, 4227, 4438, 4661, 4888

Offset: 0

N. J. A. Sloane, Robert G. Wilson v

It appears that a(n)^2 - a(n-1)^2 = A034960(n). - Gary Detlefs, Dec 20 2011
This is true. Proof: By definition we have A034960(n) = Sum_{k = (a(n-1)+1)..a(n)} (2*k-1). Since Sum_{k = 1..n} (2*k-1) = n^2, it follows A034960(n) = a(n)^2 - a(n-1)^2, for n > 1. - Hieronymus Fischer, Sep 27 2012 [formulas above adjusted to changed offset of A034960 - Hieronymus Fischer, Oct 14 2012]
Row sums of the triangle in A037126. - Reinhard Zumkeller, Oct 01 2012
Ramanujan noticed the apparent identity between the prime parts partition numbers A000607 and the expansion of Sum_{k >= 0} x^a(k)/((1-x)...(1-x^k)), cf. A046676. See A192541 for the difference between the two. - M. F. Hasler, Mar 05 2014
For n > 0: row 1 in A254858. - Reinhard Zumkeller, Feb 08 2015
a(n) is the smallest number that can be partitioned into n distinct primes. - Alonso del Arte, May 30 2017
For a(n) < m < a(n+1), n > 0, at least one m is a perfect square.
Proof: For n = 1, 2, ..., 6, the proposition is clear. For n > 6, a(n) < ((prime(n) - 1)/2)^2, set (k - 1)^2 <= a(n) < k^2 < ((prime(n) + 1)/2)^2, then k^2 < (k - 1)^2 + prime(n) <= a(n) + prime(n) = a(n+1), so m = k^2 is this perfect square. - Jinyuan Wang, Oct 04 2018
For n >= 5 we have a(n) < ((prime(n)+1)/2)^2. This can be shown by noting that ((prime(n)+1)/2)^2 - ((prime(n-1)+1)/2)^2 - prime(n) = (prime(n)+prime(n-1))*(prime(n)-prime(n-1)-2)/4 >= 0. - Jianing Song, Nov 13 2022
Washington gives an oscillation formula for |a(n) - pi(n^2)|, see links. - Charles R Greathouse IV, Dec 07 2022

E. Bach and J. Shallit, §2.7 in Algorithmic Number Theory, Vol. 1: Efficient Algorithms, MIT Press, Cambridge, MA, 1996.
H. L. Nelson, "Prime Sums", J. Rec. Math., 14 (1981), 205-206.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

R. J. Mathar, Table of n, a(n) for n = 0..100000
C. Axler, On a Sequence involving Prime Numbers, J. Int. Seq. 18 (2015) # 15.7.6.
Christian Axler, New bounds for the sum of the first n prime numbers, arXiv:1606.06874 [math.NT], 2016.
P. Hecht, Post-Quantum Cryptography: S_381 Cyclic Subgroup of High Order, International Journal of Advanced Engineering Research and Science (IJAERS, 2017) Vol. 4, Issue 6, 78-86.
R. J. Mathar, Table of 100000n, a(100000n) for n = 1..10000
Romeo Meštrović, Curious conjectures on the distribution of primes among the sums of the first 2n primes, arXiv:1804.04198 [math.NT], 2018.
Vladimir Shevelev, Asymptotics of sum of the first n primes with a remainder term
Nilotpal Kanti Sinha, On the asymptotic expansion of the sum of the first n primes, arXiv:1011.1667 [math.NT], 2010-2015.
Lawrence C. Washington, Sums of Powers of Primes II, arXiv preprint (2022). arXiv:2209.12845 [math.NT]
Eric Weisstein's World of Mathematics, Prime Sums
OEIS Wiki, Sums of powers of primes divisibility sequences

Cf. A000041, A034386, A111287, A013916, A013918 (primes), A045345, A050247, A050248, A068873, A073619, A034387, A014148, A014150, A178138, A254784, A254858.
See A122989 for the value of Sum_{n >= 1} 1/a(n).
Cf. A008472, A002110, A038107.

P:=Filtered([1..250],IsPrime);;
a:=Concatenation([0],List([1..Length(P)],i->Sum([1..i],k->P[k]))); # Muniru A Asiru, Oct 07 2018

a007504 n = a007504_list !! n
a007504_list = scanl (+) 0 a000040_list
-- Reinhard Zumkeller, Oct 01 2014, Oct 03 2011

[0] cat [&+[ NthPrime(k): k in [1..n]]: n in [1..50]]; // Bruno Berselli, Apr 11 2011 (adapted by Vincenzo Librandi, Nov 27 2015 after Hasler's change on Mar 05 2014)

s1:=[2]; for n from 2 to 1000 do s1:=[op(s1),s1[n-1]+ithprime(n)]; od: s1;
A007504 := proc(n)
    add(ithprime(i), i=1..n) ;
end proc: # R. J. Mathar, Sep 20 2015

Accumulate[Prime[Range[100]]] (* Zak Seidov, Apr 10 2011 *)
primeRunSum = 0; Table[primeRunSum = primeRunSum + Prime[k], {k, 100}] (* Zak Seidov, Apr 16 2011 *)

A007504(n) = sum(k=1,n,prime(k)) \\ Michael B. Porter, Feb 26 2010

a(n) = vecsum(primes(n)); \\ Michel Marcus, Feb 06 2021

from itertools import accumulate, count, islice
from sympy import prime
def A007504_gen(): return accumulate(prime(n) if n > 0 else 0 for n in count(0))
A007504_list = list(islice(A007504_gen(),20)) # Chai Wah Wu, Feb 23 2022

a(n) ~ n^2 * log(n) / 2. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 24 2001 (see Bach & Shallit 1996)
a(n) = A014284(n+1) - 1. - Jaroslav Krizek, Aug 19 2009
a(n+1) - a(n) = A000040(n+1). - Jaroslav Krizek, Aug 19 2009
a(A051838(n)) = A002110(A051838(n)) / A116536(n). - Reinhard Zumkeller, Oct 03 2011
a(n) = min(A068873(n), A073619(n)) for n > 1. - Jonathan Sondow, Jul 10 2012
a(n) = A033286(n) - A152535(n). - Omar E. Pol, Aug 09 2012
For n >= 3, a(n) >= (n-1)^2 * (log(n-1) - 1/2)/2 and a(n) <= n*(n+1)*(log(n) + log(log(n))+ 1)/2. Thus a(n) = n^2 * log(n) / 2 + O(n^2*log(log(n))). It is more precise than in Fares's comment. - Vladimir Shevelev, Aug 01 2013
a(n) = (n^2/2)*(log n + log log n - 3/2 + (log log n - 3)/log n + (2 (log log n)^2 - 14 log log n + 27)/(4 log^2 n) + O((log log n/log n)^3)) [Sinha]. - Charles R Greathouse IV, Jun 11 2015
G.f: (x*b(x))/(1-x), where b(x) is the g.f. of A000040. - Mario C. Enriquez, Dec 10 2016
a(n) = A008472(A002110(n)), for n > 0. - Michel Marcus, Jul 16 2020

More terms from Stefan Steinerberger, Apr 11 2006

a(0) = 0 prepended by M. F. Hasler, Mar 05 2014

A051838 Numbers k such that sum of first k primes divides product of first k primes.

Original entry on oeis.org

1, 3, 8, 13, 23, 38, 39, 41, 43, 48, 50, 53, 56, 57, 58, 66, 68, 70, 73, 77, 84, 90, 94, 98, 126, 128, 134, 140, 143, 145, 149, 151, 153, 157, 160, 164, 167, 168, 172, 174, 176, 182, 191, 194, 196, 200, 210, 212, 215, 217, 218, 219, 222, 225, 228, 229

Offset: 1

G. L. Honaker, Jr., Dec 12 1999

nonn

			Sum of first 8 primes is 77 and product of first 8 primes is 9699690. 77 divides 9699690 therefore a(3)=8.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)

Cf. A007504, A002110, A159639, A196415.

A116536 gives the quotients, A140763 the divisors and A159578 the dividends.

P:=Filtered([1..2000],IsPrime);;
Filtered([1..Length(P)],n->Product([1..n],i->P[i]) mod Sum([1..n],i->P[i])=0); # Muniru A Asiru, Dec 20 2018

import Data.List (elemIndices)
a051838 n = a051838_list !! (n-1)
a051838_list =
   map (+ 1) $ elemIndices 0 $ zipWith mod a002110_list a007504_list
-- Reinhard Zumkeller, Oct 03 2011

p = Prime@ Range@ 250; Flatten@ Position[ Mod[ First@#, Last@#] & /@ Partition[ Riffle[ Rest[ FoldList[ Times, 1, p]], Accumulate@ p], 2], 0] (* Harvey P. Dale, Dec 19 2010 *)

for(n=1,100,P=prod(i=1,n,prime(i));S=sum(i=1,n,prime(i));if(!(P%S),print1(n,", "))) \\ Derek Orr, Jul 19 2015

isok(n) = my(p = primes(n)); (vecprod(p) % vecsum(p)) == 0; \\ Michel Marcus, Dec 20 2018

From Reinhard Zumkeller, Oct 03 2011: (Start)
A002110(a(n)) mod A007504(a(n)) = 0.
A116536(n) = A002110(a(n)) / A007504(a(n)). (End)

A141092 Product of first k composite numbers divided by their sum, when the result is an integer.

Original entry on oeis.org

1, 64, 46080, 111974400, 662171811840, 310393036800000, 7230916185292800, 108238138194410864640000, 23835710455777670400935290994688000000000, 1104077556971139123493322971152384000000000

Offset: 1

Enoch Haga, Jun 01 2008

Find the products and sums of first k composites, k = 1, 2, 3, .... When the products divided by the sums produce integral quotients, add terms to sequence.

			a(3)=46080 because 4*6*8*9*10*12*14=2903040 and 4+6+8+9+10+12+14=63; 2903040/63=46080, which is an integer, so 46080 is a term.

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

Cf. A196415, A141089, A141090, A141091.
Compare with A140761, A159578, A140763, A116536.
Cf. A116536.

import Data.Maybe (catMaybes)
a141092 n = a141092_list !! (n-1)
a141092_list = catMaybes $ zipWith div' a036691_list a053767_list where
   div' x y | m == 0    = Just x'
            | otherwise = Nothing where (x',m) = divMod x y
-- Reinhard Zumkeller, Oct 03 2011

With[{cnos=Select[Range[50],CompositeQ]},Select[Table[Fold[ Times,1,Take[ cnos,n]]/ Total[Take[cnos,n]],{n,Length[cnos]}],IntegerQ]] (* Harvey P. Dale, Jan 14 2015 *)

s=0;p=1;forcomposite(n=4,100,p*=n;s+=n;if(p%s==0,print1(p/s", "))) \\ Charles R Greathouse IV, Apr 04 2013

a(n) = A036691(A196415(n)) / A053767(A196415(n)). [Reinhard Zumkeller, Oct 03 2011]

Checked by N. J. A. Sloane, Oct 02 2011.

A140763 A051838 gives numbers m such that the sum of first m primes divides the product of the first m primes. This sequence gives corresponding values of the sum of first m primes.

Original entry on oeis.org

2, 10, 77, 238, 874, 2747, 2914, 3266, 3638, 4661, 5117, 5830, 6601, 6870, 7141, 9523, 10191, 10887, 11966, 13490, 16401, 19113, 21037, 23069, 40313, 41741, 46191, 50887, 53342, 54998, 58406, 60146, 61910, 65534, 68341, 72179, 75130, 76127, 80189, 82253

Offset: 1

Enoch Haga, May 28 2008

nonn

Sums (divisors) associated with A140761.

			a(2)=10 because when 30 is divided by 10, the quotient is 3 and integral.

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

Cf. A051838, A116536, A140761, A159578, A116536.

Module[{nn=200,s,p},s=Accumulate[Prime[Range[nn]]];p=FoldList[ Times,Prime[ Range[nn]]];Select[Thread[{p,s}],Divisible[#[[1]],#[[2]]]&]][[All,2]] (* Harvey P. Dale, Jun 07 2022 *)

a(n)=A116536(n)/A159578(n) = A007504(A051838(n)). - R. J. Mathar, Jun 09 2008

Sum_{i=1..A051838(n)} prime(i).

Corrected and edited by N. J. A. Sloane, Oct 01 2011

A140761 Primes p(j) = A000040(j), j>=1, such that p(1)p(2)...*p(j) is an integral multiple of p(1)+p(2)+...+p(j).

Original entry on oeis.org

2, 5, 19, 41, 83, 163, 167, 179, 191, 223, 229, 241, 263, 269, 271, 317, 337, 349, 367, 389, 433, 463, 491, 521, 701, 719, 757, 809, 823, 829, 859, 877, 883, 919, 941, 971, 991, 997, 1021, 1033, 1049, 1091, 1153, 1181, 1193, 1223, 1291, 1301, 1319, 1327, 1361

Offset: 1

Enoch Haga, May 28 2008

			a(2) = 5 because it is the last consecutive prime in the run 2*3*5 = 30 and 2+3+5 = 10; since 30/10 = 3, it is the first integral quotient.

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

Cf. A002110, A007504, A051838, A116536, A140763, A159578.

seq = {}; sum = 0; prod = 1; p = 1; Do[p = NextPrime[p]; prod *= p; sum += p; If[Divisible[prod, sum], AppendTo[seq, p]], {200}]; seq (* Amiram Eldar, Nov 02 2020 *)

Find integral quotients of products of consecutive primes divided by their sum.

a(n) = A000040(A051838(n)). - R. J. Mathar, Jun 09 2008

Edited by R. J. Mathar, Jun 09 2008

a(1) added by Amiram Eldar, Nov 02 2020

cp's OEIS Frontend

A159578 Dividend associated with A116536.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A159580 Integral quotients occur consecutively this many times (in sequence associated with A051838 and A116536).

Views

Author

Keywords

Examples

Crossrefs

Extensions

A159581 First divisors at which integral quotients occur consecutively in A116536 (and associated with A159580).

Views

Author

Keywords

Examples

Crossrefs

Extensions

A159577 Erroneous version of A116536.

Views

Author

Keywords

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

A007504 Sum of the first n primes.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A051838 Numbers k such that sum of first k primes divides product of first k primes.

Views

Author

Keywords

Examples

Links

Crossrefs

A140761 Primes p(j) = A000040(j), j>=1, such that p(1)p(2)...*p(j) is an integral multiple of p(1)+p(2)+...+p(j).