OFFSET
1,2
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..1000
EXAMPLE
a(6) = 16 because a(5) is 9 and round(sqrt(3)*9) = 16.
MATHEMATICA
With[{c=Sqrt[3]}, NestList[Floor[c #+1/2]&, 1, 45]] (* Harvey P. Dale, May 03 2011 *)
PROG
(Python)
from gmpy2 import isqrt_rem
A134009_list = [1]
for _ in range(1000):
i, j = isqrt_rem(3*A134009_list[-1]**2)
A134009_list.append(int(i+ int(4*(j-i) >= 1))) # Chai Wah Wu, Aug 16 2016
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Ben Paul Thurston, Jan 09 2008
EXTENSIONS
More terms from Harvey P. Dale, May 03 2011
STATUS
approved