OFFSET
1,1
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..5000
MATHEMATICA
a=4/3; b=0;
Table[n+Floor[(a*n+b)^(1/2)], {n, 80}]
Table[n-1+Ceiling[(n*n-b)/a], {n, 60}]
PROG
(PARI) for(n=1, 100, print1(n + floor(sqrt(4*n/3)), ", ")) \\ G. C. Greubel, Jul 22 2017
(PARI) a(n) = n + sqrtint(4*n/3); \\ Michel Marcus, Oct 01 2024
(Python)
from math import isqrt
def A184004(n): return n+isqrt((n<<2)//3) # Chai Wah Wu, Oct 01 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Jan 08 2011
STATUS
approved