OFFSET
1,1
PROG
(PARI) print1(2, ", "); anm2=1; anm1=3; for(n=3, 1000, an=3*anm1+8*anm2; anm2=anm1; anm1=an; if(isprime(an), print1(n, ", ")))
(Perl) use ntheory ":all"; sub isok { is_prob_prime(lucasu(3, -8, $_[0])) }; for (1..700) { print "$_\n" if isok($_) }; # Daniel Suteu, Jan 12 2019
CROSSREFS
KEYWORD
nonn,more,hard
AUTHOR
Hugo Pfoertner, Jan 12 2019
EXTENSIONS
a(13) from Daniel Suteu, Jan 12 2019
STATUS
approved