|
|
A189366
|
|
a(n) = n + [n*r/t] + [n*s/t]; r=1, s=sqrt(2), t=(1+sqrt(5))/2.
|
|
3
|
|
|
1, 4, 6, 9, 12, 14, 17, 18, 21, 24, 26, 29, 32, 34, 37, 38, 41, 44, 46, 49, 51, 54, 57, 58, 61, 64, 66, 69, 71, 74, 77, 78, 81, 84, 86, 89, 91, 94, 97, 98, 101, 103, 106, 109, 111, 114, 117, 118, 121, 123, 126, 129, 131, 134, 136, 138, 141, 143, 146, 149, 151, 154, 156, 158, 161, 163, 166, 169, 171, 174, 176, 178, 181, 183, 186, 188, 191, 194, 196, 198, 201, 203, 206, 208
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
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/GoldenRatio]+Floor[(n Sqrt[2])/GoldenRatio], {n, 90}] (* Harvey P. Dale, Apr 12 2013 *)
|
|
PROG
|
(PARI) for(n=1, 100, print1(n + floor(2*n/(1+sqrt(5))) + floor(2*n*sqrt(2)/(1+sqrt(5))), ", ")) \\ G. C. Greubel, Apr 20 2018
(Magma) [n + Floor(2*n/(1+Sqrt(5))) + Floor(2*n*Sqrt(2)/(1+Sqrt(5))): n in [1..100]]; // G. C. Greubel, Apr 20 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|