phi = (1 + 5**0.5) / 2
Fib = []
for i in range(3):
    Fib.append(int(round((phi**i - (1-phi)**i) / 5**0.5)))
for j in range(3,10000):
    Fib.append(Fib[j-1]+Fib[j-2])

import math
sieve = [True] * 1300000
sieve[0] = False #To mhden kai to 1 den einai prwtoi
sieve[1] = False
for i in range(2,int(math.sqrt(130000)+1)):
    deikths = i * 2
    while deikths < 1300000:
        sieve[deikths] = False
        deikths += i
primes = []
for i in range(1300000):
    if sieve[i] == True:
        primes.append(i)
k =[]
for i in range(9999):
    if Fib[i+1] % primes[i] == 0:
        k.append(i+1)