cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A083668 Prime indices of prime Fibonacci numbers.

This page as a plain text file.
%I A083668 #33 Oct 21 2016 11:06:53
%S A083668 3,5,7,11,13,17,23,29,43,47,83,131,137,359,431,433,449,509,569,571,
%T A083668 2971,4723,5387,9311,9677,14431,25561,30757,35999,37511,50833,81839,
%U A083668 104911,130021,148091,201107,397379,433781,590041,593689,604711,931517,1049897,1285607,1636007,1803059,1968721
%N A083668 Prime indices of prime Fibonacci numbers.
%C A083668 Same as A001605 without the number 4.
%C A083668 From _V. Raman_, Oct 04 2012: (Start)
%C A083668 Also the indices of prime Fibonacci numbers which can be written as the sum of two positive squares.
%C A083668 The Fibonacci numbers F(6k+1) and F(6k+5) are congruent to 1 (mod 4).
%C A083668 (End)
%e A083668 For Fib(n) to be prime, n must be prime, except for n=4. The first 9 primes are: 2, 3, 5, 7, 11, 13, 17, 19 and 23. The corresponding Fibonacci numbers are: 1, 2, 5, 13, 89, 233, 1597, 4181 and 28657. All of these are prime except Fib(2) = 1 and Fib(19) = 4181. So the first 7 terms of this sequence are 3, 5, 7, 11, 13, 17 and 23.
%t A083668 Do[ If[ PrimeQ[ Fibonacci[ Prime[n]]], Print[ Prime[n]]], {n, 1, 1000}]
%o A083668 (PARI) pif(n) = { forprime(x=2,n, if(isprime(fibonacci(x)), print1(x" "))) }
%o A083668 (PARI) is(p)=isprime(p) & ispseudoprime(fibonacci(p))  \\ _Charles R Greathouse IV_, Sep 19 2012
%Y A083668 Cf. A001605, A075737.
%K A083668 nonn
%O A083668 1,1
%A A083668 _Cino Hilliard_, Jun 14 2003
%E A083668 More terms from _Zak Seidov_, Aug 31 2006
%E A083668 Replaced the erroneous example _Harry J. Smith_, Jan 16 2009
%E A083668 Terms a(42) to a(47) added by _V. Raman_, Oct 04 2012
%E A083668 Definition and wrong statement in example corrected by _M. F. Hasler_, Oct 08 2012