OFFSET
1,2
COMMENTS
These numbers are the positions of prime numbers in A066629.
MAPLE
with(combinat): A271729:=n->`if`(isprime(2*fibonacci(n+2)+((-1)^n-3)/2), n, NULL): seq(A271729(n), n=1..3*10^3); # Wesley Ivan Hurt, Apr 13 2016
MATHEMATICA
Select[Range[10000], PrimeQ[(2 Fibonacci[# + 2] + ((-1)^# - 3) / 2)] &]
PROG
(Magma) [n: n in [0..2000] | IsPrime(2*Fibonacci(n+2)+((-1)^n-3) div 2)];
(PARI) lista(nn) = for(n=1, nn, if(ispseudoprime(2*fibonacci(n+2)+((-1)^n-3)/2), print1(n, ", "))); \\ Altug Alkan, Apr 13 2016
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Vincenzo Librandi, Apr 13 2016
STATUS
approved