|
%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
|
