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.

A053409 Fibonacci numbers which are semiprimes.

This page as a plain text file.
%I A053409 #32 Jun 26 2020 19:01:45
%S A053409 21,34,55,377,4181,17711,121393,1346269,5702887,165580141,53316291173,
%T A053409 956722026041,2504730781961,308061521170129,806515533049393,
%U A053409 14472334024676221,1779979416004714189,19740274219868223167,573147844013817084101,10284720757613717413913
%N A053409 Fibonacci numbers which are semiprimes.
%C A053409 Subsequence of A006881, since the only square Fibonacci numbers are 1 and 144. - _Charles R Greathouse IV_, Sep 24 2012
%C A053409 Apart from a(1) = 21, all terms are of the form F(p), F(2p), or F(p^2) where F(n) is the n-th Fibonacci number. - _Charles R Greathouse IV_, Oct 06 2016
%H A053409 Vincenzo Librandi, <a href="/A053409/b053409.txt">Table of n, a(n) for n = 1..48</a>
%F A053409 a(n) = A000045(A072381(n)).
%t A053409 Select[Fibonacci@Range[120],Last/@FactorInteger[#]=={1,1}&] (* _Vladimir Joseph Stephan Orlovsky_, Jan 29 2012 *)
%t A053409 Select[Fibonacci[Range[150]],PrimeOmega[#]==2&] (* _Harvey P. Dale_, Jun 26 2020 *)
%o A053409 (PARI) issemi(n)=bigomega(n)==2
%o A053409 list(lim)=my(v=List([21]),F,t); forprime(p=2,, F=fibonacci(p); if(F>lim, break); if(issemi(F), listput(v,F))); forprime(p=2,, F=fibonacci(p^2); if(F>lim, break); if(isprime(t=fibonacci(p)) && isprime(F/t), listput(v,F))); forprime(p=2,, F=fibonacci(2*p); if(F>lim, break); if(isprime(t=fibonacci(p)) && isprime(F/t), listput(v,F))); Set(v) \\ _Charles R Greathouse IV_, Oct 06 2016
%Y A053409 Cf. A000045, A001358, A006881, A072381.
%Y A053409 Column k=2 of A303216.
%K A053409 nonn
%O A053409 1,1
%A A053409 _G. L. Honaker, Jr._, Jan 09 2000