Advertisement
Guest User

Untitled

a guest
Feb 26th, 2018
233
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.32 KB | None | 0 0
  1. from decimal import Decimal
  2. import math
  3. import numpy as np
  4. def fibgen(n):
  5. y = math.sqrt(5)
  6. x = Decimal((1+Decimal(y))**n-(1-Decimal(y))**n)//Decimal(2**n*Decimal(y))
  7. print(x)
  8. return x
  9.  
  10. n=5000
  11. while True:
  12. if(int(math.log10(fibgen(n))+1) == 1000):
  13. print(n)
  14. break
  15. n += 1
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement