Skip to content

Instantly share code, notes, and snippets.

@Hermann-SW

Hermann-SW/Fibonacci_prime_indices.gp

Created August 14, 2024 14:01
Show Gist options
  • Save Hermann-SW/8780e5ce7deae505022803a108fce8fc to your computer and use it in GitHub Desktop.
Save Hermann-SW/8780e5ce7deae505022803a108fce8fc to your computer and use it in GitHub Desktop.
Download ZIP
Fibonacci number prime indices up to largest known sofar fibonacci prime index 201107
Raw
Fibonacci_prime_indices.gp
Fibonacci_prime_indices={[
3, 4, 5, 7, 11, 13, 17, 23, 29, 43,
47, 83, 131, 137, 359, 431, 433, 449, 509, 569,
571, 2971, 4723, 5387, 9311, 9677, 14431, 25561, 30757, 35999,
37511, 50833, 81839, 104911, 130021, 148091, 201107
]}
@Hermann-SW
Copy link
Author

Hermann-SW commented Aug 14, 2024 •
edited
Loading

https://oeis.org/A001605

pi@raspberrypi5:~ $ gp -q Fibonacci_prime_indices.gp 
? #Fibonacci_prime_indices
37
? Fibonacci_prime_indices[37]
201107
? foreach(Fibonacci_prime_indices,i,if(i<1000,print1(isprime(fibonacci(i)))))
111111111111111111111
? ##
  ***   last result: cpu time 638 ms, real time 213 ms.
? vecmax([x|x<-Fibonacci_prime_indices,x<1000])
571
? #digits(fibonacci(571))
119
? #digits(fibonacci(201107))
42029
?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment