OFFSET
1,1
COMMENTS
Primes p such that both p and p+2 are happy primes.
MATHEMATICA
happyQ[n_] := NestWhile[Total[IntegerDigits[#]^2] &, n, UnsameQ, All] == 1; Select[Prime[Range[30000]], PrimeQ[# + 2] && happyQ[#] && happyQ[# + 2] &] (* Amiram Eldar, Aug 23 2025 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Anonymous, Aug 23 2025
STATUS
approved
