|
%I #13 Sep 08 2022 08:45:56
%S 3,7,10,14,18,22,25,29,33,37,40,44,48,52,55,59,63,67,71,74,78,82,86,
%T 89,93,97,101,104,108,112,116,119,123,127,131,135,138,142,146,150,153,
%U 157,161,165,168,172,176,180,183,187,191,195,198,202,206,210,214,217,221,225,229,232,236,240,244,247,251,255,259,262,266,270
%N a(n) = n + [n*r/s] + [n*t/s]; r=1, s=cos(Pi/5), t=sec(Pi/5).
%C See A190079.
%H G. C. Greubel, <a href="/A190080/b190080.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = n + [n*sec(Pi/5)] + [n*(sec(Pi/5))^2].
%t r=1; s=Cos[Pi/5]; t=Sec[Pi/5];
%t a[n_] := n + Floor[n*s/r] + Floor[n*t/r];
%t b[n_] := n + Floor[n*r/s] + Floor[n*t/s];
%t c[n_] := n + Floor[n*r/t] + Floor[n*s/t];
%t Table[a[n], {n, 1, 120}] (* A190079 *)
%t Table[b[n], {n, 1, 120}] (* A190080 *)
%t Table[c[n], {n, 1, 120}] (* A190081 *)
%o (PARI) for(n=1,100, print1(n + floor(n/cos(Pi/5)) + floor(n/(cos(Pi/5))^2), ", ")) \\ _G. C. Greubel_, Feb 15 2018
%o (Magma) R:= RealField(); [n + Floor(n/Cos(Pi(R)/5)) + Floor(n/(Cos(Pi(R)/5))^2): n in [1..100]]; // _G. C. Greubel_, Feb 15 2018
%Y Cf. A190079, A190081.
%K nonn
%O 1,1
%A _Clark Kimberling_, May 04 2011
|
