login
A356147
Lower twin primes whose sum of digits is a lower twin prime.
1
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
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
Subsequence of A046704.
Sequence in context: A141578 A327748 A272345 * A364762 A346659 A067200
KEYWORD
nonn,base
AUTHOR
J. M. Bergot and Robert Israel, Jul 27 2022
STATUS
approved