OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..1000
FORMULA
Round(sqrt(2*A000032(n-1))). - Vladeta Jovovic, Jun 14 2003
PROG
(Python)
from gmpy2 import isqrt_rem, lucas
def A057045(n):
i, j = isqrt_rem(2*lucas(n-1))
return int(i + int(4*(j-i) >= 1)) # Chai Wah Wu, Aug 16 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jul 30 2000
STATUS
approved