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.

Previous Showing 71-72 of 72 results.

A307486 a(0) = 3; a(n) = smallest k > 1 such that 1 + a(0)*a(1)*...*a(n-1)*k is composite.

Original entry on oeis.org

3, 3, 3, 2, 4, 3, 3, 3, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 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, 2, 2, 2, 2, 2, 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, 2
Offset: 0

Views

Author

Thomas Ordowski, Apr 10 2019

Keywords

Comments

a(n) = 4 for n = 4 and 39; a(n) = 3 for n = 0, 1, 2, 5, 6, 7, 14, 20, 25, 56, 90, 119, 316, 330, 1268, 1604, 1805, 1880, 1984, 2950, 3386, 3712, 4532, 4874, 8968, 18178, 19454, 23272, 45617, 51980, 52780, 60552, ...; a(n) = 2 for others n < 60552. Data from Amiram Eldar.
Similar tails have sequences with other initial terms that are natural numbers.
Conjecture: for any initial term a(0) > 0, a(n) > 3 only for finitely many n >= 0.
The question is how to prove that all these sequences are bounded, so bounded?
It seems that a(0) = 21 is the smallest initial term such that a(n) = 2 or 3 for every n > 0.
Note that if a(0) is a Sierpinski number A076336, then a(n) = 2 for every n > 0.
Carl Pomerance (in a letter to the author) wrote: Since the sequence a(0)a(1)...a(n-1) grows geometrically, the chance that any given k will give a prime is about 1/n, so the chance that both k = 2 and k = 3 will give primes is about 1/n^2, heuristically, which has a convergent sum. Thus, for this reason I would agree it's reasonable to conjecture that for any starting a(0), the sequence will eventually be made up of 2's and 3's, with 2's predominating. The number of 3's among the first n terms should be proportional to log n.

Crossrefs

Programs

  • Mathematica
    a[0] = 3; a[n_] := a[n] = Module[{k = 2, p = Product[a[i], {i, 0, n - 1}]}, While[PrimeQ[1 + p*k], k++]; k]; Array[a, 100, 0] (* Amiram Eldar, Apr 10 2019 *)
  • PARI
    p=1; for (n=0, 100, for (k=2, oo, if (!isprime(1+p*k), print1 (k", "); p*=k; break))) \\ Rémy Sigrist, Apr 23 2019

Extensions

More terms from Amiram Eldar, Apr 10 2019

A368560 a(1) = 78557 (the first Sierpinski number); thereafter a(n+1) = Od(3*5*7*13*19*37*73 - a(n)), where Od(m) is the odd part of m.

Original entry on oeis.org

78557, 34985939, 2191531, 8482363, 7696009, 31177213, 19436611, 790841, 34629797, 17710319, 13085029, 28482703, 10391933, 29829251, 78557, 34985939, 2191531, 8482363, 7696009, 31177213, 19436611, 790841, 34629797, 17710319, 13085029, 28482703, 10391933, 29829251
Offset: 1

Views

Author

Thomas Ordowski, Dec 30 2023

Keywords

Comments

Generally, if k is a Sierpinski number (or is a Riesel number) and P(k) > k is the product of all elements from the covering set for k*2^n + 1 (or for k*2^n - 1), then Od(P(k) - k) is a Riesel number (or is a Sierpinski number) with the same covering set, where Od(m) is the odd part of m.
Thus a(2n-1) is a Sierpinski number and a(2n) is a Riesel number.
This sequence is purely periodic with period P = 14.

Examples

			a(1) = 78557 is a Sierpinski number and a(2) = (3*5*7*13*19*37*73 - 78557)/2 = 34985939 is a Riesel number with the same covering set {3, 5, 7, 13, 19, 37, 73}.
		

Crossrefs

Programs

  • Mathematica
    od[n_] := n/2^IntegerExponent[n, 2]; a[1] = 78557; a[n_] := a[n] = od[70050435 - a[n-1]]; Array[a, 42] (* Amiram Eldar, Dec 30 2023 *)

Formula

a(n + 14) = a(n).
Previous Showing 71-72 of 72 results.