A356147 Lower twin primes whose sum of digits is a lower twin prime.

3, 5, 29, 41, 137, 179, 191, 197, 227, 269, 281, 311, 461, 641, 809, 821, 827, 881, 1019, 1031, 1091, 1277, 1301, 1451, 1619, 1721, 1871, 2027, 2081, 2087, 2111, 2267, 2339, 2591, 2711, 2801, 2999, 3167, 3251, 3257, 3329, 3527, 3581, 3671, 3851, 4001, 4049, 4157, 4229, 4241, 4337, 4421, 4481

Offset

1,1

Comments

Members k of A001359 such that A007953(k) is in A001359.

Links

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

Example

a(3) = 29 is a term because 29, 29 + 2 = 31, 2 + 9 = 11 and 2 + 9 + 2 = 13 are primes.

Maple

filter:= proc(n) local s;
  if not(isprime(n) and isprime(n+2)) then return false fi;
   s:= convert(convert(n, base, 10), `+`);
   isprime(s) and isprime(s+2)
end:
select(filter, [3, seq(i, i=5..5000, 6)]);

Mathematica

Select[Prime[Range[600]], AllTrue[{# + 2, (d = Plus @@ IntegerDigits[#]), d + 2}, PrimeQ] &] (* Amiram Eldar, Jul 28 2022 *)

Crossrefs

Cf. A001359, A007953.

Subsequence of A046704.

Sequence in context: A141578 A327748 A272345 * A364762 A346659 A067200

Adjacent sequences: A356144 A356145 A356146 * A356148 A356149 A356150

Keyword

nonn,base

Author

J. M. Bergot and Robert Israel, Jul 27 2022

Status

approved