|
%I #11 Sep 08 2022 08:45:56
%S 1,4,6,9,12,14,17,18,21,24,26,29,32,34,37,38,41,44,46,49,51,54,57,58,
%T 61,64,66,69,71,74,77,78,81,84,86,89,91,94,97,98,101,103,106,109,111,
%U 114,117,118,121,123,126,129,131,134,136,138,141,143,146,149,151,154,156,158,161,163,166,169,171,174,176,178,181,183,186,188,191,194,196,198,201,203,206,208
%N a(n) = n + [n*r/t] + [n*s/t]; r=1, s=sqrt(2), t=(1+sqrt(5))/2.
%C See A189364.
%H G. C. Greubel, <a href="/A189366/b189366.txt">Table of n, a(n) for n = 1..10000</a>
%t r = 1; s = Sqrt[2]; t = (1 + Sqrt[5])/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}] (* A189364 *)
%t Table[g[n], {n, 1, 120}] (* A189365 *)
%t Table[h[n], {n, 1, 120}] (* A189366 *)
%t Table[n+Floor[n/GoldenRatio]+Floor[(n Sqrt[2])/GoldenRatio],{n,90}] (* _Harvey P. Dale_, Apr 12 2013 *)
%o (PARI) for(n=1,100, print1(n + floor(2*n/(1+sqrt(5))) + floor(2*n*sqrt(2)/(1+sqrt(5))), ", ")) \\ _G. C. Greubel_, Apr 20 2018
%o (Magma) [n + Floor(2*n/(1+Sqrt(5))) + Floor(2*n*Sqrt(2)/(1+Sqrt(5))): n in [1..100]]; // _G. C. Greubel_, Apr 20 2018
%Y Cf. A189364, A189365.
%K nonn
%O 1,2
%A _Clark Kimberling_, Apr 20 2011
|
