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

A321058 Indices of records in A067849.

Original entry on oeis.org

1, 2, 44, 561329, 9549959, 42932384, 13044904289, 277344139214, 2045466215756534, 47702521115271164
Offset: 1

Views

Author

Torlach Rush, Oct 26 2018

Keywords

Comments

The corresponding record values are 2,4,6,7,8,9,10,12,13,14.
From David A. Corneth, Nov 10 2018: (Start)
Terms a(n) are of the form 3*k+2 for n > 1.
If 2^k - 1 is composite then a(n) is not divisible by any prime factor of 2^k-1 for n > k. So for example, gcd(a(n), 105) = 1 for n > 5. (End)
From Glen Whitney, Sep 14 2022: (Start)
Similarly to Corneth's observations, modulo any prime p, any residue for a(n) of the form 2^k - 1 mod p is forbidden for n greater than or equal to the number of such residues; for example a(n) may not be congruent to 0, 1, or 3 mod 7 for n >= 3.
For n > 2, if a(n) appears in this sequence, 2a(n) + 1 must appear in A057331. (End)

Examples

			2 is a term because A067849(2) = 4 > A067849(1) = 2.
44 is a term because A067849(44) = 6 > A067869(2) = 4.
		

Crossrefs

Programs

  • PARI
    b(n) = {my(nb = 0, newn); while (isprime(newn=2*n+1), nb++; n = newn); nb;} \\ A067849
    lista(nn) = {my(mmax = -1, mm); for (n=1, nn, if ((mm=b(n)) > mmax, mmax = mm; print1(n, ", ")););} \\ Michel Marcus, Nov 10 2018

Extensions

a(7) from Amiram Eldar, Nov 10 2018
a(8)-a(10) from A057331 by Glen Whitney, Sep 14 2022

A057331 a(n) = smallest prime p such that the first n iterates of p under x->2x+1 are all primes.

Original entry on oeis.org

2, 2, 2, 2, 2, 89, 1122659, 19099919, 85864769, 26089808579, 554688278429, 554688278429, 4090932431513069, 95405042230542329
Offset: 0

Views

Author

Patrick De Geest, Aug 15 2000

Keywords

Comments

For n > 10, a(n) == -1 (mod 2*3*5*11*13). - Farideh Firoozbakht, Apr 24 2004
From Glen Whitney, Sep 14 2022: (Start)
Extending Firoozbakht's observation, modulo any prime p, all residues of a(n) of the form 2^k - 1 mod p are forbidden for n greater than or equal to the number of such residues, e.g., a(n) may not be congruent to 1 or 3 mod 7 for n >= 2.
A067849(a(n)) >= n and for each odd a(n) that occurs in this sequence, (a(n)-1)/2 occurs in A321058. (End)

Examples

			a(5) = 89 because the numbers 89, 179, 359, 719, 1439, 2879 are all primes and 89 is the first number to have this property.
		

Crossrefs

See A067849 (number of prime iterates starting from any n) and A321058 (starting points that yield record numbers of iterates).
See also A005602.

Programs

  • Mathematica
    f[n_] := 2n + 1; k = 1; Do[ While[ Union[ PrimeQ[ NestList[ f, Prime[k], n]]] != {True}, k++ ]; Print[ Prime[k]], {n, 1, 9} ]
  • PARI
    has(p,n)=for(k=1,n,if(!isprime(p), return(0)); p=2*p+1); isprime(p)
    a(n)=forprime(p=2,, if(has(p,n), return(p))) \\ Charles R Greathouse IV, Apr 29 2015

Extensions

More terms from Farideh Firoozbakht, Apr 24 2004
a(11) (from the Caldwell link) sent by Peter Deleu, Hulste, Belgium, Nov 22 2004
a(13) added from A005602, Paul Zimmermann, Mar 09 2018
Showing 1-2 of 2 results.