login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A189371 a(n) = n + [n*r/s] + [n*t/s]; r=1, s=sqrt(2), t=sqrt(5). 3

%I #10 Sep 08 2022 08:45:56

%S 2,6,9,12,15,19,22,25,29,32,35,38,42,45,48,52,55,58,62,65,68,71,75,77,

%T 81,85,88,91,94,98,101,104,108,111,114,117,121,124,127,131,133,137,

%U 140,144,147,150,154,156,160,164,167,170,173,177,179,183,187,190,193,196,200,203,206,210,212,216,219,223,226,229,233,235,239,243,246,249,252,256,258,262,266,268,272

%N a(n) = n + [n*r/s] + [n*t/s]; r=1, s=sqrt(2), t=sqrt(5).

%C See A189370.

%H G. C. Greubel, <a href="/A189371/b189371.txt">Table of n, a(n) for n = 1..10000</a>

%t r = 1; s = Sqrt[2]; t = Sqrt[5];

%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}] (* A189370 *)

%t Table[g[n], {n, 1, 120}] (* A189371 *)

%t Table[h[n], {n, 1, 120}] (* A189372 *)

%o (PARI) for(n=1,100, print1(n + floor(n/sqrt(2)) + floor(n*sqrt(5)/sqrt(2)), ", ")) \\ _G. C. Greubel_, Apr 20 2018

%o (Magma) [n + Floor(n/Sqrt(2)) + Floor(n*Sqrt(5)/Sqrt(2)): n in [1..100]]; // _G. C. Greubel_, Apr 20 2018

%Y Cf. A189370, A189372.

%K nonn

%O 1,1

%A _Clark Kimberling_, Apr 20 2011

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 09:59 EDT 2024. Contains 371268 sequences. (Running on oeis4.)