OFFSET
1,1
FORMULA
a(n) = n + floor(1/2+(2n/7)^(1/2)).
MATHEMATICA
Table[n+Floor[1/2+(2n/7)^(1/2)], {n, 100}]
PROG
(Python)
from math import isqrt
def A183221(n): return n+(isqrt((n<<3)//7)+1>>1) # Chai Wah Wu, Oct 05 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jan 01 2011
STATUS
approved