%I #21 Sep 08 2022 08:45:58
%S 9,18,28,37,47,56,66,75,85,94,104,113,123,132,142,151,161,170,179,189,
%T 198,208,217,227,236,246,255,265,274,284,293,303,312,322,331,340,350,
%U 359,369,378,388,397,407,416,426,435,445,454,464,473,483,492,502
%N Beatty sequence for 5+2*sqrt(5); complement of A194151.
%H Vincenzo Librandi, <a href="/A194152/b194152.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = 5*n + A022839(2*n). - _R. J. Mathar_, Aug 25 2011
%F a(n) = floor(n*(5+2*sqrt(5))). - _Vincenzo Librandi_, Oct 25 2011
%t r=Sqrt[5]/2;
%t c[k_]:=Floor[k*r];
%t Table[c[k],{k,1,90}] (* A194151 *)
%t s=5+2*Sqrt[5];
%t d[k_]:=Floor[k*s];
%t Table[d[k],{k,1,90}] (* A194152 *)
%o (Magma) [Floor(n*(5+2*Sqrt(5))): n in [1..80]]; // _Vincenzo Librandi_, Oct 25 2011
%o (Python)
%o from sympy import integer_nthroot
%o def A194152(n): return 5*n+integer_nthroot(20*n**2,2)[0] # _Chai Wah Wu_, Mar 16 2021
%Y Cf. A022839, A194151.
%K nonn
%O 1,1
%A _Clark Kimberling_, Aug 17 2011
|