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-10 of 28 results. Next

A052340 Record numbers of iterations reached in A050412.

Original entry on oeis.org

1, 2, 3, 4, 7, 24, 2552, 800516
Offset: 1

Views

Author

Christian G. Bower, Dec 23 1999

Keywords

Crossrefs

Formula

a(n) = A050412(A052339(n)). - Amiram Eldar, Jul 25 2019

Extensions

a(8) from the b-file at A050412 added by Amiram Eldar, Jul 25 2019

A101036 Riesel numbers n (n*2^k-1 is composite for all k>0, n odd) that have a covering set.

Original entry on oeis.org

509203, 762701, 777149, 790841, 992077, 1106681, 1247173, 1254341, 1330207, 1330319, 1715053, 1730653, 1730681, 1744117, 1830187, 1976473, 2136283, 2251349, 2313487, 2344211, 2554843, 2924861, 3079469, 3177553, 3292241, 3419789, 3423373, 3580901
Offset: 1

Views

Author

David W. Wilson, Jan 17 2005

Keywords

Comments

Conjecture: there are infinitely many Riesel numbers that do not arise from a covering system. See page 16 of the Filaseta et al. reference. - Arkadiusz Wesolowski, Nov 17 2014
a(1) = 509203 is also the smallest odd n for which either n^p*2^k - 1 or abs(n^p - 2^k) is composite for every k > 0 and every prime p > 3. - Arkadiusz Wesolowski, Oct 12 2015
Theorem 11 of Filaseta et al. gives a Riesel number which is thought to violate the assumption of a periodic sequence of prime divisors mentioned in the title of this sequence. - Jeppe Stig Nielsen, Mar 16 2019
If the Riesel number mentioned in the previous comment does in fact not have a covering set, then this sequence is different from A076337, because then that number, 3896845303873881175159314620808887046066972469809^2, is a term of A076337, but not of this sequence. - Felix Fröhlich, Sep 09 2019
Named after the Swedish mathematician Hans Ivar Riesel (1929-2014). - Amiram Eldar, Jun 20 2021
Conjecture: if R is a Riesel number (that has a covering set), then there exists a prime P such that R^p is also a Riesel number for every prime p > P. - Thomas Ordowski, Jul 12 2022
Problem: are there numbers K such that K + 2^m is a Riesel number for every m > 0? If so, then (K + 2^m)*2^n - 1 is composite for every pair of positive integers m,n. Also, by the dual Riesel conjecture, |K + 2^m - 2^n| are always composite. Note that, by the dual Riesel conjecture, if p is an odd prime and n is a positive integer, then there exists n such that (p + 2^m)*2^n - 1 is prime. So if such a number K exists, it must be composite. - Thomas Ordowski, Jul 20 2022

References

  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 238.

Crossrefs

Main sequences for Riesel problem: A038699, A040081, A046069, A050412, A052333, A076337, A101036, A108129.
See A076337 for references and additional information. Cf. A076336.

Extensions

Up to 3292241, checked by Don Reble, Jan 17 2005, who comments that up to this point each n*2^k-1 has a prime factor <= 241.
New name from Felix Fröhlich, Sep 09 2019

A076337 Riesel numbers: odd numbers n such that for all k >= 1 the numbers n*2^k - 1 are composite.

Original entry on oeis.org

509203
Offset: 1

Views

Author

N. J. A. Sloane, Nov 07 2002

Keywords

Comments

509203 has been proved to be a member of the sequence, and is conjectured to be the smallest member. However, as of 2009, there are still several smaller numbers which are candidates and have not yet been ruled out (see links).
Riesel numbers are proved by exhibiting a periodic sequence p of prime divisors with p(k) | n*2^k-1 and disproved by finding prime n*2^k-1. It is conjectured that numbers that cannot be proved Riesel in this way are non-Riesel. However, some numbers resist both proof and disproof.
Others conjecture the opposite: that there are infinitely many Riesel numbers that do not arise from a covering system, see A101036. The word "odd" is needed in the definition because otherwise for any term n, all numbers n*2^m, m >= 1, would also be Riesel numbers, but we don't want them in this sequence (as is manifest from A101036). Since 1 and 3 obviously are not in this sequence, for any n in this sequence n-1 is an even number > 2 and therefore composite, so one could replace "k >= 1" equivalently by "k >= 0". - M. F. Hasler, Aug 20 2020
Named after the Swedish mathematician Hans Ivar Riesel (1929-2014). - Amiram Eldar, Apr 02 2022

