|
%I #12 Sep 08 2022 08:45:56
%S 1,2,3,5,6,7,9,10,11,13,14,15,17,18,19,21,22,23,25,26,27,28,30,31,32,
%T 34,35,36,38,39,40,42,43,44,46,47,48,50,51,52,54,55,56,57,59,60,61,63,
%U 64,65,67,68,69,71,72,73,75,76,77,79,80,81,83,84,85,86,88,89,90,92,93,94,96,97,98,100,101,102,104,105
%N a(n) = floor(s*n), where s=1+sqrt(3)-sqrt(2); complement of A187385.
%C A187385 and A187386 are the Beatty sequences based on r=1+sqrt(3)+sqrt(2) and s=1+sqrt(3)-sqrt(2); 1/r+1/s=1.
%H G. C. Greubel, <a href="/A187386/b187386.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Be#Beatty">Index entries for sequences related to Beatty sequences</a>
%F a(n) = floor(s*n), where s=1+sqrt(3)-sqrt(2).
%t k=3; r=1+k^(1/2)+(k-1)^(1/2); s=1+k^(1/2)-(k-1)^(1/2);
%t Table[Floor[r*n],{n,1,80}] (* A187385 *)
%t Table[Floor[s*n],{n,1,80}] (* A187386 *)
%o (PARI) vector(120, n, floor(n*(sqrt(3) - sqrt(2) + 1))) \\ _G. C. Greubel_, Aug 19 2018
%o (Magma) [Floor(n*(Sqrt(3) - Sqrt(2) +1)): n in [1..120]]; // _G. C. Greubel_, Aug 19 2018
%Y Cf. A187385.
%K nonn
%O 1,2
%A _Clark Kimberling_, Mar 09 2011
|
