OFFSET
0,3
COMMENTS
b(n) = a(n)/a(n-1) begins: 1, 3, 4, 1, 5, 1, 2, 4, 9, 10, 1, 12, 7, 2, 15, 16, ...
LINKS
Max Alekseyev, Table of n, a(n) for n = 0..1000
FORMULA
a(n) = A001177(n!)
EXAMPLE
Least k such that 2! divides Fibonacci(k) is 3: Fibonacci(3)=2, so a(2)=3.
Least k such that 3! divides Fibonacci(k) is 12: Fibonacci(12)=144, so a(3)=12.
PROG
(Python)
n = f = c = d = 1 # f = (n-1)!
fc1 = fd1 = 0 # Fib[c-1], Fib[d-1]
fc = fd = 1 # Fib[c], Fib[d]
while 1:
if fc % f:
if c==d:
fd, fd1 = fc, fc1
t = fc*fc
fc, fc1 = (2*fc*fc1+t), (fc1*fc1+t)
else:
fc, fc1 = (fc*(fd1+fd) + fc1*fd), (fc*fd + fc1*fd1)
c += d
#print '.',
else:
print c,
d = c
f *= n
n += 1
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Aug 08 2012
EXTENSIONS
Terms a(17) onward from Max Alekseyev, Jan 30 2014
STATUS
approved