A189468 - OEIS

%I #11 Sep 08 2022 08:45:56

%S 2,5,7,11,13,16,18,22,24,27,29,33,35,38,40,44,46,49,51,55,57,60,64,66,

%T 69,71,75,77,80,82,86,88,91,93,97,99,102,104,108,110,113,115,119,121,

%U 124,128,130,133,135,139,141,144,146,150,152,155,157,161,163,166,168,172,174,177,179,183,186,188,192,194,197,199,203,205,208,210,214,216,219,221,225,227,230,232

%N a(n) = n + [n*r/t] + [n*s/t]; r=1, s=1+sqrt(2), t=1/2+sqrt(2).

%C See A189459.

%H Ivan Panchenko, <a href="/A189468/b189468.txt">Table of n, a(n) for n = 1..10000</a>

%t r=1; s=1+2^(1/2); t=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}]  (* A189459 *)

%t Table[g[n], {n, 1, 120}]  (* A189467 *)

%t Table[h[n], {n, 1, 120}]  (* A189468 *)

%o (PARI) for(n=1,100, print1(n + floor(2*n/(1 + 2*sqrt(2))) + floor(2*n*(1+sqrt(2))/(1+2*sqrt(2))), ", ")) \\ _G. C. Greubel_, Apr 13 2018

%o (Magma) [n + Floor(2*n/(1 + 2*Sqrt(2))) + Floor(2*n*(1+Sqrt(2))/(1 + 2*Sqrt(2))) : n in [1..100]]; // _G. C. Greubel_, Apr 13 2018

%Y Cf. A189459, A189467.

%K nonn

%O 1,1

%A _Clark Kimberling_, Apr 22 2011