login
A189471
a(n) = n + [n*r/t] + [n*s/t]; r=1, s=1+sqrt(2), t=1+sqrt(3).
3
1, 3, 6, 8, 10, 13, 15, 17, 19, 21, 24, 26, 28, 31, 33, 35, 38, 39, 41, 44, 46, 49, 51, 53, 56, 57, 59, 62, 64, 66, 69, 71, 74, 76, 77, 80, 82, 84, 87, 89, 92, 94, 95, 98, 100, 102, 105, 107, 109, 112, 114, 116, 118, 120, 123, 125, 127, 130, 132, 134, 136, 138, 141, 143, 145, 148, 150, 152, 154, 156, 158, 161, 163, 166, 168, 170, 173, 174, 176, 179, 181, 184, 186, 188
OFFSET
1,2
COMMENTS
See A189469.
LINKS
MATHEMATICA
r=1; s=1+Sqrt[2]; t=1+Sqrt[3];
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}] (* A189469 *)
Table[g[n], {n, 1, 120}] (* A189470 *)
Table[h[n], {n, 1, 120}] (* A189471 *)
PROG
(PARI) for(n=1, 100, print1(n + floor(n/(1+sqrt(3))) + floor(n*(1+sqrt(2))/(1+sqrt(3))), ", ")) \\ G. C. Greubel, Apr 13 2018
(Magma) [n + Floor(n/(1+Sqrt(3))) + Floor(n*(1+Sqrt(2))/(1+Sqrt(3))): n in [1..100]]; // G. C. Greubel, Apr 13 2018
CROSSREFS
Sequence in context: A128420 A099135 A323028 * A169863 A304500 A047282
KEYWORD
nonn
AUTHOR
Clark Kimberling, Apr 22 2011
STATUS
approved