References

  • R. K. Guy, Unsolved Problems in Number Theory, Section B21.
  • Paulo Ribenboim, The Book of Prime Number Records, 2nd ed., 1989, p. 282.

Crossrefs

Main sequences for Riesel problem: A040081, A046069, A050412, A052333, A076337, A101036, A108129.

Extensions

Normally we require at least four terms but we will make an exception for this sequence in view of its importance. - N. J. A. Sloane, Nov 07 2002. See A101036 for the most likely extension.
Edited by N. J. A. Sloane, Nov 13 2009
Definition corrected ("odd" added) by M. F. Hasler, Aug 23 2020

A040081 Riesel problem: a(n) = smallest m >= 0 such that n*2^m-1 is prime, or -1 if no such prime exists.

Original entry on oeis.org

2, 1, 0, 0, 2, 0, 1, 0, 1, 1, 2, 0, 3, 0, 1, 1, 2, 0, 1, 0, 1, 1, 4, 0, 3, 2, 1, 3, 4, 0, 1, 0, 2, 1, 2, 1, 1, 0, 3, 1, 2, 0, 7, 0, 1, 3, 4, 0, 1, 2, 1, 1, 2, 0, 1, 2, 1, 3, 12, 0, 3, 0, 2, 1, 4, 1, 5, 0, 1, 1, 2, 0, 7, 0, 1, 1, 2, 2, 1, 0, 3, 1, 2, 0, 5, 6, 1, 23, 4, 0, 1, 2, 3, 3, 2, 1, 1, 0, 1, 1, 10, 0, 3
Offset: 1

Views

Author

Keywords

Crossrefs

Main sequences for Riesel problem: A038699, A040081, A046069, A050412, A052333, A076337, A101036, A108129.

Programs

  • Haskell
    a040081 = length . takeWhile ((== 0) . a010051) .
                           iterate  ((+ 1) . (* 2)) . (subtract 1)
    -- Reinhard Zumkeller, Mar 05 2012
    
  • Mathematica
    Table[m = 0; While[! PrimeQ[n*2^m - 1], m++]; m, {n, 100}] (* Arkadiusz Wesolowski, Sep 04 2011 *)
  • PARI
    a(n)=for(k=0,2^16,if(ispseudoprime(n*2^k-1), return(k))) \\ Eric Chen, Jun 01 2015
    
  • Python
    from sympy import isprime
    def a(n):
      m = 0
      while not isprime(n*2**m - 1): m += 1
      return m
    print([a(n) for n in range(1, 88)]) # Michael S. Branicky, Feb 01 2021

A052333 Riesel problem: start with n; repeatedly double and add 1 until reach a prime. Sequence gives a(n) = prime reached, or 0 if no prime is ever reached.

Original entry on oeis.org

3, 5, 7, 19, 11, 13, 31, 17, 19, 43, 23, 103, 223, 29, 31, 67, 71, 37, 79, 41, 43, 367, 47, 199, 103, 53, 223, 463, 59, 61, 127, 131, 67, 139, 71, 73, 151, 311, 79, 163, 83, 5503, 738197503, 89, 367, 751, 191, 97, 199, 101, 103, 211, 107, 109, 223, 113, 463
Offset: 1

Views

Author

Christian G. Bower, Dec 19 1999

Keywords

Comments

