OFFSET
1,1
COMMENTS
All terms greater than a(11)=25 are prime so far.
a(80) > 10^6. - Paul Vanderveen, Oct 02 2022
LINKS
Paul Vanderveen, Table of n, a(n) for n = 1..79
FORMULA
MATHEMATICA
a[0] = 0; a[1] = a[2] = a[3] = a[5] = 1; a[4] = 5; a[n_] := a[n] = -a[n - 1] + a[n - 2] + 3a[n - 3] + a[n - 4] - a[n - 5] - a[n - 6]; lst = {}; Do[ If[ PrimeQ[a [ n]], AppendTo[lst, n]], {n, 21900}]; lst
PROG
(Haskell)
import Data.List (findIndices)
a104499 n = a104499_list !! (n-1)
a104499_list = findIndices ((== 1) . a010051) a001945_list
-- Reinhard Zumkeller, Jan 11 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Mar 10 2005
EXTENSIONS
a(50) and a(52) corrected and more terms from Paul Vanderveen, Jul 27 2022
STATUS
approved