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

1, 4, 7, 9, 12, 15, 18, 20, 22, 25, 28, 30, 33, 36, 38, 41, 44, 46, 48, 51, 54, 57, 59, 62, 65, 67, 69, 72, 75, 77, 80, 83, 86, 88, 90, 93, 96, 98, 101, 104, 107, 109, 111, 114, 116, 119, 122, 125, 127, 130, 133, 135, 137, 140, 143, 145, 148, 151, 154, 156, 158, 161, 164, 166, 169, 172, 175, 177, 179, 182, 184, 187, 190, 193, 195, 198, 201, 203, 205, 208, 211, 214, 216, 219

Offset

1,2

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(3)+1)) + floor(n*(sqrt(2)+1)/(1+sqrt(3))), ", ")) \\ G. C. Greubel, Apr 20 2018
(Magma) [n + Floor(2*n/(Sqrt(3)+1)) + Floor(n*(1+Sqrt(2))/(1+Sqrt(3))): n in [1..100]]; // G. C. Greubel, Apr 20 2018

Crossrefs

Cf. A189524, A189525.

Sequence in context: A047537 A247985 A190438 * A184589 A062458 A180122

Adjacent sequences: A189523 A189524 A189525 * A189527 A189528 A189529

Keyword

nonn

Author

Clark Kimberling, Apr 23 2011

Status

approved