%I #9 Sep 08 2022 08:45:56
%S 13,27,41,55,69,83,98,112,126,139,153,167,181,196,210,224,238,252,266,
%T 279,294,308,322,336,350,364,378,393,407,420,434,448,462,476,491,505,
%U 519,533,547,560,574,589,603,617,631,645,659,674,688,701,715,729,743,757,772,786,800,814,828,841,855,870,884,898,912
%N a(n) = n + [n*r/s] + [n*t/s]; r=1, s=sqrt(3)-sqrt(2), t=sqrt(3)+sqrt(2).
%C See A189527.
%H G. C. Greubel, <a href="/A189528/b189528.txt">Table of n, a(n) for n = 1..10000</a>
%t r=1; s=3^(1/2)-2^(1/2); t=3^(1/2)+2^(1/2);
%t f[n_] := n + Floor[n*s/r] + Floor[n*t/r];
%t g[n_] := n + Floor[n*r/s] + Floor[n*t/s];
%t h[n_] := n + Floor[n*r/t] + Floor[n*s/t];
%t Table[f[n], {n, 1, 120}] (* A189527 *)
%t Table[g[n], {n, 1, 120}] (* A189528 *)
%t Table[h[n], {n, 1, 120}] (* A189529 *)
%o (PARI) for(n=1,100, print1(n + floor(n/(sqrt(3) - sqrt(2))) + floor(n*(sqrt(3)+sqrt(2))/(sqrt(3)-sqrt(2))), ", ")) \\ _G. C. Greubel_, Apr 20 2018
%o (Magma) [n + Floor(n/(Sqrt(3)-Sqrt(2))) + Floor(n*(Sqrt(3) +Sqrt(2))/( Sqrt(3) -Sqrt(2))): n in [1..100]]; // _G. C. Greubel_, Apr 20 2018
%Y Cf. A189527, A189529.
%K nonn
%O 1,1
%A _Clark Kimberling_, Apr 23 2011
|