login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A057045
Let R(i,j) be the rectangle with antidiagonals 1; 2,3; 4,5,6; ...; the n-th Lucas number is in antidiagonal a(n).
4
2, 1, 2, 3, 4, 5, 6, 8, 10, 12, 16, 20, 25, 32, 41, 52, 66, 85, 107, 137, 174, 221, 281, 358, 455, 579, 737, 937, 1192, 1516, 1929, 2454, 3121, 3970, 5050, 6424, 8171, 10394, 13221, 16818, 21393, 27212
OFFSET
1,1
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