A189525 a(n) = n + [n*r/s] + [n*t/s]; r=2, s=1+sqrt(2), t=1+sqrt(3).

2, 5, 8, 11, 14, 16, 19, 23, 26, 29, 32, 34, 37, 40, 43, 47, 50, 52, 55, 58, 61, 64, 68, 70, 73, 76, 79, 82, 85, 87, 91, 94, 97, 100, 102, 105, 108, 112, 115, 118, 120, 123, 126, 129, 132, 136, 138, 141, 144, 147, 150, 153, 155, 159, 162, 165, 168, 171, 173, 176, 180, 183, 186, 189, 191, 194, 197, 200, 204, 206, 209, 212, 215, 218, 221, 224, 227, 230, 233, 236, 239, 241, 244

Offset

1,1

Comments

See A189524.

Links

G. C. Greubel, Table of n, a(n) for n = 1..10000

Mathematica

r=2; s=1+2^(1/2); t=1+3^(1/2);
f[n_] := n + Floor[n*s/r] + Floor[n*t/r];
g[n_] := n + Floor[n*r/s] + Floor[n*t/s];
h[n_] := n + Floor[n*r/t] + Floor[n*s/t];
Table[f[n], {n, 1, 120}]  (* A189524 *)
Table[g[n], {n, 1, 120}]  (* A189525 *)
Table[h[n], {n, 1, 120}]  (* A189526 *)

Prog

(PARI) for(n=1, 100, print1(n + floor(2*n/(sqrt(2)+1)) + floor(n*(1+sqrt(3))/(1+sqrt(2))), ", ")) \\ G. C. Greubel, Apr 20 2018
(Magma) [n + Floor(2*n/(Sqrt(2)+1)) + Floor(n*(1+Sqrt(3))/(1+Sqrt(2))): n in [1..100]]; // G. C. Greubel, Apr 20 2018

Crossrefs

Cf. A189524, A189526.

Sequence in context: A059560 A022842 A356057 * A189369 A055048 A045930

Adjacent sequences: A189522 A189523 A189524 * A189526 A189527 A189528

Keyword

nonn

Author

Clark Kimberling, Apr 23 2011

Status

approved