OFFSET
1,2
MATHEMATICA
Select[Range[10000], PrimeQ[2 Fibonacci[#] + 1]&] (* Vincenzo Librandi, Aug 13 2013 *)
PROG
(PARI) isok(n) = isprime(2*fibonacci(n)+1) \\ Michel Marcus, Jun 03 2013
(Magma) [n: n in [1..1200] | IsPrime(2*Fibonacci(n)+1)]; // Vincenzo Librandi, Aug 13 2013
(PARI) for(n=1, 10^9, if(ispseudoprime(2*fibonacci(n)+1), print1(n, ", "))); \\ Joerg Arndt, Aug 13 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Teofilatto, Dec 12 2006
EXTENSIONS
Inserted a(1)=1 and extended by Michel Marcus, Jun 03 2013
More terms from Vincenzo Librandi, Aug 13 2013
a(48) from Jorge Coveiro, Sep 05 2022
a(49)-a(55) from Michael S. Branicky, Jun 17 2023
STATUS
approved