OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..336 (terms below 10^11)
EXAMPLE
The prime 3631 is in the sequence because reversed it is 1363 = 29*47 and abs(29*phi-47) = 0.077014... < 1.
MATHEMATICA
revDigits[n_] := FromDigits @ Reverse @ IntegerDigits[n]; goldQ[n_] := Module[{f = FactorInteger[n]}, If[Length[f] != 2, False, If[Max[f[[;; , 2]]] != 1, False, Abs[f[[2, 1]] - f[[1, 1]] * GoldenRatio] < 1]]]; seqQ[n_] := PrimeQ[n] && goldQ @ revDigits[n]; Select[Range[4*10^6], seqQ] (* Amiram Eldar, Nov 29 2019 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Jason Earls, Jun 20 2005
EXTENSIONS
a(19)-a(25) from Donovan Johnson, Nov 11 2008
Offset corrected by Amiram Eldar, Nov 29 2019
STATUS
approved