fork download
  1. # prime fibonacci numbers
  2.  
  3. def isPrime(n, k=5): # miller-rabin
  4. from random import randint
  5. if n < 2: return False
  6. for p in [2,3,5,7,11,13,17,19,23,29]:
  7. if n % p == 0: return n == p
  8. s, d = 0, n-1
  9. while d % 2 == 0:
  10. s, d = s+1, d/2
  11. for i in range(k):
  12. x = pow(randint(2, n-1), d, n)
  13. if x == 1 or x == n-1: continue
  14. for r in range(1, s):
  15. x = (x * x) % n
  16. if x == 1: return False
  17. if x == n-1: break
  18. else: return False
  19. return True
  20.  
  21. a, b, f = 1, 1, 2
  22. while f < 1000000000:
  23. if isPrime(f): print f
  24. a, b, f = b, f, b+f
Success #stdin #stdout 0.01s 10080KB
stdin
Standard input is empty
stdout
2
3
5
13
89
233
1597
28657
514229
433494437