OFFSET
0,3
FORMULA
a(n) = if Hofstadter[n] is prime then Hofstadter[n] else zero
MATHEMATICA
digits = 200 Hofstadter[n_Integer? Positive] := Hofstadter[n] = Hofstadter[n - Hofstadter[n - 1]] + Hofstadter[n - Hofstadter[n - 2]] Hofstadter[0] = Hofstadter[1] = 1 b = Table[If[PrimeQ[ Hofstadter[n]] == True, Hofstadter[n], 0], {n, 0, digits - 1}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Mar 10 2005
STATUS
approved