A271981 Primes p such that p + 40 is also prime.

3, 7, 13, 19, 31, 43, 61, 67, 73, 97, 109, 127, 139, 151, 157, 193, 199, 211, 223, 229, 241, 271, 277, 307, 313, 349, 379, 409, 421, 439, 463, 523, 547, 577, 601, 607, 613, 619, 643, 661, 733, 757, 769, 787, 823, 907, 937, 991, 1009, 1021, 1051, 1063, 1069

Offset: 1

Author

Karl V. Keller, Jr., Apr 17 2016

Keywords

• nonn

Comments

A126721 is a subsequence of this sequence.

Examples

			3 is a term since 3 + 40 = 43 is also prime.
7 is a term since 7 + 40 = 47 is also prime.
		

Links

• Karl V. Keller, Jr., Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A000040, A126721.

Programs

Maple

q:= n-> andmap(isprime, [n, n+40]):
select(q, [$2..2000])[];  # Alois P. Heinz, Jul 21 2022

Mathematica

Select[Prime@ Range@ 180, PrimeQ[# + 40] &] (* Michael De Vlieger, Apr 18 2016 *)

PARI

lista(nn) = forprime(p=2, nn, if (isprime(p+40), print1(p, ", "))); \\ Michel Marcus, Apr 19 2016

Python

from sympy import isprime
for i in range(3, 2001,2):
     if isprime(i) and isprime(i+40): print (i,end=', ')