%I #9 Sep 08 2022 08:45:56
%S 5,11,17,23,30,35,41,47,53,60,65,71,77,83,90,95,102,108,113,120,125,
%T 132,138,143,150,155,162,168,174,180,186,192,198,204,210,216,222,228,
%U 234,240,246,252,258,264,270,277,282,288,295,300,307,312,318,325,330,337,342,349,355,360,367,373,379,385,390,397,403,409,415,420
%N a(n) = n + [n*r/s] + [n*t/s]; r=1, s=sin(Pi/8), t=cos(Pi/8).
%C See A189937.
%H G. C. Greubel, <a href="/A189938/b189938.txt">Table of n, a(n) for n = 1..10000</a>
%F A189937: a(n) = n + [n*sin(pi/8)] + [n*cos(pi/8)].
%F A189938: b(n) = n + [n*csc(pi/8)] + [n*cot(pi/8)].
%F A189939: c(n) = n + [n*sec(pi/8)] + [n*tan(pi/8)].
%t r=1; s=Sin[Pi/8]; t=Cos[Pi/8];
%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}] (*A189937*)
%t Table[b[n], {n, 1, 120}] (*A189938*)
%t Table[c[n], {n, 1, 120}] (*A189939*)
%o (PARI) for(n=1,100, print1(n + floor(n/sin(Pi/8)) + floor(n/tan(Pi/8)), ", ")) \\ _G. C. Greubel_, Jan 13 2018
%o (Magma) C<i> := ComplexField(); [n + Floor(n/Sin(Pi(C)/8)) + Floor(n/Tan(Pi(C)/8)): n in [1..100]]; // _G. C. Greubel_, Jan 13 2018
%Y Cf. A189937, A189939, A189926.
%K nonn
%O 1,1
%A _Clark Kimberling_, May 01 2011
|