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-5 of 5 results.

A047802 Smallest abundant number (sigma(x) > 2x) which is not divisible by any of the first n primes.

Original entry on oeis.org

12, 945, 5391411025, 20169691981106018776756331, 49061132957714428902152118459264865645885092682687973, 7970466327524571538225709545434506255970026969710012787303278390616918473506860039424701
Offset: 0

Views

Author

Ulrich Schimke (ulrschimke(AT)aol.com)

Keywords

Comments

a(n) exists for every n, since the sum of the inverses of the primes is infinite.
From Fred Schneider, Sep 20 2006; edited by Danny Rorabaugh, Nov 26 2018: (Start)
Heuristic: Add the squares of several successive primes and then add successive primes until the number is abundant.
a(2) = 5^2 * 7 * 11 * 13 * 17 * 19 * 23 * 29;
a(3) = 7^2 * 11^2 * 13 * 17 * ... * 61 * 67;
a(4) = 11^2 * 13^2 * 17 * 19 * ... * 131 * 137;
a(5) = 13^2 * 17^2 * 19 * 23 * ... * 223 * 227. (End)
a(6) = 17^2 * 19^2 * 23^2 * 29 * 31 * ... * 347 * 349;
a(7) = 19^2 * 23^2 * 29^2 * 31 * 37 * ... * 491 * 499 (both coming from the D. Iannucci paper). - Michel Marcus, May 01 2013
The known terms of this sequence provide Egyptian decompositions of unity in which all the denominators lack the first n primes, as follows: Every term listed in this sequence is a semiperfect number, which means that a subset of its divisors add up to the number itself. The decomposition 1 = 1/a + 1/b + ... + 1/m, where the denominators are a(n) divided by those divisors, is the desired decomposition. - Javier Múgica, Nov 15 2017
a(n) is the product of consecutive primes starting from prime(n+1) raised to nonincreasing powers. - Jianing Song, Apr 10 2021
From Jianing Song, Apr 14 2021: (Start)
By definition, Omega(a(n)) >= A108227(n+1) for all n, where Omega = A001222. For 0 <= n <= 12 we have Omega(a(n)) = A108227(n+1), but this is not true for n = 13, where Omega(a(13)) = 335 > A108227(14) = 334.
We also have omega(a(n)) >= A001276(n+1) for all n, where omega = A001221. The differences for known terms are 0, 0, 1, 1, 2, 3, 2, 3, 4, 4, 5, 6, 6, 6 respectively.
Conjecture: other than a(1) = 945, all terms are cubefree. (End)

Examples

			a(0) = 12, the first abundant number; a(1) = 945, the first odd abundant number; a(5) is the first abundant number not divisible by 2,3,5,7 or 11.
		

References

  • M. T. Whalen and C. L. Miller, Odd abundant numbers: some interesting observations, Journal of Recreational Mathematics 22 (1990), pp. 257-261.

Crossrefs

Subsequence of A005101 and A133812; cf. A005231.

Formula

Iannucci shows that log a(n) = (n log n)^(2 + o(1)). - Charles R Greathouse IV, Feb 16 2011

Extensions

2 more terms from Fred Schneider, Sep 20 2006

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

A001275 Smallest prime p such that the product of q/(q-1) over the primes from prime(n) to p is greater than 2.

Original entry on oeis.org

3, 7, 23, 61, 127, 199, 337, 479, 677, 937, 1193, 1511, 1871, 2267, 2707, 3251, 3769, 4349, 5009, 5711, 6451, 7321, 8231, 9173, 10151, 11197, 12343, 13487, 14779, 16097, 17599, 19087, 20563, 22109, 23761, 25469, 27259, 29123, 31081, 33029
Offset: 1

Views

Author

Keywords

Comments

A perfect (or abundant) number with prime(n) as its lowest prime factor must be divisible by a prime greater than or equal to a(n).

References

  • 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).

Crossrefs

Cf. A001276.

Programs

  • Mathematica
    a[n_] := Module[{p = If[n == 1, 1, Prime[n - 1]], r = 1}, While[r <= 2, p = NextPrime[p]; r *= p/(p - 1)]; p]; Array[a, 50]  (* Amiram Eldar, Jul 12 2019 *)
  • PARI
    a(n)=my(pr=1.);forprime(p=prime(n),default(primelimit),pr*=p/(p-1);if(pr>2,return(p))) \\ Charles R Greathouse IV, May 09 2011

Formula

a(n) = prime(n)^2 + O(n^2/exp((log n)^(4/7 - e))) for any e > 0.
a(n) = prime(A001276(n) + n - 1). - Amiram Eldar, Jul 12 2019

Extensions

Comment, formula, program, and new definition from Charles R Greathouse IV, May 09 2011

A340194 a(n) = A337544(A003961(n)).

Original entry on oeis.org

1, 3, 5, 15, 9, 13, 11, 75, 35, 25, 15, 65, 17, 31, 43, 375, 21, 91, 27, 125, 53, 43, 29, 325, 99, 49, 245, 155, 35, 95, 39, 1875, 73, 61, 97, 455, 41, 79, 83, 625, 45, 121, 51, 215, 301, 85, 57, 1625, 143, 275, 103, 245, 59, 637, 133, 775, 133, 103, 65, 475, 69, 115, 371, 9375, 151, 173, 71, 305, 143, 245, 77, 2275, 81, 121
Offset: 1

Views

Author

Antti Karttunen, Dec 31 2020

Keywords

Comments

First negative term is a(510510) = -686785. See A001276.

Crossrefs

Programs

Formula

A343357 7-rough abundant numbers.

Original entry on oeis.org

20169691981106018776756331, 21373852696395930345517903, 21975933054040886129898689, 23476198863254546445077041, 23782174126975753483041047, 23836908704943476736166573, 24137500239684251978741183, 24272002214551310731350839, 24955720586792192723783257, 24986334842265665051802619
Offset: 1

Views

Author

David A. Corneth, Apr 12 2021

Keywords

Comments

Each term has at least A001276(4) = 15 distinct prime factors and A108227(4) = 18 prime factors counted with multiplicity. - Jianing Song, Apr 13 2021
The smallest term with exactly 15 distinct prime factors is a(830) = 465709156638373299218537971 = 7^3 * 11^2 * 13^2 * 17^2 * 19 * 23 * ... * 61. - Jianing Song, Apr 14 2021

Examples

			k = 20169691981106018776756331 is in the sequence as its smallest prime factor is at least 7 and it is abundant as sigma(k) > 2*k.
		

Crossrefs

Programs

  • PARI
    is(n) = gcd(n, 30) == 1 && sigma(n) > 2*n
Showing 1-5 of 5 results.