OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..2595 (terms below 10^11)
EXAMPLE
1122 is a term because 1122/6=187 and 187=11*17 and 11*phi-17 = 0.798373... < 1.
MATHEMATICA
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]]]; sumDigits[n_] := Plus @@ IntegerDigits[n]; seqQ[n_] := Divisible[n, (sd = sumDigits[n])] && goldQ[n/sd]; Select[Range[360000], seqQ] (* Amiram Eldar, Nov 29 2019 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Jason Earls, Jun 26 2005
STATUS
approved