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.

A051653 Primes p such that 2310*p + 1 is also prime.

Original entry on oeis.org

2, 5, 11, 19, 47, 71, 73, 97, 109, 113, 137, 139, 151, 163, 167, 199, 229, 239, 263, 271, 311, 313, 317, 331, 347, 353, 379, 431, 433, 457, 461, 479, 503, 509, 523, 593, 599, 613, 617, 691, 701, 727, 761, 797, 811, 823, 853, 863, 883, 929, 937, 941, 947, 953
Offset: 1

Views

Author

Keywords

Examples

			p = 2 and 2310*p + 1 = 4621 are both primes.
		

Crossrefs

Programs

  • Mathematica
    Select[Prime[Range[200]],PrimeQ[2310#+1]&] (* Harvey P. Dale, Oct 14 2014 *)
  • PARI
    isok(k) = isprime(k) && isprime(2310*k+1); \\ Amiram Eldar, Feb 24 2025

Formula

a(n) = (A051649(n)-1)/2310. - Amiram Eldar, Feb 24 2025

A051902 Minimal primorial safe primes: p and primorial*p + 1 are both primes.

Original entry on oeis.org

5, 13, 61, 421, 4621, 150151, 8678671, 106696591, 2454021571, 71166625531, 401120980261, 170676977100631, 2129751844690471, 562558737261811291, 11682905869181336791, 97767475431570134191, 9613801750771063195351, 234576762718813941966541, 55008250857561869391153631
Offset: 1

Views

Author

Labos Elemer, Dec 16 1999

Keywords

Comments

In A051888, 13 of the first 25 values are distinct, while here all corresponding min-primorial-safe-primes are different: {2,5,17,11,23,43,19,3,7,61,53,41,47}.

Examples

			The first five terms of A051887 are 2, so the first 5 terms have the form 1 + 2*A002110(n): 5, 13, 61, 421, 4621, which are smallest terms in A005385, A051644, A051646, A051648, A051649. The 6th term here is A051651(1) = A051887(6)*A002110(6) + 1 = 5*30030 + 1.
		

Crossrefs

Programs

  • Mathematica
    a[n_] := Module[{p = 2, r =Times @@ Prime[Range[n]]}, While[!PrimeQ[p * r + 1], p = NextPrime[p]]; p * r + 1]; Array[a, 20] (* Amiram Eldar, Feb 25 2025 *)
  • PARI
    a(n) = {my(p = 2, r = vecprod(primes(n))); while(!isprime(p * r + 1), p = nextprime(p+1)); p * r + 1;} \\ Amiram Eldar, Feb 25 2025

Formula

a(n) = 1 + A002110(n)*A051887(n).
Showing 1-2 of 2 results.