OFFSET
1,1
COMMENTS
Complement of A138235; a(n) = A138236(A138235(n)) and A138236(a(n)) = A138235(n). - Reinhard Zumkeller, Mar 07 2008
Numbers k such that A248515(k+1) = A248515(k) + 1 = 1 + least number h such that 1 - h*sin(1/h) < 1/n^2. The difference sequence of A248515 is (0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, ...), so that A138235 = (1, 3, 5, 6, 8, ...) and A022840 = (2, 4, 7, 9, 12, 14, ...). - Clark Kimberling, Jun 16 2015
LINKS
Iain Fox, Table of n, a(n) for n = 1..10000 (first 1000 terms from Reinhard Zumkeller)
Clark Kimberling, Beatty sequences and trigonometric functions, Integers 16 (2016), #A15.
Eric Weisstein's World of Mathematics, Beatty Sequence
MATHEMATICA
Table[Floor[n Sqrt[6]], {n, 70}] (* Vincenzo Librandi, Jun 17 2015 *)
PROG
(Haskell)
a022840 = floor . (* sqrt 6) . fromIntegral
-- Reinhard Zumkeller, Sep 14 2014
(Magma) [Floor(n*Sqrt(6)): n in [1..60]]; // Vincenzo Librandi, Jun 17 2015
(PARI) a(n) = floor(n*sqrt(6)) \\ Iain Fox, Nov 20 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved