cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-7 of 7 results.

A007741 a(n) = prime(n)*...*prime(m), the least product of consecutive primes which is abundant.

Original entry on oeis.org

30, 15015, 33426748355, 1357656019974967471687377449, 7105630242567996762185122555313528897845637444413640621
Offset: 1

Views

Author

Keywords

Comments

Essentially, i.e., except for a(1), identical to A007702. All terms are primitive abundant numbers (A091191) and thus, except for the first term, odd primitive abundant (A006038). The next term is too large to be displayed here, see A007707 (and formula) for many more terms, using a more compact encoding. - M. F. Hasler, Apr 30 2017

Crossrefs

Programs

  • Mathematica
    a[n_] := Module[{p = Prime[n]}, r = 1; prod = 1; While[r <= 2, r *= 1 + 1/p; prod *= p; p = NextPrime[p]]; prod]; Array[a, 5] (* Amiram Eldar, Jun 29 2019 *)
  • PARI
    a(n) = {p = prime(n); sig = p+1; prd = p; while (sig <= 2*prd, p = nextprime(p+1); sig *= p+1; prd *= p;); return (prd);} \\ Michel Marcus, Mar 10 2013

Formula

a(n) = Product_{k=n..A007707(n)} prime(k) = Product_{0 <= i < A108227(n)} prime(n+i). - M. F. Hasler, Apr 30 2017 and Jun 15 2017

Extensions

More terms from Don Reble, Nov 10 2005

A007708 Prime(n)*...*a(n) is the least product of consecutive primes which is abundant.

Original entry on oeis.org

5, 13, 31, 73, 149, 233, 367, 521, 733, 991, 1249, 1579, 1949, 2341, 2791, 3343, 3881, 4481, 5147, 5849, 6619, 7499, 8387, 9341, 10321, 11411, 12517, 13709, 15013, 16363, 17881, 19381, 20873, 22369, 24007, 25763, 27611, 29399, 31357
Offset: 1

Views

Author

Keywords

Comments

Differs from A007686 only for n=1. - Michel Marcus, Mar 10 2013

Crossrefs

Programs

  • Mathematica
    a[n_] := Module[{p = Prime[n]}, r = 1 + 1/p; While[r <= 2,  p = NextPrime[p]; r *= 1 + 1/p]; p]; Array[a, 39] (* Amiram Eldar, Jun 29 2019 *)
  • PARI
    a(n) = {p = prime(n); sig = p+1; prd = p; while (sig <= 2*prd, p = nextprime(p+1); sig *= p+1; prd *= p;); return (p);} \\ Michel Marcus, Mar 10 2013

Extensions

More terms from Don Reble, Nov 10 2005

A108227 a(n) is the least number of prime factors for any abundant number with p_n (the n-th prime) as its least factor.

Original entry on oeis.org

3, 5, 9, 18, 31, 46, 67, 91, 122, 158, 194, 238, 284, 334, 392, 456, 522, 591, 668, 749, 835, 929, 1028, 1133, 1242, 1352, 1469, 1594, 1727, 1869, 2019, 2163, 2315, 2471, 2636, 2802, 2977, 3157, 3342, 3534, 3731, 3933, 4145, 4358, 4581, 4811
Offset: 1

Views

Author

Hugo van der Sanden, Jun 17 2005

Keywords

Comments

