|
|
A318294
|
|
Prime-indexed primes p such that 2*p + 1 is also a prime-indexed prime.
|
|
1
|
|
|
5, 41, 1031, 3299, 3761, 6323, 7649, 7841, 10589, 12959, 18443, 37853, 38039, 44729, 58511, 59723, 60509, 77813, 105971, 134489, 137999, 190523, 194933, 200351, 239633, 250751, 266093, 288653, 315011, 366161, 376241, 377633, 393401, 410999, 420419, 420779, 450431, 474983
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Prime-indexed Sophie Germain primes whose associated safe primes are also prime-indexed. Subsequence of A005384.
|
|
LINKS
|
|
|
EXAMPLE
|
5 = prime(3) is a prime-indexed prime, and 2*5 + 1 = 11 = prime(5) is also a prime-indexed prime. Therefore 5 is in the sequence.
3 = prime(2) is a prime-indexed prime, but 2*3 + 1 = 7 = prime(4) is not a prime-indexed prime since 4 is composite. Therefore 3 is not in the sequence.
|
|
MATHEMATICA
|
Select[Table[Prime[p], {p, Prime[Range[5000]]}], AllTrue[{2#+1, PrimePi[2#+1]}, PrimeQ]&] (* Harvey P. Dale, Jun 18 2023 *)
|
|
PROG
|
(PARI) isok(p) = isprime(p) && isprime(primepi(p)) && isprime(2*p+1) && isprime(primepi(2*p+1)); \\ Michel Marcus, Aug 28 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|