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

A309130 Smallest prime factor of A077586(n).

Original entry on oeis.org

7, 127, 2147483647, 170141183460469231731687303715884105727, 47, 338193759479, 231733529, 62914441, 2351, 1399, 295257526626031, 18287, 106937, 863, 4703, 138863, 22590223644617
Offset: 1

Views

Author

Richard N. Smith, Jul 13 2019

Keywords

Comments

A263686 is a subsequence.
Agrees with A263686 in the first four terms, but then the two sequences differ for the first time at n = 5, because prime(5) = 11 is not in A000043.
a(18) = A263686(9) is greater than 1.56*10^17*(2^61-1), see link.
a(n) = A077586(n) iff A077586(n) is prime, A077586(n) is prime for 1 <= n <= 4, but composite for 5 <= n <= 17. The status of A077586(18) = 2^(2^61-1)-1 is unknown. It is conjectured that A077586(n) is composite for all n >= 5.
a(20) = 456959, a(21) = 18384329, a(22) = 198839, a(23) = 2349023, a(24) = A263686(10) is greater than 1.25*10^16*(2^89-1).
Conjecture: All terms are in A122094 (all terms in A263686 are in A122094).
For examples related to that conjecture, see A322568. - Jeppe Stig Nielsen, Aug 29 2019
a(30) = 46559, a(32) = 23671, a(36) = 7151489, a(39) = 4698047, a(41) = 719, a(43) = 1440847, a(45) = 179689, a(47) = 11759383, a(48) = 23602441, a(50) = 9024439, a(51) = 28875361, a(52) = 6301423, a(54) = 2493983, a(56) = 33518137, a(59) = 6727783, a(66) = 95111, a(72) = 1439, a(73) = 99833, a(78) = 38119, a(81) = 26849, a(83) = 8258911, a(86) = 16173559, a(89) = 625343, a(93) = 9743. - Chai Wah Wu, Oct 16 2019

Crossrefs

Programs

  • PARI
    A309130(n)=A020639(2^(2^prime(n)-1)-1) \\ For efficiency, use addprimes([large terms of this sequence]). - M. F. Hasler, Mar 01 2025

Formula

a(n) = A020639(A077586(n)).
a(n) = A049479(A001348(n)). - M. F. Hasler, Mar 01 2025

A077585 a(n) = 2^(2^n - 1) - 1.

Original entry on oeis.org

0, 1, 7, 127, 32767, 2147483647, 9223372036854775807, 170141183460469231731687303715884105727, 57896044618658097711785492504343953926634992332820282019728792003956564819967
Offset: 0

Views

Author

Henry Bottomley, Nov 07 2002

Keywords

Comments

a(2), a(3), a(5) and a(7) are prime; a(11) is not.
Let S be a set of n elements. First we perform a set partition on S. Let SU be the set of all nonempty subsets of S. As is well known, 2^n - 1 is the number of nonempty subsets of a set with n elements (see A000225). That is, 2^n - 1 is the number of elements |SU| of SU. In the second step, we select k elements from SU. We want to know how many different selections are possible. Let W be the resulting set of selections formed from SU. Then the number of elements |W| of W is |W| = Sum_{k=1..2^n-1} binomial(2^n-1,k) = 2^(2^n-1) - 1 = A077585. Example: |W(n)| = a(n=2) = 7, because W = {[[1, 2]], [[1]], [[1, 2], [1]], [[1, 2], [2], [1]], [[1, 2], [2]], [[2], [1]], [[2]]}. - Thomas Wieder, Nov 08 2007
These are so-called double Mersenne numbers, sometimes denoted MM(n) where M = A000225. MM(n) is prime iff M(n) = 2^n-1 is a Mersenne prime exponent (A000043), which isn't possible unless n itself is also in A000043. Primes of this form are called double Mersenne primes MM(p). For all Mersenne exponents between 7 and 61, factors of MM(p) are known. MM(61) is far too large for any currently known primality test, but a distributed search for factors of this and other MM(p) is ongoing, see the doublemersenne.org web site. - M. F. Hasler, Mar 05 2020
This is a strong divisibility sequence, that is, gcd(a(n),a(m)) = a(gcd(n,m)) for n, m >= 1. - Peter Bala, Dec 05 2022

Examples

			a(5) = 2^(2^5 - 1) - 1 = 2^31 - 1 = 2147483647.
		

Crossrefs

