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