OFFSET
1,4
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = n/A054494(n). [Corrected by Charles R Greathouse IV, Nov 05 2014]
EXAMPLE
a(10)=2 because 10/1=10 is not a Fibonacci number but 10/2=5 is.
PROG
(PARI) A010056(n)=my(k=n^2); k+=(k+1)<<2; issquare(k) || (n>0 && issquare(k-8))
a(n)=fordiv(n, d, if(A010056(n/d), return(d))) \\ Charles R Greathouse IV, Nov 05 2014
(Python)
from sympy import divisors
from sympy.ntheory.primetest import is_square
def A054495(n): return next(d for d in divisors(n) if is_square(m:=5*(n//d)**2-4) or is_square(m+8)) # Chai Wah Wu, May 06 2024
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Henry Bottomley, Apr 04 2000
EXTENSIONS
a(34), a(55), a(68) corrected by Charles R Greathouse IV, Nov 06 2014
STATUS
approved