%I #10 Sep 08 2022 08:45:56
%S 2,6,9,12,15,19,22,25,29,32,35,38,42,45,48,52,55,58,62,65,68,71,75,77,
%T 81,85,88,91,94,98,101,104,108,111,114,117,121,124,127,131,133,137,
%U 140,144,147,150,154,156,160,164,167,170,173,177,179,183,187,190,193,196,200,203,206,210,212,216,219,223,226,229,233,235,239,243,246,249,252,256,258,262,266,268,272
%N a(n) = n + [n*r/s] + [n*t/s]; r=1, s=sqrt(2), t=sqrt(5).
%C See A189370.
%H G. C. Greubel, <a href="/A189371/b189371.txt">Table of n, a(n) for n = 1..10000</a>
%t r = 1; s = Sqrt[2]; t = Sqrt[5];
%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}] (* A189370 *)
%t Table[g[n], {n, 1, 120}] (* A189371 *)
%t Table[h[n], {n, 1, 120}] (* A189372 *)
%o (PARI) for(n=1,100, print1(n + floor(n/sqrt(2)) + floor(n*sqrt(5)/sqrt(2)), ", ")) \\ _G. C. Greubel_, Apr 20 2018
%o (Magma) [n + Floor(n/Sqrt(2)) + Floor(n*Sqrt(5)/Sqrt(2)): n in [1..100]]; // _G. C. Greubel_, Apr 20 2018
%Y Cf. A189370, A189372.
%K nonn
%O 1,1
%A _Clark Kimberling_, Apr 20 2011
|