Programs

  • Maple
    a:= n-> 2^(2^n-1)-1:
    seq(a(n), n=0..8); # Thomas Wieder, Nov 08 2007
  • Mathematica
    2^(2^Range[0, 9] - 1) - 1 (* Vladimir Joseph Stephan Orlovsky, Jun 22 2011 *)
  • PARI
    a(n)=if(n<1,0,-1+2*(1+a(n-1))^2)
    
  • PARI
    apply( {A077585(n)=1<<(1<M. F. Hasler, Mar 05 2020
    
  • Python
    def A077585(n): return (1<<(1<Chai Wah Wu, Mar 14 2023

Formula

a(n) = A000225(A000225(n)).
a(n) = A058891(n+1) - 1. - corrected by Maurizio De Leo, Feb 25 2015
a(n) = (A001146(n) - 2)/2.
a(n) = A056220(1+a(n-1)).
a(n) = Sum_{k=1..2^n-1} binomial(2^n-1,k). - Thomas Wieder, Nov 08 2007
a(n) = 2*a(n-1)^2 + 4*a(n-1) + 1. - Roderick MacPhee, Oct 05 2012

Extensions

Corrected by Lekraj Beedassy, Jan 02 2007

A103901 Mersenne primes p such that M(p) = 2^p - 1 is also a (Mersenne) prime.

Original entry on oeis.org

3, 7, 31, 127
Offset: 1

Views

Author

Jonathan Sondow, Feb 20 2005

Keywords

Comments

Same as exponents of double Mersenne primes. Only four terms are known.

Examples

			2^2 - 1 = 3 and 2^3 - 1 = 7 are Mersenne primes, so 3 is a member.
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, 3rd ed., Springer-Verlag, NY, 2004, Sec. A3.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 3rd ed., Oxford Univ. Press, 1954, p. 16.
  • P. Ribenboim, The New Book of Prime Number Records, Springer-Verlag, NY, 1996, Chap. 2, Sec. VII.

Crossrefs

A263686 Smallest prime factor of double Mersenne numbers.

Original entry on oeis.org

7, 127, 2147483647, 170141183460469231731687303715884105727, 338193759479, 231733529, 62914441, 295257526626031
Offset: 1

Views

Author

Arkadiusz Wesolowski, Oct 23 2015

Keywords

Comments

A double Mersenne number is a Mersenne number of the form 2^(2^p - 1) - 1, where p is a Mersenne exponent (A000043).
From M. F. Hasler, Feb 28 2025: (Start)
The prime factors of Mersenne numbers 2^q - 1 must be of the form 2*q*k + 1.
The four smallest double Mersenne numbers (p = 2, 3, 5, 7 => q = 3, 7, 31, 127) are prime, so their smallest prime factor is equal to themselves, a(n) = M(q). This is equivalent to k = (2^(q-1)-1)/q, which is almost as large as M(q) itself: k = 1, 9 and 34636833 for the first three terms, and for q = 127, k has just three digits less than M(q) = a(4) itself. The prime p = 11 is not a Mersenne exponent.
The fifth term, a(5) = 2*(2^13-1)*k + 1 with k = 20644229 (which is prime) is the first proper divisor of the respective M(q), as are the next three, corresponding to p = 17, 19 and 31.
For p = 61, M(q) has 694127911065419642 digits, and so far no factor is known, but it is known that it has no factor less than 10^36. (End)

Crossrefs

Cf. A000043, A000668, A001348, A020639, A049479, A077586, A122094. Subsequence of A016047. Subsequence of A309130.

Programs

  • PARI
    forprime(p=2,,q=2^p-1; !ispseudoprime(q) && next(); if(ispseudoprime(2^q-1), print1(2^q-1,", ");next()); forstep(r=2*q+1,+oo,2*q, !ispseudoprime(r) && next(); if(Mod(2,r)^q-1 == 0, print1(r,", ");next(2)))) \\ Jeppe Stig Nielsen, Aug 28 2019

Formula

a(n) = spf(MM(A000043(n))) = A049479(A000668(n)), where spf = A020639 is the smallest prime factor, A049479 = spf o M, M(p) = 2^p-1 = A000225(p), MM = M o M = A077585, A000668(n) = M(A000043(n)), A000043 are the Mersenne prime exponents. - M. F. Hasler, Mar 01 2025

A103902 Mersenne primes p such that the Mersenne number M(p) = 2^p - 1 is composite.

Original entry on oeis.org

8191, 131071, 524287, 2147483647
Offset: 1

Views

Author

Jonathan Sondow, Feb 20 2005

Keywords

Comments

Only four terms are known.
The first four Mersenne primes (p=2^q-1 in A000668) are double Mersenne primes, i.e., in A103901. The next four yield a composite M(p) and therefore are in this sequence. The next larger Mersenne prime p = A000668(9) has already 19 digits and is much too large to enable us, as of today, to test the primality of 2^p-1 (which would require over 10^8 gigabytes just to be stored in binary). This explains that only 4 terms are known of this sequence and of A103901; for all the 30+ remaining members of A000668 it is not known whether they belong to A103901 or to this sequence A103902. - M. F. Hasler, Jan 21 2015

Examples

			M(13) = 8191 is a Mersenne prime and M(1891) is composite, so 1891 is a member.
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, 3rd ed., Springer-Verlag, NY, 2004, Sec. A3.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 3rd ed., Oxford Univ. Press, 1954, p. 16.
  • P. Ribenboim, The New Book of Prime Number Records, Springer-Verlag, NY, 1996, Chap. 2, Sec. VII.

Crossrefs

Programs

A277634 a(n) = 2^A000668(n) - 1.

Original entry on oeis.org

7, 127, 2147483647, 170141183460469231731687303715884105727
Offset: 1

Views

Author

Felix Fröhlich, Oct 24 2016

Keywords

Comments

The terms of this sequence are sometimes called "Double Mersenne numbers" (cf. A263686).
Agrees with A077586 in the first four terms, but then the two sequences differ for the first time at n = 5, because prime(5) = 11 is not in A000043.
a(5) is too large to include in data section (see A276641).
a(n) = A263686(n) iff a(n) is prime, which is the case iff A000668(n) is in A103901.
Agrees with A263686 at least in the first four terms. - Omar E. Pol, Oct 24 2016

Crossrefs

Programs

  • PARI
    forprime(p=1, 11, if(ispseudoprime(2^p-1), print1(2^(2^p-1)-1, ", ")))

Formula

a(n) = A000225(A000668(n)).
Showing 1-6 of 6 results.