A155916 - OEIS

%I #7 Nov 21 2016 04:18:53

%S 3,5,7,11,13,29,37,41,43,47,59,67,71,73,89,97,101,103,113,127,131,139,

%T 149,151,157,163,179,191,193,199,223,233,239,251,263,269,277,281,283,

%U 307,311,313,331,337,347,353,359,367,373,389,397,401,419,431,433,449

%N Prime factors of odd Fibonacci numbers.

%C Primes not in A133247.

%H Robert Israel, <a href="/A155916/b155916.txt">Table of n, a(n) for n = 1..10000</a>

%p filter:= proc(p) local a,b,i;

%p if not isprime(p) then return false fi;

%p a:= 0: b:= 1;

%p for i from 2 do

%p   a,b:= b, (a+b) mod p;

%p   if b = 0 then

%p     if i mod 3 <> 0 then return true

%p     elif a = 1 then return false

%p     fi

%p   fi

%p od:

%p end proc:

%p select(filter, [seq(i,i=3..1000,2)]); # _Robert Israel_, Nov 20 2016

%Y Cf. A014437, A133247.

%K nonn

%O 1,1

%A _Tanya Khovanova_, Jan 30 2009