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

2, 5, 7, 11, 13, 16, 18, 22, 24, 27, 29, 33, 35, 38, 40, 44, 46, 49, 51, 55, 57, 60, 64, 66, 69, 71, 75, 77, 80, 82, 86, 88, 91, 93, 97, 99, 102, 104, 108, 110, 113, 115, 119, 121, 124, 128, 130, 133, 135, 139, 141, 144, 146, 150, 152, 155, 157, 161, 163, 166, 168, 172, 174, 177, 179, 183, 186, 188, 192, 194, 197, 199, 203, 205, 208, 210, 214, 216, 219, 221, 225, 227, 230, 232

Offset

1,1

Comments

See A189459.

Links

Ivan Panchenko, Table of n, a(n) for n = 1..10000

Mathematica

r=1; s=1+2^(1/2); t=1/2+2^(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}]  (* A189459 *)
Table[g[n], {n, 1, 120}]  (* A189467 *)
Table[h[n], {n, 1, 120}]  (* A189468 *)

Prog

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

Crossrefs

Cf. A189459, A189467.

Sequence in context: A057922 A243047 A113543 * A004134 A191406 A322969

Adjacent sequences: A189465 A189466 A189467 * A189469 A189470 A189471

Keyword

nonn

Author

Clark Kimberling, Apr 22 2011

Status

approved