Equivalently, a(n) = smallest prime of form (n+1)*2^k-1 for k >= 1, or 0 if no such prime exists.
a(509202) = 0 (i.e. never reaches a prime) - Chris Nash (chris_nash(AT)hotmail.com). (Of course this is related to the entry 509203 of A076337.)
a(73) is a 771-digit prime reached after 2552 iterations - Warut Roonguthai. This was proved to be a prime by Paul Jobling (Paul.Jobling(AT)WhiteCross.com) using PrimeForm and by Ignacio Larrosa Cañestro using Titanix (http://www.znz.freesurf.fr/pages/titanix.html). [Oct 30 2000]

Examples

			a(4)=19 because 4 -> 9 (composite) -> 19 (prime).
		

Crossrefs

CMain sequences for Riesel problem: A038699, A040081, A046069, A050412, A052333, A076337, A101036, A108129.

Programs

  • Mathematica
    Table[NestWhile[2#+1&,2n+1,!PrimeQ[#]&,1,1000],{n,60}] (* Harvey P. Dale, May 08 2011 *)
  • PARI
    a(n)=while(!isprime(n=2*n+1),);n \\ oo loop when a(n) = 0. - Charles R Greathouse IV, May 08 2011

A038699 Riesel problem: Smallest prime of form n*2^m-1, m >= 0, or 0 if no such prime exists.

Original entry on oeis.org

3, 3, 2, 3, 19, 5, 13, 7, 17, 19, 43, 11, 103, 13, 29, 31, 67, 17, 37, 19, 41, 43, 367, 23, 199, 103, 53, 223, 463, 29, 61, 31, 131, 67, 139, 71, 73, 37, 311, 79, 163, 41, 5503, 43, 89, 367, 751, 47, 97, 199, 101, 103, 211, 53, 109, 223, 113, 463, 241663, 59, 487, 61
Offset: 1

Views

Author

N. J. A. Sloane, Dec 30 1999

Keywords

Crossrefs

Primes arising in A040081 (or 0 if no prime exists).
Main sequences for Riesel problem: A038699, A040081, A046069, A050412, A052333, A076337, A101036, A108129.

Programs

  • Haskell
    a038699 = until ((== 1) . a010051) ((+ 1) . (* 2)) . (subtract 1)
    -- Reinhard Zumkeller, Mar 05 2012
  • Mathematica
    getm[n_]:=Module[{m=0},While[!PrimeQ[n 2^m-1],m++];n 2^m-1]; Array[getm,80]  (* Harvey P. Dale, Apr 24 2011 *)

Extensions

More terms from Henry Bottomley, Apr 24 2001

A067760 a(n) is the least positive k such that (2n+1) + 2^k is prime, or 0 if no such k exists.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 3, 2, 1, 1, 4, 2, 1, 2, 1, 1, 2, 1, 5, 2, 1, 3, 2, 1, 1, 8, 2, 1, 2, 1, 1, 4, 2, 1, 2, 1, 7, 2, 1, 3, 4, 2, 1, 2, 1, 1, 2, 1, 1, 2, 1, 7, 4, 5, 3, 4, 2, 1, 2, 1, 3, 2, 1, 1, 10, 3, 3, 2, 1, 1, 4, 2, 1, 4, 2, 1, 2, 1, 5, 2, 1, 3, 2, 1, 1, 4, 3, 3, 2, 1, 1, 2, 1, 1, 6, 5, 3, 6
Offset: 0

Views

Author

Don Reble, Feb 05 2002

Keywords

Comments

From Phil Moore (moorep(AT)lanecc.edu), Dec 14 2009: (Start)
It is known that a(39278) = 0, since no such prime exists for the Sierpiński number 78557 (cf. A076336).
It has recently been discovered that 2131+2^4583176 and 41693+2^5146295 are probable primes, so a(1065) is probably 4583176 and a(20846) is probably 5146295.
At present, the only odd value less than 78557 for which no prime or strong probable prime of the form t+2^k is known is t = 40291, so a(20145) is completely unknown. In addition, for 25 values of t < 78557, only strong probable primes are known. (End)
The last case was resolved in 2011 when the probable prime 40291+2^9092392 was found as a part of a distributed project "Five or Bust". See links. - Jeppe Stig Nielsen, Mar 29 2019

Examples

			a(15)=4 because (2*15+1)+2^k is composite for k=1,2,3 and prime for k=4.
		

Crossrefs

Programs

  • PARI
    a(n) = {my(k=1); while (! isprime((2*n+1)+2^k), k++); k;} \\ Michel Marcus, Feb 26 2018

A046069 Riesel Problem: Smallest m >= 0 such that (2n-1)2^m-1 is prime, or -1 if no such value exists.

Original entry on oeis.org

2, 0, 2, 1, 1, 2, 3, 1, 2, 1, 1, 4, 3, 1, 4, 1, 2, 2, 1, 3, 2, 7, 1, 4, 1, 1, 2, 1, 1, 12, 3, 2, 4, 5, 1, 2, 7, 1, 2, 1, 3, 2, 5, 1, 4, 1, 3, 2, 1, 1, 10, 3, 2, 10, 9, 2, 8, 1, 1, 12, 1, 2, 2, 25, 1, 2, 3, 1, 2, 1, 1, 2, 5, 1, 4, 5, 3, 2, 1, 1, 2, 3, 2, 4, 1, 2, 2, 1, 1, 8, 3, 4, 2, 1, 3, 226, 3, 1, 2
Offset: 1

Views

Author

Keywords

Comments

There exist odd integers 2k-1 such that (2k-1)2^n-1 is always composite.

References

  • Ribenboim, P., The New Book of Prime Number Records. New York: Springer-Verlag, pp. 357-359, 1996.

Crossrefs

Main sequences for Riesel problem: A038699, A040081, A046069, A050412, A052333, A076337, A101036, A108129.
Bisection of A040081.

Programs

  • Mathematica
    max = 10^6; (* this maximum value of m is sufficient up to n=1000 *) a[1] = 2; a[2] = 0; a[n_] := For[m = 1, m <= max, m++, If[PrimeQ[(2*n - 1)*2^m - 1], Return[m]]] /. Null -> -1; Reap[ Do[ Print[ "a(", n, ") = ", a[n]]; Sow[a[n]], {n, 1, 100}]][[2, 1]] (* Jean-François Alcover, Nov 15 2013 *)

A108129 Riesel problem: let k=2n-1; then a(n)=smallest m >= 1 such that k*2^m-1 is prime, or -1 if no such prime exists.

Original entry on oeis.org

2, 1, 2, 1, 1, 2, 3, 1, 2, 1, 1, 4, 3, 1, 4, 1, 2, 2, 1, 3, 2, 7, 1, 4, 1, 1, 2, 1, 1, 12, 3, 2, 4, 5, 1, 2, 7, 1, 2, 1, 3, 2, 5, 1, 4, 1, 3, 2, 1, 1, 10, 3, 2, 10, 9, 2, 8, 1, 1, 12, 1, 2, 2, 25, 1, 2, 3, 1, 2, 1, 1, 2, 5, 1, 4, 5, 3, 2, 1, 1, 2, 3, 2, 4, 1, 2, 2, 1, 1, 8, 3, 4, 2, 1, 3, 226, 3, 1, 2, 1, 1, 2
Offset: 1

Views

Author

Jorge Coveiro, Jun 04 2005

Keywords

Comments

It is conjectured that the integer k = 509203 is the smallest Riesel number, that is, the first n such that a(n) = -1 is 254602.
Browkin & Schinzel, having proved that 509203*2^k - 1 is composite for all k > 0, ask for the first such number with this property, noting that the question is implicit in Aigner 1961. - Charles R Greathouse IV, Jan 12 2018
Record values begin a(1) = 2, a(7) = 3, a(12) = 4, a(22) = 7, a(30) = 12, a(64) = 25, a(96) = 226, a(330) = 800516; the next record appears to be a(1147), unless a(1147) = -1. (The value for a(330), i.e., for k = 659, is from the Ballinger & Keller link, which also lists k = 2293, i.e., n = (k+1)/2 = (2293+1)/2 = 1147, as the smallest of 50 values of k < 509203 for which no prime of the form k*2^m-1 had yet been found.) - Jon E. Schoenfield, Jan 13 2018
Same as A046069 except for a(2) = 1. - Georg Fischer, Nov 03 2018

References

  • Hans Riesel, Några stora primtal, Elementa 39 (1956), pp. 258-260.

Crossrefs

Main sequences for Riesel problem: A038699, A040081, A046069, A050412, A052333, A076337, A101036, A108129.

Programs

  • Mathematica
    Array[Function[k, SelectFirst[Range@300, PrimeQ[k 2^# - 1] &]][2 # - 1] &, 102] (* Michael De Vlieger, Jan 12 2018 *)
    smk[n_]:=Module[{m=1,k=2n-1},While[!PrimeQ[k 2^m-1],m++];m]; Array[smk,120] (* Harvey P. Dale, Dec 26 2023 *)
  • PARI
    forstep(k=1,301,2,n=1;while(!isprime(k*2^n-1),n++);print1(n,","))

Extensions

Edited by Herman Jamke (hermanjamke(AT)fastmail.fm), Oct 25 2006
Name corrected by T. D. Noe, Feb 13 2011

A374965 a(n) = 2*a(n-1) + 1 for a(n-1) not prime, otherwise a(n) = prime(n) - 1; with a(1)=1.

Original entry on oeis.org

1, 3, 4, 9, 19, 12, 25, 51, 103, 28, 57, 115, 231, 463, 46, 93, 187, 375, 751, 70, 141, 283, 82, 165, 331, 100, 201, 403, 807, 1615, 3231, 6463, 12927, 25855, 51711, 103423, 156, 313, 166, 333, 667, 1335, 2671, 192, 385, 771, 1543, 222, 445, 891, 1783, 238, 477
Offset: 1

Views

Author

Bill McEachen, Jul 25 2024

Keywords

Comments

Sequence is clearly infinite and not monotonic. Primes are sparse.
When is the next prime after n=10016 ? [Answer from N. J. A. Sloane, Aug 01 2024: The point of Bill's question is that a(10016) is the prime 838951, which is in fact the 289th prime in this sequence, as can be seen from A375028 and A373799. Thanks to the work of Lucas A. Brown (see A050412), we now know that the answer to Bill's question is that the 290th prime is the 102410-digit prime 104917*2^340181 - 1 = 5079...8783, which is a(350198). It was a very good question!]
It appears that the trajectories for different initial conditions a(1) converge to a few attractors. For all prime values and most nonprime values of a(1), the trajectories converge to the same attractor with prime 838951 at n=10016. For a(1) = 147, 295, 591, 1183, ... the trajectories converge to prime 85796863 at n=4390. For a(1) = 658, the trajectory reaches a prime with 240983 digits after 800516 steps. For a(1) = 509202, the trajectory never reaches a prime (see A050412, A052333). - Chai Wah Wu, Jul 29 2024

Examples

			a(1) = 1 is not a prime, so a(2) = 2*1+1 = 3. a(2) is a prime, so a(3) = prime(3)-1 = 4. a(4) = 2*4+1 = 9.
		

Crossrefs

The primes are listed in A375028 (see also A373798 and A373804).
Cf. A050412 and A052333.

Programs

  • Mathematica
    a[n_] := a[n] = If[!PrimeQ[a[n-1]], 2*a[n-1] + 1, Prime[n]-1]; a[1] = 1; Array[a, 60] (* Amiram Eldar, Jul 25 2024 *)
    nxt[{n_,a_}]:={n+1,If[!PrimeQ[a],2a+1,Prime[n+1]-1]}; NestList[nxt,{1,1},60][[;;,2]] (* Harvey P. Dale, Jul 28 2024 *)
  • Python
    from itertools import islice
    from sympy import isprime, nextprime
    def A374965_gen(): # generator of terms
        a, p = 1, 3
        while True:
            yield a
            a, p = p-1 if isprime(a) else (a<<1)+1, nextprime(p)
    A374965_list = list(islice(A374965_gen(),30)) # Chai Wah Wu, Jul 29 2024
Showing 1-10 of 28 results. Next