|
%I #9 Sep 08 2022 08:45:56
%S 1,2,3,5,6,7,9,10,11,14,15,16,18,19,20,22,23,24,26,28,29,30,32,33,34,
%T 36,37,38,40,42,43,45,46,47,49,50,51,53,54,56,58,59,60,61,63,64,65,67,
%U 68,70,72,73,74,76,77,78,80,81,82,85,86,87,89,90,91,92,94,95,96,99,100,101,103,104,105,107,108,109,111,113,114,116,117,118,120,121
%N a(n) = n + [n*r/t] + [n*s/t]; r=1, s=sqrt(3)-sqrt(2), t=sqrt(3)+sqrt(2).
%C See A189527.
%H G. C. Greubel, <a href="/A189529/b189529.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, A189528.
%K nonn
%O 1,2
%A _Clark Kimberling_, Apr 23 2011
|
