%I #9 Sep 08 2022 08:45:56
%S 2,5,8,11,14,16,19,23,26,29,32,34,37,40,43,47,50,52,55,58,61,64,68,70,
%T 73,76,79,82,85,87,91,94,97,100,102,105,108,112,115,118,120,123,126,
%U 129,132,136,138,141,144,147,150,153,155,159,162,165,168,171,173,176,180,183,186,189,191,194,197,200,204,206,209,212,215,218,221,224,227,230,233,236,239,241,244
%N a(n) = n + [n*r/s] + [n*t/s]; r=2, s=1+sqrt(2), t=1+sqrt(3).
%C See A189524.
%H G. C. Greubel, <a href="/A189525/b189525.txt">Table of n, a(n) for n = 1..10000</a>
%t r=2; s=1+2^(1/2); t=1+3^(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}] (* A189524 *)
%t Table[g[n], {n, 1, 120}] (* A189525 *)
%t Table[h[n], {n, 1, 120}] (* A189526 *)
%o (PARI) for(n=1,100, print1(n + floor(2*n/(sqrt(2)+1)) + floor(n*(1+sqrt(3))/(1+sqrt(2))), ", ")) \\ _G. C. Greubel_, Apr 20 2018
%o (Magma) [n + Floor(2*n/(Sqrt(2)+1)) + Floor(n*(1+Sqrt(3))/(1+Sqrt(2))): n in [1..100]]; // _G. C. Greubel_, Apr 20 2018
%Y Cf. A189524, A189526.
%K nonn
%O 1,1
%A _Clark Kimberling_, Apr 23 2011
|