|
|
A189365
|
|
a(n) = n + [n*r/s] + [n*t/s]; r=1, s=sqrt(2), t=(1+sqrt(5))/2.
|
|
3
|
|
|
2, 5, 8, 10, 13, 16, 19, 22, 25, 28, 30, 33, 36, 39, 42, 45, 48, 50, 53, 56, 59, 62, 65, 67, 70, 73, 76, 79, 82, 85, 87, 90, 93, 96, 99, 102, 105, 107, 110, 113, 115, 119, 122, 125, 127, 130, 133, 135, 139, 142, 145, 147, 150, 153, 155, 159, 162, 165, 167, 170, 173, 175, 179, 182, 184, 187, 190, 193, 195, 199, 202, 204, 207, 210, 213, 215, 219, 222, 224, 227, 230, 232, 235, 239
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MATHEMATICA
|
r = 1; s = Sqrt[2]; t = (1 + Sqrt[5])/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}] (* A189364 *)
Table[g[n], {n, 1, 120}] (* A189365 *)
Table[h[n], {n, 1, 120}] (* A189366 *)
Table[n+Floor[n/Sqrt[2]]+Floor[(n*GoldenRatio)/Sqrt[2]], {n, 90}] (* Harvey P. Dale, Aug 21 2015 *)
|
|
PROG
|
(PARI) for(n=1, 100, print1(n + floor(n/sqrt(2)) + floor(n*(1+sqrt(5))/(2*sqrt(2))), ", ")) \\ G. C. Greubel, Apr 20 2018
(Magma) [n + Floor(n/Sqrt(2)) + Floor(n*(1+Sqrt(5))/(2*Sqrt(2))): n in [1..100]]; // G. C. Greubel, Apr 20 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|