OFFSET
1,1
COMMENTS
MATHEMATICA
s[1] = s[2] = s[3] = 1; s[n_] := s[n] = s[n - 1] + s[n - 2] - s[n - 3] + GCD[s[n - 1], s[n - 3]]; Select[Array[s, 180], PrimeQ] (* Amiram Eldar, Feb 24 2023 *)
CROSSREFS
KEYWORD
nonn,fini
AUTHOR
Jack Braxton, Feb 12 2023
STATUS
approved