OFFSET
1,1
COMMENTS
Primes not in A133247.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(p) local a, b, i;
if not isprime(p) then return false fi;
a:= 0: b:= 1;
for i from 2 do
a, b:= b, (a+b) mod p;
if b = 0 then
if i mod 3 <> 0 then return true
elif a = 1 then return false
fi
fi
od:
end proc:
select(filter, [seq(i, i=3..1000, 2)]); # Robert Israel, Nov 20 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Tanya Khovanova, Jan 30 2009
STATUS
approved