Guest User

Untitled

a guest
Dec 9th, 2018
57
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.75 KB | None | 0 0
  1. def Fib(n):
  2.   F=[[1,1],
  3.      [1,0]]
  4.   if (n==0):
  5.     return
  6.   power(F,n-1)
  7.   return F[0][0]
  8.  
  9. def multiply(F,M):
  10.  
  11.   w = (F[0][0] * M[0][0])+(F[0][1] * M[1][0])
  12.   x = (F[0][0] * M[0][1])+(F[0][1] * M[1][1])
  13.   y = (F[1][0] * M[0][0])+(F[1][1] * M[1][0])
  14.   z = (F[1][0] * M[0][1])+(F[1][0] * M[1][1])
  15.    
  16.   F[0][0] = w
  17.   F[0][1] = x
  18.   F[1][0] = y
  19.   F[1][1] = z
  20.  
  21. def power(F,n):
  22.   M = [[1,1],
  23.        [1,0]]
  24.   for i in range(2,n+1):
  25.     multiply(F,M)
  26.  
  27. #print(Fib(2))
  28. for i in range(10000):
  29.   if(Fib(i) > 4000000 ):
  30.     k = i
  31.     print(i)
  32.     print(Fib(i))
  33.     break
  34.  
  35. EvenSum = 0
  36. last_numb = i
  37. for i in range(2,last_numb,3):
  38.   if(Fib(i)<4000000):
  39.     EvenSum = EvenSum + Fib(i)
  40.   elif(Fib(i)>4000000):
  41.     break
  42.  
  43. print(EvenSum)
Add Comment
Please, Sign In to add comment