A351728 Primes p such that if q is the next prime, p+A004086(q) and q+A004086(p) are prime.

2, 61, 83, 433, 677, 2351, 2399, 2441, 4397, 4457, 4673, 6257, 6367, 6961, 8263, 8713, 8761, 20627, 21391, 21649, 22721, 22871, 23227, 23761, 25111, 25321, 25589, 25609, 25741, 25841, 26597, 26731, 26981, 27179, 27271, 27299, 27367, 27409, 27481, 27961, 28559, 29881, 40609, 40927, 40933, 42197

Offset: 1

Author

J. M. Bergot and Robert Israel, Mar 20 2022

Keywords

• nonn
• base

Comments

For each term p except 2, A013632(p) is divisible by 6.

Examples

			a(3) = 83 is a term because it is prime, the next prime is 89, and 83+98 = 181 and 38+89 = 127 are both prime.
		

Links

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

Crossrefs

Cf. A004086, A013632.

Programs

Maple

revdigs:= proc(n) local L,i; L:= convert(n,base,10);
  add(L[-i]*10^(i-1),i=1..nops(L))
end proc:
Primes:= select(isprime, [2,seq(i,i=3..10000,2)]):
RPrimes:= map(revdigs,Primes):
Primes[select(i -> isprime(Primes[i]+RPrimes[i+1]) and isprime(RPrimes[i]+Primes[i+1]), [$1..nops(Primes)-1])]: