A194832 Triangular array (and fractal sequence): row n is the permutation of (1,2,...,n) obtained from the increasing ordering of fractional parts {r}, {2r}, ..., {nr}, where r= -tau = -(1+sqrt(5))/2.
1, 1, 2, 3, 1, 2, 3, 1, 4, 2, 3, 1, 4, 2, 5, 3, 6, 1, 4, 2, 5, 3, 6, 1, 4, 7, 2, 5, 8, 3, 6, 1, 4, 7, 2, 5, 8, 3, 6, 1, 9, 4, 7, 2, 5, 8, 3, 6, 1, 9, 4, 7, 2, 10, 5, 8, 3, 11, 6, 1, 9, 4, 7, 2, 10, 5, 8, 3, 11, 6, 1, 9, 4, 12, 7, 2, 10, 5, 8, 3, 11, 6, 1, 9, 4, 12, 7, 2, 10, 5, 13, 8, 3, 11
Offset: 1
Examples
Fractional parts: {-r}=-0.61..;{-2r}=-0.23..;{-3r}=-0.85..;{-4r}=-0.47..; thus, row 4 is (3,1,4,2) because {-3r} < {-r} < {-4r} < {-2r}. [corrected by _Michel Dekking_, Nov 30 2020] First nine rows: 1 1 2 3 1 2 3 1 4 2 3 1 4 2 5 3 6 1 4 2 5 3 6 1 4 7 2 5 8 3 6 1 4 7 2 5 8 3 6 1 9 4 7 2 5
References
- C. Kimberling, Fractal sequences and interspersions, Ars Combinatoria 45 (1997), 157-168.
Links
- Wikipedia, Fractal sequence
Programs
-
Mathematica
r = -GoldenRatio; t[n_] := Table[FractionalPart[k*r], {k, 1, n}]; f = Flatten[Table[Flatten[(Position[t[n], #1] &) /@ Sort[t[n], Less]], {n, 1, 20}]] (* A194832 *) TableForm[Table[Flatten[(Position[t[n], #1] &) /@ Sort[t[n], Less]], {n, 1, 15}]] row[n_] := Position[f, n]; u = TableForm[Table[row[n], {n, 1, 20}]] g[n_, k_] := Part[row[n], k]; p = Flatten[Table[g[k, n - k + 1], {n, 1, 13}, {k, 1, n}]] (* A194833 *) q[n_] := Position[p, n]; Flatten[Table[q[n], {n, 1, 80}]] (* A194834 *)
Extensions
Table in overview corrected by Georg Fischer, Jul 30 2023
Comments