If we replace "abundant" in the definition with "non-deficient", we get the same sequence with an initial 2 instead of 3, barring an astronomically unlikely coincidence with some as-yet-undiscovered odd perfect number. [This is sequence A107705. - M. F. Hasler, Jun 14 2017]
It appears that all terms >= 5 correspond to the odd primitive abundant numbers (A006038) which are products of consecutive primes (cf. A285993), i.e., of the form N = Product_{0<=iM. F. Hasler, May 08 2017
From Jianing Song, Apr 21 2021: (Start)
Let x_1 < x_2 < ... < x_k < ... be the numbers of the form p of p^2 + p, where p is a prime >= prime(n). Then a(n) is the smallest N such that Product_{i=1..N} (1 + 1/x_i) > 2. See my link below for a proof.
For example, for n = 3, we have {x_1, x_2, ..., x_k, ...} = {5, 7, 11, 13, 17, 19, 23, 29, 5^2 + 5, ...}, we have Product_{i=1..8} (1 + 1/x_i) < 2 and Product_{i=1..9} (1 + 1/x_i) > 2, so a(3) = 9. (End)

Examples

			a(2) = 5 since 945 = 3^3*5*7 is an abundant number with p_2 = 3 as its smallest prime factor, and no such number exists with fewer than 5 prime factors.
		

Crossrefs

Cf. A107705.
Cf. A001276 (least number of prime factors for a (p_n)-rough abundant number, counted without multiplicity).

Programs

  • PARI
    A108227(n, s=1+1/prime(n))=for(a=1, 9e9, if(2M. F. Hasler, Jun 15 2017
    
  • PARI
    isform(k,q) = my(p=prime(k)); if(isprime(q) && (q>=p), 1, if(issquare(4*q+1), my(r=(sqrtint(4*q+1)-1)/2); isprime(r) && (r>=p), 0))
    a(n) = my(Prod=1, Sum=0); for(i=prime(n), oo, if(isform(n,i), Prod *= (1+1/i); Sum++); if(Prod>2, return(Sum))) \\ Jianing Song, Apr 21 2021

Formula

a(n) = A007684(n)-n+1, for n>1. A007741(n) = Product_{0<=iM. F. Hasler, Jun 15 2017

Extensions

Data corrected by Amiram Eldar, Aug 08 2019

A007684 Prime(n)*...*prime(a(n)) is the least product of consecutive primes that is non-deficient.

Original entry on oeis.org

2, 6, 11, 21, 35, 51, 73, 98, 130, 167, 204, 249, 296, 347, 406, 471, 538, 608, 686, 768, 855, 950, 1050, 1156, 1266, 1377, 1495, 1621, 1755, 1898, 2049, 2194, 2347, 2504, 2670, 2837, 3013, 3194, 3380, 3573, 3771, 3974, 4187, 4401, 4625, 4856
Offset: 1

Views

Author

Keywords

Comments

Subscript of the smallest primorial number that when divided by the (n-1)-th primorial number gives an abundant number.
Products of consecutive primes started with prime(a) up to prime(b) result in abundant squarefree numbers if b is large enough and provides perhaps the least squarefree solutions to Rivera Puzzle 329 and its generalization.
Adding a new prime p to the product increases the relative abundancy sigma(N)/N by a factor 1+1/p. This leads to a simple and fast algorithm, see the PARI code. - M. F. Hasler, Jul 30 2016

Examples

			n=1: a(1)=2 means that primorial(2)=6 divided by primorial(1-1)=1 gives the quotient 6/1=6 which is just non-deficient (being a perfect number);
n=3: a(n)=11 because prime(3)=5, primorial(11) = 2*3*5*...*29*31, primorial(3-1) = 2*3 = 6.
p#(11)/p#(2) = 3*5*7*11*13*17*19*23*29*31 = 33426748355 = q and sigma(q)/q = 2.00097 > 2 so q is an abundant number. Also p#(10)/p#(3-1) is not yet abundant.
		

Crossrefs

Cf. A005100, A007686, A007702, A007707 (an essentially identical sequence).

Programs

  • Mathematica
    spr[x_, y_] :=Apply[Times, Table[(Prime[w]+1)/(Prime[w]), {w, x, y}]];
    Table[Min[Flatten[Position[Table[Floor[spr[n, w]], {w, 1, 1000}], 2]]], {n, 1, 20}] (* Labos Elemer, Sep 19 2005 *)
  • PARI
    a=1;i=0;for(n=1,99,while(2>a*=1+1/prime(i++),);print1(i",");a/=1+1/prime(n)) \\ M. F. Hasler, Jul 30 2016

Formula

a(n) is the minimal x such that floor(sigma(p#(x)/p#(n-1)) / (p#(x)/p#(n-1))) = 2, where p#(w) is the w-th primorial number, the product of first w prime numbers. For a>b, the p#(a)/p#(b)=A002110(a)/A002110(b) quotients are prime(b+1)*prime(b+2)*...*prime(a).

Extensions

Additional comments from Labos Elemer, Sep 19 2005
More terms from Don Reble, Nov 10 2005
Edited by N. J. A. Sloane, Dec 22 2006

A107705 a(n) is the least number of prime factors in any non-deficient number that has the n-th prime as its least prime factor.

Original entry on oeis.org

2, 5, 9, 18, 31, 46, 67, 91, 122, 158, 194, 238, 284, 334, 392, 456, 522, 591, 668, 749, 835, 929, 1028, 1133, 1242, 1352, 1469, 1594, 1727, 1869, 2019, 2163, 2315, 2471, 2636, 2802, 2977, 3157, 3342, 3534, 3731, 3933, 4145, 4358, 4581, 4811, 5053, 5293
Offset: 1

Views

Author

Hugo van der Sanden, Jun 10 2005

Keywords

Comments

Barring unforeseen odd perfect numbers (which it has been proved must have at least 29 prime factors if they exist at all), if we replace "non-deficient" in the description with "abundant", the value of a(1) becomes 3 and all other values stay the same.
The above mentioned sequence is A108227, see there for a comment on the relation of this sequence to that of primitive abundant numbers (A006038) which are products of consecutive primes, i.e., of the form N = Product_{0<=iA007702. - M. F. Hasler, Jun 15 2017

Examples

			a(2) is 5 since 1) there are abundant numbers with a(2)=5 prime factors of which p_2=3 is the least prime factor (such as 945 = 3^3.5.7); 2) there are no non-deficient numbers with fewer than 5 prime factors, of which 3 is the least prime factor.
		

Crossrefs

Programs

  • PARI
    A107705(n,s=1+1/prime(n))=for(a=1,9e9,2>(s*=1+1/prime(n+a))||return(a+1)) \\ M. F. Hasler, Jun 15 2017

Formula

a(n) = A007684(n)-n+1. A007702(n) = Product_{0<=iM. F. Hasler, Jun 15 2017

Extensions

Data corrected by Amiram Eldar, Aug 08 2019

A285993 Largest odd abundant number (A005231) equal to the product of n consecutive primes.

Original entry on oeis.org

15015, 255255, 4849845, 111546435, 33426748355, 1236789689135, 50708377254535, 2180460221945005, 102481630431415235, 5431526412865007455, 320460058359035439845, 19548063559901161830545, 1309720258513377842646515, 1357656019974967471687377449, 107254825578022430263302818471
Offset: 5

Views

Author

M. F. Hasler, Apr 30 2017

Keywords

Comments

The smallest term is a(5) = 3*5*7*11*13, there is no odd abundant number (A005231) equal to the product of less than 5 consecutive primes.
The smallest odd abundant number (A005231) equal to the product of n consecutive primes is equal (when it exists, i.e., for n >= 5) to the least odd number with n (distinct) prime divisors, equal to the product of the first n odd primes = A070826(n+1) = A002110(n+1)/2.
See A188342 = (945, 3465, 15015, 692835, 22309287, ...) for the least odd primitive abundant number (A006038) with n distinct prime factors, and A275449 for the least odd primitive abundant number with n prime factors counted with multiplicity.
The terms are in general not primitive abundant numbers (A091191), in particular this cannot be the case when a(n) is a multiple of a(n-1), as is the case for most of the terms, for which a(n) = a(n-1)*A117366(a(n-1)). In the other event, spf(a(n)) = nextprime(spf(a(n-1))), and a(n) is in A007741(2,3,4...). These are exactly the primitive terms in this sequence.

Examples

			For n < 5, there is no odd abundant number equal to the product of n distinct primes.
For 5 <= n <= 8, the largest odd abundant number equal to the product of n consecutive primes is 3*...*prime(n+1).
For 9 <= n <= 17, the largest odd abundant number equal to the product of n consecutive primes is 5*...*prime(n+2).
For 18 <= n <= 30, the largest odd abundant number equal to the product of n consecutive primes is 7*...*prime(n+3).
For 31 <= n <= 45, the largest odd abundant number equal to the product of n consecutive primes is 11*...*prime(n+4).
For 46 <= n <= 66, the largest odd abundant number equal to the product of n consecutive primes is 13*...*prime(n+5).
		

Crossrefs

A subsequence of A112643 (odd squarefree abundant numbers); see also A108227 (~ A107705) which give indices of primitive terms = those with smallest prime factor larger than that of earlier terms.

Programs

  • PARI
    a(r,f=vector(r,i,prime(i+1)),o)={ while(sigma(factorback(f),-1)>2, o=f; f=concat(f[^1],nextprime(f[r]+1)));factorback(o)} \\ Intentionally throws an error when n < 5.

Formula

a(n) >= a(n-1)*p where p = A117366(a(n-1)) = A151800(A006530(a(n-1))) = nextprime(gpf(a(n-1))), an odd abundant number equal to the product of n consecutive primes. We have strict inequality for n = 9, 18, 31, 46, 67, ..., in which case a(n) = a(n-1)*p*p'/q, where p' = nextprime(p), q = least prime factor of a(n-1). This is the case if a(n) is in A007741.

A286042 Largest prime factor of A285993(n), the largest odd abundant number (A005231) equal to the product of n consecutive primes.

Original entry on oeis.org

13, 17, 19, 23, 31, 37, 41, 43, 47, 53, 59, 61, 67, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293, 307, 311, 313, 317, 331, 337, 347, 349, 353
Offset: 5

Views

Author

M. F. Hasler, May 01 2017

Keywords

Comments

The smallest term is a(5), there is no odd abundant number (A005231) equal to the product of less than 5 consecutive primes.
The corresponding abundant numbers are A285993(n) = prime(k-n+1)*...*prime(k), with prime(k) = a(n).

Examples

			For n < 5, there is no odd abundant number equal to the product of n distinct primes.
For 5 <= n <= 8, the largest odd abundant number equal to the product of n consecutive primes is 3*...*a(n) with a(n) = prime(n+1).
For 9 <= n <= 17, the largest odd abundant number equal to the product of n consecutive primes is 5*...*a(n) with a(n) = prime(n+2).
For 18 <= n <= 30, the largest odd abundant number equal to the product of n consecutive primes is 7*...*a(n) with a(n) = prime(n+3).
For 31 <= n <= 45, the largest odd abundant number equal to the product of n consecutive primes is 11*...*a(n) with a(n) = prime(n+4).
For 46 <= n <= 66, the largest odd abundant number equal to the product of n consecutive primes is 13*...*a(n) with a(n) = prime(n+5).
		

Crossrefs

Programs

  • PARI
    a(r,f=vector(r,i,prime(i+1)),o)={ while(sigma(factorback(f),-1)>2, o=f; f=concat(f[^1],nextprime(f[r]+1)));o[#o]} \\ Intentionally throws an error when n < 5.

Formula

a(n) = A006530(A285993(n)) >= A151800(a(n-1)) = nextprime(a(n-1)), with strict inequality for n = 9, 18, 31, 46, 67, ..., in which case a(n) = nextprime(nextprime(a(n-1))). This is the case if A285993(n) is in A007741.

Extensions

a(66) corrected by Amiram Eldar, Sep 24 2019
Showing 1-7 of 7 results.