login
Nearest integer to n*sqrt(2).
18

%I #41 Feb 10 2023 14:47:26

%S 0,1,3,4,6,7,8,10,11,13,14,16,17,18,20,21,23,24,25,27,28,30,31,33,34,

%T 35,37,38,40,41,42,44,45,47,48,49,51,52,54,55,57,58,59,61,62,64,65,66,

%U 68,69,71,72,74,75,76,78,79,81,82,83,85,86,88,89,91,92,93,95,96

%N Nearest integer to n*sqrt(2).

%C Let R(i,j) be the rectangle with antidiagonals 1; 2,3; 4,5,6; ...; n^2 is in antidiagonal number a(n). Proof: n^2 is in antidiagonal m iff A000217(m-1)< n^2 <=A000217(m), where A000217(m)=m*(m+1)/2. So m = A002024(n^2) = round(n*sqrt(2)) = a(n). - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Mar 07 2003

%C In the rectangle R(i,j), n^2 is the number in row i=A057049(n) and column j=A057050(n), so that for n >= 1, a(n) = -1 + A057049(n) + A057050(n). - _Clark Kimberling_, Jan 31 2011

%C Number of triangular numbers less than n^2. - _Philippe Deléham_, Mar 08 2013

%H Vincenzo Librandi, <a href="/A022846/b022846.txt">Table of n, a(n) for n = 0..10000</a>

%H Clark Kimberling, <a href="https://www.emis.de/journals/INTEGERS/papers/q15/q15.Abstract.html">Beatty sequences and trigonometric functions</a>, Integers 16 (2016), #A15.

%F a(n) = A002024(n^2).

%F a(n+1) - a(n) = 1 or 2. - _Philippe Deléham_, Mar 08 2013

%e n = 4, n^2 = 16; 0, 1, 3, 6, 10, 15 are triangular numbers in interval [0, 16); a(4) = 6. - _Philippe Deléham_, Mar 08 2013

%t Round[Sqrt[2]Range[0,70]] (* _Harvey P. Dale_, Jun 18 2013 *)

%o (PARI) a(n)=round(n*sqrt(2))

%o (Magma) [Round(n*Sqrt(2)): n in [0..60]]; // _Vincenzo Librandi_, Oct 22 2011

%o (Haskell)

%o a022846 = round . (* sqrt 2) . fromIntegral

%o -- _Reinhard Zumkeller_, Mar 03 2014

%o (Python)

%o from math import isqrt

%o def A022846(n): return isqrt(n**2<<3)+1>>1 # _Chai Wah Wu_, Feb 10 2023

%Y Cf. A063957 (complement of this set).

%Y Cf. A214848 (first differences), also A006338.

%K nonn,easy

%O 0,3

%A _Clark Kimberling_