|
%I #15 Sep 08 2022 08:44:47
%S 0,3,5,8,11,13,16,19,21,24,26,29,32,34,37,40,42,45,48,50,53,56,58,61,
%T 63,66,69,71,74,77,79,82,85,87,90,93,95,98,101,103,106,108,111,114,
%U 116,119,122,124,127,130,132,135,138,140,143,146,148,151,153,156
%N Integer nearest n*x, where x = sqrt(7).
%H Vincenzo Librandi, <a href="/A022850/b022850.txt">Table of n, a(n) for n = 0..10000</a>
%t Table[Round[n*Sqrt[7]], {n,0,60}] (* _G. C. Greubel_, Sep 29 2018 *)
%o (Magma) [Round(n*Sqrt(7)): n in [0..60]]; // _Vincenzo Librandi_, Oct 24 2011
%o (PARI) vector(60, n, n--; round(n*sqrt(7))) \\ _G. C. Greubel_, Sep 29 2018
%o (Python)
%o from math import isqrt
%o def A022850(n): return (m:=isqrt(k:=7*n*n))+int(k-m*(m+1)>=1) # _Chai Wah Wu_, Jul 31 2022
%Y Cf. A022841.
%K nonn,easy
%O 0,2
%A _Clark Kimberling_
|
