|
|
A155011
|
|
Sophie Germain Fibonacci prime numbers.
|
|
3
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
2*2+1=5, 2*3+1=7, 2*5+1=11, ...
No additional terms up through Fibonacci(10000). - Harvey P. Dale, Nov 26 2013
No additional terms up through Fibonacci(50000). - Chai Wah Wu, Nov 04 2015
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
a={}; Do[f=Fibonacci[n]; If[PrimeQ[f], If[PrimeQ[2*f+1], AppendTo[a, f]]], {n, 3*6!}]; a
Select[Fibonacci[Range[2000]], And@@PrimeQ[{#, 2#+1}]&] (* Harvey P. Dale, Nov 26 2013 *)
|
|
PROG
|
(Python)
from gmpy2 import is_prime
a, b, a2, b2 = 0, 1, 1, 3
for _ in range(10**6):
if is_prime(b) and is_prime(b2):
a, b, a2, b2 = b, a+b, b2, a2+b2-1 # Chai Wah Wu, Nov 04 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|