A094069 Prime(p)-4 for primes p such that prime(p) - 4 is prime.

7, 13, 37, 79, 349, 397, 457, 613, 769, 1213, 1429, 1783, 2347, 2377, 2473, 3907, 4129, 4513, 4783, 5437, 5647, 7477, 7603, 7879, 8803, 9829, 10429, 10453, 10627, 11443, 11863, 11923, 12109, 13033, 13099, 13327, 14173, 14779, 15679, 16057, 16069

Offset

1,1

Comments

Primes q such that q+4 is in A006450. - Robert Israel, Nov 09 2020

Links

Robert Israel, Table of n, a(n) for n = 1..10000

Example

Prime(13) = 41. 41 - 4 = 37.

Maple

R:= NULL: p:= 7:
for i from 5 to 10000 do
  q:= p; p:= nextprime(p);
  if isprime(i) and q=p-4 then R:= R, q fi;
od:
R; # Robert Israel, Nov 09 2020

Mathematica

Prime[#]&/@Select[Prime[Range[2, 300]], PrimeQ[Prime[#]-4]&]-4 (* Harvey P. Dale, May 21 2021 *)

Prog

(PARI) a(n) = { forprime(x=2, n, y=prime(x)- 4; if(isprime(y), print1(y", ")) ) }

Crossrefs

Cf. A006450.

This is a subsequence of A172367.

Sequence in context: A118819 A118525 A213537 * A338812 A052378 A090607

Adjacent sequences: A094066 A094067 A094068 * A094070 A094071 A094072

Keyword

nonn

Author

Cino Hilliard, May 31 2004

Extensions

Definition corrected by Robert Israel, Nov 09 2020

Status

approved